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Robert W. Bussard is on the 
staff of the Division Office of the 
Nuclear Propulsion Division of the 
Los Alamos Scientific Laboratory. 
He received the B.S. and M.S. de- 
grees in engineering from the 
University of California, Los An- 
geles. He entered rocket work in 
1949, on the Falcon missile pro- 
gram at Hughes Aircraft Com- 
pany. Following this he spent 
several years in high power 
density reactor development on 
the ANP Project at the Oak Ridge 
National Laboratory. Mr. Bussard has been active in the nuclear 
rocket program since its inception and is the author of numerous 
classified papers in this and allied fields. He is currently at 
Princeton University. 


HE fundamentals, problems and potentialities of rocket 

propulsion systems powered by heat exchanger nuclear 
reactors have often been discussed in the technical litera- 
ture (la, b, c).!. In such systems a reactor with solid fuel 
structure is used as an energy source to heat a working fluid 
which is then expelled through an exhaust nozzle, in the con- 
ventional fashion. Fission energy thus supplants chemical 
energy in the rocket motor combustion chamber. This 
straightforward approach is in the present tradition of reac- 
tor development—adapting the atom to the boiler, so to 
speak, rather than the other way around. Unfortunately 
this approach places an artificial upper limit on estimates of 
potential nuclear rocket performance, since all conventional 
heat engines, whether chemical- or nuclear-powered, are 
fundamentally limited by the energy of the molecular rather 
than the nuclear bond. 

The energy available from fission of U* is about 10’ times 
that from chemical reaction of an equal mass of high explo- 
sive or other combustible mixture. In addition, fission 
energy appears principally as kinetic energy of the two fission 
fragme nts, each with an energy of 60 to 100 Mev, correspond- 
ing to a kinetic “temperature” of the order of 10'2 R. The 
energy released as a result of other nuclear processes, such as 
fusion, beta decay and alpha emission, can also give rise to 
very energetic (several Mev) particles. It is clear that nu- 
clear interaction processes offer a tremendous potential advan- 
tage over chemical reactions, which can never yield more than 
a few ev per particle. 


Fusion 


First in interest, and certainly most speculative, is the ap- 
plication of fusion energy to rocket propulsion. Fusion reac- 
tions in deuterium gas yield a variety of low mass particles 
(T, He*, He‘, p and n) with energies of several Mev each. If 
it were possible to expel the products of these fusion reactions 
rearward from a rocket vehicle, the effective specific impulse 
of the “propellant” would be about 3 X 106 lb, sec/Ib,,,? some 
10‘ larger than presently attainable from chemical rocket pro- 
pellants. More likely, the products themselves could not 


Received Feb. 7, 1958. 

' Numbers in parentheses indicate References at end of paper. 

? Subscripts f and m here denote “force’’ and ‘‘mass,’’ respec- 
tively. 
given simply as sec. 
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practically be constrained to move unidirectionally; however 
large masses of nonfusioned gas could be heated to very high 
temperature by collisions with the energetic fusion products, 
resulting in high specific impulse performance even though 
only a small fraction of the fusionable fuel has been ‘‘burned.” 
For example, a gas mixture in which only one D-D fusion reac- 
tion has taken place for each 104 D nuclei present will reach 
a bulk temperature of about 7 X 106 R (neglecting all losses) 
and will give a specific impulse of roughly 3 X 10‘ sec 
when expelled through a proper nozzle configuration. With 
this specific impulse a 100-ton rocket could take off from the 
sarth, land on the moon and return, all under power, with 
the expenditure of only 5 tons of propellant. 

Of course, gases at temperatures of the order of 106 R and 
higher cannot be contained by solid walls, at least at present, — 
even with use of advanced methods of liquid cooling. How 
ever, since gases are highly ionized at these temperatures: it 
is possible, in theory, to contain them within appropriately 
shaped magnetic fields. This approach is under investiga- — 
tion in the current thermonuclear power program (2), but no — 


satisfactory solution to the containment problem has yet | 


been announced. 


It is evident that thermonuclear fusion reactions can yield — 


propellant temperature and performance far beyond that — 
available conventionally today if such reactions can be made > 
to “go” in a controllable manner. It is equally evident that 
speculation on the possible forms (weight, size, performance, 

etc.) of thermonuclear-powered propulsion systems is — 

less until the controlled release of fusion energy is an ac- 


Radioisotope Decay‘ 


Many radioisotopes are found in nature (e.g., Ra, Rd, ete.). 
others result from the fission of uranium or plutonium nuclei — 
(e.g., Sr®, etc.) and still others can be manufactured by neu- | 


tron irradiation of stable elements placed in nuclear reactors. _ 


The prime difficulty in considering the use of fission product — 
isotopes is that they do not appear singly, but are formed as — 
part of a large group of “mixed” fission products. Chemical 
separation processes are necessary to isolate any one desired — 
fission product from the mixture. On the other hand, radio- | 
isotope production by irradiation of a stable element results — 
in the formation of the desired isotope mixed only with the — 
parent target element. No chemical separation plants are | 
required as for the extraction of single fission products; 
however, it is necessary to utilize reactors especially designed ; 
for irradiation usage. Some characteristics of potentially — 
useful radioisotopes of both types are shown in Table 1.8 

The decay of a radioactive isotope is generally accompanied — 
by the release of fairly large amounts of energy carried by beta | 
particles (electrons), gamma photons or alpha particles. — 
Energetic photons are highly penetrating 
through large masses of material in order to lose their energy. | 
Alpha particles and high energy electrons have mass and are 
charged, and lose their energy quickly by ionization proc-_ 
esses in passage through matter. As a consequence gamma ~ 
emitters are not as useful for radioisotope heat sources as are 
alpha and beta emitters. The decay energy for potentially 


3 Adapted by permission from material in ‘Nuclear Rocket 
Propulsion,’ by R. W. Bussard and R. D. DeLauer, McGraw- 
Hill, New York, 1958. 
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useful beta emitters averages about 1 Mev per disintegration, 
while alpha emitters yield about 4 Mev apiece. Each fission 
process yields two radioisotopes and about one excess neutron. 
Thus about 6 Mev of decay energy can be obtained per fis- 
sion if the two fission product isotopes are beta emitters 
yielding 1 Mev each, and the excess neutron is used to pro- 
duce a 4-Mev alpha emitter from some stable parent ele- 
ment, such as bismuth (yields alpha-active Po). Since the 
initial kinetic energy of the two fission fragments is about 
160 Mev, only some 4 per cent of the total energy of the fis- 
sion process can be converted for use as decay energy. Prac- 
tical considerations reduce this to less than 3 per cent for 
most cases (3). 

The most obvious use of decaying radioisotopes is as heat 
sources, replacing fission reactors or chemical combustion 
used to heat a propellant gas to high temperature. One dis- 
advantage of this application is that no control is possible 
over the rate of energy production; thus auxiliary cooling — 
systems (heat dumps) are required to prevent melting or 
vaporization of the source while not in use. Another dis-— 
advantage is the present limited production capability (8) 
for radioisotope sources. This is a result of the low energy 
conversion ratio obtained when using the fission process to 
produce active isotopes. For an assumed conversion effi- 
ciency as high as 1 per cent an installed production reactor 
capacity of the order of 10* Mw would be required to provide 
the heat sources for one large rocket vehicle per month. This 
is about twenty times the present powerplant capacity of the 
United States. Still another disadvantage is that the specific 
power output of almost any of the useful radioisotopes is very 
low by rocket motor standards. It has been shown (1b) 
that a rocket reactor specific power of the order of 0.5 Mw/lb 
or higher is required for satisfactory missile propulsion by 
nuclear reactor heat-exchanger rocket motors. In contrast, 
the best isotope listed in Table 1, Tm!”, yields only 0.01 


Table 1 Characteristics of some radioisotopes as heat 


sources® 


Artificially Produced Radioisotcpes 
Thermal power available 
at initial use* 


Parent Daughter — 
element isotope Kw/lb Kw/10? Mw yr 

Tm'® Tm!” 10 1400 

Cs!33 Cs!34 6.9 550 

TI 0.17 140 


Production conditions: 

* Material weights include both parent and daughter. 
Production by irradiation in large thermal reactors for two 
half-lives. Initial use assumed at end of irradiation period. 


Chemically Separated Fission Products 
Thermal power available 


Production rate at initial use? 


product Mw yr Kw/lb Kw/108 Mw yr 
Mixed 880 
Sr%/y% 25 
Ru!%/Rh10 40 
Cs'87/Ba!37 93 
Ce!44/Pri44 84 


Production conditions: 

Weights include weight of all chemically similar isotopes. 
Fission products produced in large thermal reactors with 
a fuel cycle time of 180 days (assumed). Initial use assumed 
at 90 days following removal from production reactors. 
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Mw/lb at initial use. The use of radioisotopes as heat sources 
for rocket propulsion thus does not appear very attractive. 

The possibility remains of using directly the momentum 
of the energetic decay particles (4). Unfortunately these 
are emitted isotropically, and at least half of the total decay 
energy will be deposited within the vehicle if half is assumed 
to go rearward. This latter assumption implies a thin layer 
of radioisotope “‘painted”’ on the flat base of the vehicle to be 
propelled. For such a configuration it has been shown (5) 
that about 100 Mw must be dissipated by the vehicle for 
every pound of thrust produced by alpha emission. For bet: 
emission the heat dissipation must be about 3500 Mw pe 
pound of thrust. Clearly, such a scheme is impractical. 

In general, it does not appear possible to make practica 


use of radioisotopes for rocket propulsion. =| 
bik 


‘ission Energy 


Many possible applications of fission energy and reacto 
systems have been suggested for propulsion aside from th 
obvious direct-heat-exchanger approach mentioned _ pri 
viously (la, b, ¢). 

The possibility of using fission fragment momentum direct 
suffers from the objections given previously to use of radio 
isotope decay-product momentum. The ratio of kineti: 
energy to momentum is about the same for fission fragment 
as for energetic (several Mev) alpha particles, and the specific 
heat dissipation which must be achieved is about 80 Mw pe 
pound of thrust produced. 

Fission fragments lose their energy by ionization of th: 
material through which they pass while slowing down. I 
fissions can be made to occur in a reactor in such a way that 
the fragment slowing down takes place in the propellant 
gas rather than in structure (i.e., solid fuel elements) it would 
be possible to heat the gas above the limiting temperatures 
of containing structural materials. In principle, this can be 
done if fissionable fuel can be applied in thin films to the outer 
surfaces of fuel elements in a reactor core. If the fuel film 
is thin enough (<10~4 em) very little self-absorption will 
take place, and half of the fission fragments will travel away 
from the fuel element to heat the surrounding gas directly, 
while half will lose their energy in the element structure, 
which thus requires internal cooling. In the limit, if the pro- 
pellant gas is used both as the internal coolant and the ex- 
ternal flowing gas (in series flow),‘ an increase in specific im- 
pulse of about 40 per cent appears possible over that at- 
tainable from the conventional high temperature heat ex- 
changer system. Consideration of the practical reactor de- 
sign problems associated with the exploitation of this pheno- 
menon leads to the general conclusion that the potential 
propellant performance gain is outweighed by increased com- 
plexity, size and weight of the reactor system. 

Another line of attack® is by making most of the fissions 
occur in the gas phase itself in a reflector-moderated “cavity”’ 
reactor (6). Here an intimate mixture of fissionable material 
and diluent (propellant) is fed into a large void space sur- 
rounded by a neutron moderating material such as DO. 
Fissions take place in the mixture within the void, principally 
by thermal neutrons returned from the reflector-moderator, 
and heat the mixture to a temperature limited only by pres- 
sure-temperature-stress limitations of the container. A 
convergent-divergent nozzle must be located at one end of the 
void core to allow the escape of hot gas. An annular nozzle 
is preferable in order to minimize neutron leakage or stream- 
ing from the core. Fig. 1 shows a schematic outline of such 
a device, similar to those proposed as early as 1949 (7). 

In order to achieve superperformance, temperatures of in- 
terest in this system must be much higher than those reached 
in the combustion process in conventional rocket motors. 


‘ Following a suggestion by W. C. Cooley. Private communi- 
cation, 1955. 
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Since rocket motor combustion temperatures are within a 
few thousand degrees of the boiling points of most structural 
materials, it is clear that all superperformance reactors op- 
erating on the principal of interest here can be analyzed as 
gaseous reactors, whether the fissionable fuel and diluent are 
introduced in solid, liquid or gaseous form. 

For a gaseous reactor operated at a temperature suffi- 
ciently high to assure that all core gases are monatomic, a 
simple, approximate relation can be found between various 
parameters describing the reactor and vehicle operating con- 


ditions. Neglecting ionization effects and assuming perfect 
gases this is 


Pe = 


Here P,, is core gas pressure in lb/in.’, p, is critical fuel density 
in |!) /ft,* Z,, is propellant specific impulse in seconds, and W, 
is the weight in pounds of fissionable material expended dur- 
ing operation; Ito: is the system total impulse, and is simply 
the product of rocket motor thrust F in pounds and operat- 
ing time ¢, in seconds. It is certainly desirable to retain as 
many of the unfissioned atoms of fuel as possible to prevent 
their escape from the system with the outflowing diluent gas. 
The factor S in Equation [1] is a measure of this retention 
or weight separation ability, defined here as the ratio of fuel 
mass expelled during operation to that which would have 
been expelled if no separation had taken place. 

The only nuclear requirement on Equation [1] is that the 
fuel density be sufficient to ensure reactor criticality. Fig. 2 
shows the critical fuel density, taken from the work of Safo- 
nov (6), required for cavity reactors within reflector modera- 
tors of neutronically “infinite” thickness. In practice ‘‘in- 
finite’ means 5 to 10 neutron slowing down lengths, so that 
reflector thicknesses of several feet are of interest. 

Consideration of two examples serves to illustrate the prac- 
tical difficulties which confront useful gaseous reactor propul- 
sion systems: 

| Assume no separation takes place, so that S = 1, and 
that criticality can be achieved with p,; = 0.1 lb/ft’. Assume 
further that it is desired to fly a vehicle which requires 100,000 
lb thrust for 200 see with propellant of performance com- 
parable to that from present chemical rockets; hence Jtot = 
2x 10’ lb sec and J,, = 300 sec. For these conditions Equa- 
tion [1] becomes 


= 2.7 


Here a system pressure of 1000 Ib/in.? will result in the loss of 
27,000 lb of fissionable fuel; this is clearly impractical. 
However, for an allowable fuel loss of 300 lb the system pres- 
sure must be 90,000 lb/in.? which is equally impractical. 

2 Approaching the problem from another viewpoint, as- 
sume an allowable fuel expense of 300 lb and a specific im- 
pulse of 3000 sec, some tenfold better than for present con- 
ventional rockets. With the vehicle impulse previously 
postulated, the system pressure is related to the separation 
ratio by 


For a system pressure of 1000 Ib/in.? the factor S must be 
1.11 X 10~%, implying retention of all but one fuel atom in 900 
of those which would normally escape by being swept out with 
the diluent gas. 

Steady-state gaseous reactors thus appear practical for 
rocket propulsion only if separation ratios of the order of 
10-8 can be achieved. Although it is not clear at present how 
this can be done, the enormous potential value of a rocket 
propulsion system capable of producing both very high 
specific impulse and high thrust-to-weight ratio makes it 
imperative that serious consideration be given to the problem 
of gas-phase separation of atomic species. 


Apna. 1958 
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Fig. 2 Critical fuel density in ‘‘cavity’’ reactors (data from 
Reference 6) 


Reactor Systems 


Aside from the more or less exotic and novel schemes pre- 
viously discussed, fission energy can play another role in the 
rocket propulsion systems of the future. For this, interest is 
in the application of conventional nuclear-powered heat en- 
gine equipment for the production of relatively small amounts 
of very hot gas, on a pulsed or steady-state basis. 

Two general classes of such systems are the “thermo- 
mechanical” and ‘electrical’ open cycles. The first of these 
typically makes use of a succession of nuclear-powered heat 
engines to heat a working fluid to successively higher and 
higher temperatures by use of conventional thermodynamic 
cycles. The resulting hot gas is expelled through a nozzle 
to produce the thrust output of the system. Of course, a 
great deal of waste heat must be dumped in the primary reac- 
tor circuit in order to produce a small amount of high tem- 
perature, high specific impulse exhaust gas from the second- 
ary circuit. The cost of dumping this waste heat is reflected 
only in fixed equipment weight, not in increased propellant 
flow rate. Thus it appears possible to achieve high propel- 
lant performance at the expense of low over-all propulsion 
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GAS COMPRESSORS 


system thrust-to-weight ratio. An illustrative example of a 
possible thermomechanical cycle is shown in Fig. 3. Here a 
conventional closed-cycle gas turbine system drives a series 
set of positive displacement gas compressors, each one fed 
from the preceding unit. The propellant is vaporized and 


heated by cooling the compressor bank and is then supplied “a 
to the first stage pumps. The high temperature outlet gas is I, = 
exhausted to space to drive the vehicle. Waste heat in the — =. 


turbine drive circuit is dumped to space by thermal radiation. 
Although not shown, an auxiliary circuit may be needed to 
dump exhaust circuit waste heat in excess of that which can 
be absorbed by the propellant within structural temperature 
limitations of the compressor bank. 

Electrical open cycle systems follow a pattern similar to 
that just described for the thermomechanical cycles. Here, — 


Woe? 
however, two basic methods of utilizing electrical energy are 


evident. First is that of heating gases in the bulk by passing 
a current through a flowing gas stream, causing dissociation 
and ionization, followed by recombination in the expansion 
process. This process is exemplified by the system shown in 
schematic outline in Fig. 4. Here a nuclear-electric generat- 
ing system is used to drive a d-c are maintained within a flow- 
ing, centrifugally stabilized cylinder of liquid propellant. 
Propellant vaporized from the liquid-walled cylinder feeds 
the arc plasma and is exhausted from the plasma core through 
a hole in the cathode. Russian experimental work (8) on such 
ares has produced plasma temperatures as high as 90,000 R, 
capable of yielding exhaust gas specific impulse an order of 
magnitude higher than that attainable today from chemical 
rocket motors. 

The second method of interest (first proposed by H. Oberth 
in 1929) is that of accelerating charged particles, using elec- 
trical energy to produce directed motion of the individual 
particles of propellant gas. This acceleration of ionized gases 
can be accomplished by use of static electric fields or by mov- 
ing magnetic fields. Here, as for the arc system, a nuclear- 
electric generator set can be used to drive the electrostatic 
or electromagnetic accelerator. Waste heat ~inthe ‘reactor 
circuit must be dumped to space by radiation. 
the propellant exhaust velocity from such a device is limited J 
only by the velocity of light. : 

However, power requirements are proportional to the — 
square of the exhaust velocity while thrust is only linearly 
proportional; thus the specific power (per unit thrust) in- 
creases linearly with increasing exhaust velocity. Since the 
weights of power handling components such as _ turbines, 
compressors, radiators, electrical generators, pumps, and other 
heat engine or electrical plant equipment are generally rather 
directly related to their power capacities, it is clear that in- 
creasing exhaust velocity results in roughly linearly decreas- 
ing thrust-to-weight ratio for any of the reactor-powered sys- 
tems discussed earlier. For these or similar systems yielding 
superperformance from the propellant gas, it is easy to show 
that an optimum exhaust velocity exists for each specified 
vehicle burnout velocity. 
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Consider a vehicle of gross weight Wo, composed of fixed 
weight (payload, guidance, etc.) W,, propellant weight WV, 
and propulsion system weight W,. For a constant propellant 
weight flow rate the gross specific impulse of the complete 
vehicle is given by 


Ft 
[2] 


W. .W.atW,t+W, 


A minimum weight system (arbitrarily chosen here as “‘opti- 
mum’’) results from operation at maximum J, for any given 
value of W,. This fixed load is related through the familiar 
mass-ratio equation to the gross weight and propellant weiglit 


by 


In principle, 


The exponent ¢ is the ratio of vehicle burnout velocity in free 
space to propellant exhaust velocity, § = Av/g_J,,. 

For no net energy losses in the equipment of the propellant 
exhaust circuit, the power required to produce the high tem- 
perature exhaust gas is given by 


where P, is in Mw, and 7, is the exhaust gas expansion effi- 
ciency, generally greater than 0.75 for large nozzles operating 
at low ambient pressure. 

By definition, the weight of the complete propulsion system 
is simply related to the propellant exhaust circuit power re- 
quirements by 


where ®, is the propulsion system gross specific weight in 
lb/exhaust-circuit-Mw. The detailed functional form of 
®, depends upon the choice of system type, since different 
components may be used to make up each different system. 
Combining Equations [2 through 5] to eliminate direct 
weight terms, the vehicle gross specific impulse becomes 


2.18 X 10*1,,@ 


I sp 


where AX = 1 + (W./Wo) 
occurs when 


The maximum value of /, 


Note that the optimum propellant specific impulse or ex- 
haust velocity thus depends principally upon the operating 


time and the propulsion system gross specific weight. For 

this condition Equation [6] reduces to A ee 

I, = 107 

= 


= = 2.3 X 104 Px 
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Now, the vehicle thrust-to-weight ratio is found from Equa- 
tions [2 and 8] to be 


P/Wo = 2.3 X 10! [9] 


for the units given previously. Fig. 5 shows this relation 
graphically for several arbitrary values of ®,, for an expan- 
sion efficiency of 7, = 0.8. Note that one gravity accelera- 
tion (i.e., F/Wo = 1) requires a propulsion system specific 
weight of 18.4 lb/exhaust-circuit-Mw for a propellant specific 
impulse of 1000 sec. 

or the thermomechanical gas compression system pre- 
viously discussed, the gross specifie weight is made up of 
terms describing the turbine-compressor set, reactor, radiator, 
and the secondary (exhaust) circuit gas compressors. A 
“reasonable” estimate of a minimum value of ®, can be 
made for this system, assuming operation at a radiator tem- 
perature of 3000 R, reactor exit gas temperature of 4000 R, 
and the use of very lightweight rotating machinery. Similar 
estimates can be made for the d-c are and ion accelerator elec- 
trical systems, using data on presently available lightweight 
electrical generators. Typical estimated minimum values 
for the principal components of each system are shown in 
Tuble 2, assuming operation at the conditions given above. 
In general it appears that the minimum system specific weight 
attainable with present-day equipment is the order of 1000 
lb ‘exhaust-circuit-Mw. For this value Fig. 5 shows that 
the vehicle thrust-to-weight ratio will be less than 0.07 for 
propellant specific impulse greater than 250 sec. Such sys- 
tems thus appear to be of potential use only in free-fall condi- 
tions where no artificial acceleration requirements (such as 
overcoming the earth’s field) exist. 

A reduction in system specific weight of about a factor of 
100 seems necessary in order to permit operation with high 
specific impulse and high thrust-to-weight ratio. Unfortu- 
nately the outlook is not bright for order of magnitude reduc- 
tions in the specific weight of heat engine rotating machinery. 
Great decreases in weight require corresponding increases in 
the strength of materials used in construction. These do 
not appear on the research horizon at present. Significant 
reduction in the weight of electrical generators may possibly 
be made if electrostatic fields can be used as the basis for 
generator design (9), rather than electromagnetic fields with 
their attendant massive magnetic flux guides. However, 
reduction of the generator weight in electrical cycle systems 
simply reasserts the problem of reducing the turbine-com- 
pressor-radiator weight in the reactor circuit. What is really 
needed here is a conceptually new, and lightweight, method 
of producing shaft power or electrical power from fission. The 
efficient production of electricity directly from nuclear proc- 
esses (10) would at last provide the key to space travel and 
the practical exploration of our solar system. 
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Fig. 5 Generalized performance of reactor system powered 
vehicles 


Table 2 Estimated minimum specific weights for reactor 
system components 


Thermomechanical Gas Compression System 
Specific weight 


| 

| 

Estimated specific per unit power 


weight per unit in exhaust 
Component power handled circuit 
Turbine-compressor 300 lb/Mw = 300 lb/Mw 
Radiator 50 lb/Mw, 20% = 200 lb/Mw 
reactor circuit 
efficiency 
Reactor 5 lb/Mw, 20% = 25 lb/Mw 


reactor circuit 

efficiency 

100 lb/Mw, 5 = 500 lb/Mw 
compressors as- 

sumed in series 


Gas compressors 


Electric d-c Arc or Ion Acceleration System 
Reactor circuit (see above) = 525 lb/Mw 
| Electrical generator 1000 lb/Mw = 1000 lb/Mw 
| Are electrodes and 
container, or ac- 25 lb/Mw = 25 lb/Mw 
celerator structure 
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This paper reviews some of the effects of the presence of 
free electrons in a hypersonic shock layer. Shock front 
structure and boundary layer structure are briefly dis- 

cussed and some of the differential plasma peopertios s are 
mentioned. 
Introduction 
HEN an object travels through air at altitudes below 
something like 200,000 ft and at Mach numbers greater 
than 12 to 15, the production of electrons in the shocked layer 
of air which surrounds the object has a significant effect on 
some of the properties of the shocked air. Abovesome altitude 
like 200,000 ft the influence of the electrons is felt at even 
lower Mach numbers, particularly in the ionosphere where the 
ambient electron density is no longer negligible. 

The influence of electrons is felt not only on aerodynamic 
quantities such as heat transfer, drag and flow field, but also 
on physical quantities such as transport properties, radiative 
emission and absorption, and electromagnetic signal interac- 
tion. The main thesis of the present paper is to show how, 
and under what conditions, the electrons produced by the 
passage of a high speed missile through the atmosphere sig- 
nificantly affect both aerodynamic and physical quantities. 

A plasma may be defined as a gas containing electrons in 
sufficient quantity to seriously affect the physical properties 
of the gas. Plasma physics then is simply the physics of such 
a gas. For example, the air in the shock layer surrounding a 
hypersonic missile may be in a plasma state if the missile is 
traveling at an extremely high velocity. 

We will first discuss in general terms the structure of a hy- 
personic axially symmetric flow field. This will be followed 
by a more detailed discussion of the shock front and boundary 
lay er with emphasis on the electron distribution. We will 
then show an example of an analysis of the Boltzmann equa- 
tion and the relationship of the electron scattering cross sec- 
tion and the elec ‘tron velocity distribution function. 


At hypersonic velocities a detached shock surrounds a blunt 
nosed object as it travels through the atmosphere, provided 
the ambient density is sufficiently high. For Mach numbers 
up to about 30 this density corresponds to altitudes below 
about 200,000 ft for hypersonic missiles. A schematic dia- 
gram of the flow of air about a hemisphere cylinder is shown 
in Fig. 1 

The shock front is a nonequilibrium region in which energy 
is being transferred from the translational degrees of freedom 
of the air molecules to various internal degrees of freedom. 
A more detailed discussion of the shock front is given below. 

Immediately adjacent to the surface of the missile is the 
boundary layer. It is a narrow region in which the tangential 
gas velocity increases from zero at the surface to the local flow 
velocity in the shock layer. 

In between the shock front and the boundary layer is the 
shock layer. At the lower altitudes it is a region of ther- 
modynamic equilibrium; but at higher altitudes it can 


be a region of nonequilibrium due to the increase in relaxation 
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times with decreasing density. In addition, at the higher alti- 
tudes the shock front and boundary layer thicknesses in- 
crease; and at some altitude they may overlap in such a way 
that the equilibrium shock layer disappears. 

The wake, or afterflow region, contains the rocket exhaust; 
or, in the case of a ballistic missile, the wake is a region of low 
density, high temperature, turbulent flow. 

The conditions which exist in the shock layer in the vicinity 
of the stagnation point in front of a blunt object traveling 
through the atmosphere at hypersonic speed are approxi- 
mately the same as those which exist behind a normal shock 
front. Thus the one-dimensional shock relations may be 
used to compute the equilibrium conditions in the stagnation 
region. This has been done for a sphere entering the eart):’s 
atmosphere at about 24,000 fps and at an angle of about 20 
deg to the horizontal (1).2. The velocity and Mach number 
trajectories along with the corresponding values of the 
thermodynamic variables in the stagnation region of thie 
shock layer are shown in Fig. 2. The shock compression 7 is 
the ratio of the density behind the shock to that in front of 
the shock, and the shock strength £ is the corresponding pres- 
sure ratio. 

The composition of the air in the stagnation region is shown 
for the trajectory in Fig. 3 along with corresponding values for 
the effective gamma (ratio of enthalpy to internal energy) 
and compressibility factor g = pV/RT. (Note that g is 
equivalent to the average number of particles into which an 
average initial air molecule has been dissociated or ionized.) 

Before considering some of the details of the boundary layer 
and shock front the flow field will be discussed briefly. In 
the first approximation the shock front can be considered as 
negligibly thick compared to the shock layer. This is a good 
approximation at low altitudes. The flow field is then one of 
axial symmetry and varies from subsonic in the stagnation 
region to supersonic after about 40 deg about the hemisphere 
and to hypersonic along the cylinder. An analytic solution to 
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Fig. 1 Flow about a hemisphere cylinder moving through the 
atmosphere at hypersonic speeds 
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2 Velocity, Mach number, shock strength, compression and temperature trajectories for a sphere entering the earth’s atmosphere 
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Fig. 3 Electron density, effective gamma, compressibility and composition trajectories for the sphere (stagnation region values) 
Nore: Electron densities were computed for a slightly different trajectory 
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such a flow field has not been obtained; however, numerous 
attempts have been made to obtain a good numerical solu- 
tion. Most of these attempts have failed, primarily because 
of problems involving numerical instability near the sonic 
line. Even for the case of a constant gamma gas such as 
helium (below about Mach 25) a good numerical solution is 
difficult. 

It is difficult to obtain a numerical solution because both the 
shock front and the sonic line are hard to locate. Also, when 
the shock front is assumed to be a discontinuity, it is difficult 
to integrate across it by ordinary methods. This difficulty can 
be circumvented by smearing the shock front out with the use 
of an artificial viscosity (2). Several mass points can then be 
located within the shock front and the integration carried out. 
The trouble with this method is inadequate machine storage 
space. Recently, a promising method has been developed by 
Harlow at the Los Alamos Scientific Laboratory (3, 4). It is 
called the ‘‘particle-in-the cell’? method. The name follows 
from the representation of the hydrodynamic flow field by dis- 
crete particles that move in fixed cells. 


Shock Front Structure in Air 


We are concerned here primarily with the determination of 
electron densities rather than atom or molecule densities. 
However, since electrons are produced by ionizing collisions in- 
volving one or more of the heavier particles, it becomes neces- 
sary to determine the particle density of the atoms or mole- 
cules which are most directly concerned in the production of 
free electrons. In the case of air the situation is complicated 
by the formation of nitric oxide, NO. At the ambient con- 
ditions which exist in the atmosphere in front of the shock, the 
concentration of NO is completely negligible. However, be- 
hind the normal shock front the concentration may be as high 
as 10 per cent (cf. Fig. 4). The rate at which NO is formed in 
a shock front is difficult to calculate because the equilibrium 
concentration of NO is a maximum at about 5000 K and drops 
off rapidly for either higher or lower temperatures. In addi- 
tion, there is a strong density dependence. Also plotted on 
Fig. 4 are the equilibrium values for y and g as functions of 
density and temperature. 

Corresponding to the previous definition, the shock front is 
the nonequilibrium region between the undisturbed air in 
front of the shock and the region in the shock layer behind the 
shock where thermodynamic equilibrium can first be assumed 
to exist. The “local” gas temperature inside the shock front 
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Fig. 4 Nitric oxide (NO) concentration, effective gamma and 


compressibility factor as functions of density and temperature te = + At [3] 
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Fig. 5 Schematic diagram showing temperature distribution 
through a shock front in air 


is considerably higher than the equilibrium temperature be 
hind the shock front. A schematic diagram of the temper:- 
ture distribution is shown in Fig. 5. 

The ionization potential of NO is lower than any oth: 
major component of air in the temperature range of interest. 
Hence, most of the electrons present in the shock layer pro!) 
ably come from reactions involving NO or NO*; for example 


is the most influential reaction in controlling the electron 
density. 

In view of the complexity introduced into the structure of a 
shock front in air by the formation of NO, the structure of 
shock fronts, first in a monatomic gas, argon, and second in a 
pure diatomic gas, oxygen, will be considered. Let A, be 
the distance between the front of the shock and the point 
behind the shock at which the electron density can be assumed 
to have (essentially) reached its equilibrium value. An esti- 
mate of A, for argon, oxygen and air has been made for normal 
shocks as a function of density in front of the shock. The 
electron and plasma temperatures will be assumed to be 
equal for the three gases. 

Consider a shock wave in pure argon. If an argon atom is 
initially at density and temperature p,; and 7, where the 
subscript 1 refers to conditions in front of the shock, then 
upon being struck by a strong shock the translational de- 
grees of freedom are excited to conditions p’ and 7” within a 
distance of about two mean free paths 2Z;. If 7’ is sufficiently 
high, energy will be transferred to other degrees of freedom 
such as electron excitation and ionization with resulting de- 
crease in temperature and increase in density. At the peak 
of the shock front, represented by the conditions p’, 7’, the 
only particles which exist are argon atoms (in the ground 
state); hence the ionization reaction 


or the excitation reactions A + A + Aion \ 


A+A*2e+At+AJ 


are the only possible electron producing reactions. Argon 
atoms in their ground states are represented by A, excited 
states by A*, and singly ionized states by A+. The two-step 
process probably requires about the same order of time as the 
one-step process; hence, it is sufficient to assume that the 
latter is the predominant process initially. 

The production coefficient is denoted by 8, and-the recom- 
bination coefficient by a. The coefficient 8; has been evalu- 
ated by assuming Maxwellian velocity distributions (5). 

Once a sufficient number of electrons has been produced by 
reaction [2] the following reaction takes over because of the 

_ much larger cross section 
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The point where this reaction becomes predominant is 
called the onset of ionization; after this point is reached elec- 
tions are produced at a much higher rate. In fact, the time 
required for the onset of ionization to occur probably repre- 
sents over half of the shock front thickness; hence, it will be 
assumed that A, is given to a sufficient degree of accuracy by 
the onset thickness. 

If it is assumed that reactions [2 and 3] are the two pre- 
dominant reactions up to the onset position, then the source 
term may be written as 


O( une) 
ox 


dn. on. 


dt a 


In this relation n, is the electron density, wu the material 
velocity, X the Eulerian distance coordinate and 


where (na is the atom density.) Recombination terms as well 
as other processes should be included in the expression for S; 
however, it is a good approximation to omit these terms in the 
region including and preceding the onset point. 

In a Lagrangian system, where 2’ is the distance coordinate, 
Equation [5] becomes 


dne On, p ou rs ‘las 
dt ot x’ 1 Or 
If the equation of continuity 
OV Ou 
V = 
ot or p 
is inserted, the result is 
din(n,V) 


dt Ne 


Let the degree of ionization x be given by 


In this relation it is assumed that the electron and ion densi- 

ties are equal. Since p = m(na + n,) where m is the mass of an 

argon atom, then n,V = x/m, and Equation [6] becomes 
dx Bip 


+ (8, — 26,) — (8, B)......18) 
t m m m 


Up to the onset point the density p is fairly constant and 
the degree of ionization extremely small; hence Equation [8] 
may be integrated to give 


At the onset point the first two terms on the right of [8] are 
equal. From this and the assumption that the density does 
not change much up to the onset point, the (Eulerian) onset 
distance is given by 

29 X 10-4 
From this relation it is apparent that the shock front thick- 
ness goes, at least to a first approximation, inversely as the 
ambient density in front of the shock. The coefficients 6, and 
8» are, of course, functions of the shock velocity u. 

The production coefficients 6, and 82 depend on the tem- 
perature and not on the density. The dependence of {; is on 
the atom, or plasma, temperature 7’, whereas 6, depends on 
the electron temperature 7. If the electron temperature is 
about the same as the plasma temperature, then 6. >> 6, and 
X” becomes 


2.9 X U 
x’ = 
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A few approximate values of XY” are given in Table 1. 


Table 1 Approximate values for the thickness of an 


argon shock front (assuming T, = T) 


X"(p. = 1.3 X = 1:3 

U, cm/sec 10-* gm/cm*) 10-7 gm/cm?) 
5 X 108 3.4 X 107-4 cm 3 cm 
6 X 108 1.6 X 10-4cm 2 cm 


Table 1 shows that if the density is decreased by a factor of a 
10~4 the shock front thickness becomes of the order of cen- — 
timeters. If these values are extrapolated to the case of 
air, then at an altitude of about 200,000 ft the shock front 
thickness may be of the same magnitude as the shock — 
layer thickness. 

For argon the effective electron temperature inside the — 
shock front is different from the plasma temperature (5). In | 
order to estimate the electron temperature an additional rela- — 
tion is needed. This may be taken as an energy balance be- | 
tween ionizing collisions and Coulomb collisions. If this is — 
done it is found that the shock front thickness for argon is — 
considerably larger than the values noted in the table. How- 
ever, for air, correspondence principle arguments may be used | 
to show that the electron and plasma temperatures =i 
probably about the same, although it should be pointed out 
that the concept of a “local” temperature inside the shock | 
front is highly questionable. 

A shock front in pure oxygen is considerably more difficult 
to analyze than one in pure argon. In the first place there are 
additional degrees of freedom of rotation, vibration, dissocia- — 
tion and molecular electronic excitation and ionization. In — 
addition, both the oxygen molecule and, particularly, the — 
oxygen atom have large electron affinities. Processes in — 
which these various quantities are concerned involve cross | 
sections which, for oxygen, are largely unknown. 

It is obvious that the equilibration process occurs in such a — 
way that the excitation of the various degrees of freedom — 
occurs conjointly; however, the rates of excitation differ. 
In order to approximate the equilibrium process we may as-— 
sume that vibrational equilibrium is achieved in times small — 
compared to those for dissociation equilibration. The — 
corresponding gamma used to compute the thermodynamic 
quantities is 


For this value of y 


and (denoting onset values by a double prime) 


2 1 1 
+ 1) no 
The corresponding compression is given by 
- Na =F 1 
and the temperature by 
n 


We have used the concept of onset for oxygen because of its 
usefulness in the argon analysis. Actually, the analogy with 
argon is not very good. 

In addition to the internal degrees of freedom which oxygen 
has, there are many low-lying electronic states. Also, impact 
ionization is relatively unimportant in an oxygen shock front. 
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The predominant reactions are probably 


O. + 0O.+~20+0, (AE ~ 5.1 ev) 
O + O02: > 30 (AE ~ 5.lev)}....... [16] 
0+0—-0,+ +e (AE ~ 7.4 ev) 


Ps production coefficients in the expression for the source 


— 


terms S in Equation [5] are simply the Maxwellian averages 
ov; hence, from Equation [11] 


This relation may be used to obtain an order-of-magnitude 
estimate for the oxygen shock front thickness. (It would be 
more correct to use the threshold energies which are some- 
what higher than the ionization energies. Also, it is assumed 
in [17] that the difference between the ionization energy and 
the threshold energy may be scaled.) 

Instead of doing a detailed analysis for oxygen we will con- 
sider the more complicated case of air. The complexity in- 
troduced into the computation of the shock front structure 
by the formation of nitric oxide has already been mentioned. 
The presence of NO will, because of its low ionization potential 
(9.5 ev), tend to make the onset distance less. However, 
there are processes other than impact ionization from which 
electrons may be obtained. For example, reaction [1] can 
easily go to the right but NO* is likely to dissociate before 
colliding with an electron; hence the recombination process 
[1] going to the left can be neglected prior to onset. Some 
other possible electron-producing reactions might be 


In addition to the displacement type of reactions [1, 18 and 
19], there are various impact ionization reactions. 


my Table 2 Some energy levels for air 
First 
Dissociation Ionization electronic 
potential, potential, level, 

Specie ev ev ev 
No 9.7 15.51 6.0 
Oz 5 5.08 12.5 0.9 
14.54 2.38 
13.61 1.96 
NO 6.48 9.5 3.8 
NOt 10.6 
Or 35.15 
N+ 29.61 
N,* 8.72 
O- —1.48 
A 15.76 11.5 


In order to obtain a feeling for the magnitude of the shock 
front thickness in air it seems reasonable to compare the re- 
laxation times for the two (competing) processes 


Ba 
[20] 
Bs 
N40 
In these reactions represents any particle, 8; is the produc- 
tion coefficient, and 73 the relaxation time for [20]; 64, 7s are 
the corresponding quantities for [21]. If one of the relaxation 


times is much smaller than the other then the process cor- 
responding to the longer time can be neglected. However, if 
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Fig. 6 Details of an argon shock front (V ~ 6 X 10° cm/sec, 
p: ~ 0.001 gm/cm’, 7, ~300 K) 


73 and 7, are about equal, some sort of average relaxation time 

must be used. The energy balance for [20] is simply th 

ionization potential for NO, i.e., AE; = 9.5ev. For [21], the 

energy balance is AF, = 3 ev. 

The ratio of the relaxation times is, very roughly, given by 
rs AEW/kT 


If the temperatures are equal then 


Since kT is about 1 ev or less, it follows that 73 >> 74; hence 
reaction [21] should be the predominant reaction for the de- 
termination of A, for air. 

The analysis of the rate at which [21] proceeds inside a 
shock front in air depends upon the rates at which N and O 
are separately produced. The rate for the production of 
atomic oxygen can, to a good approximation, be computed 
from the coupled reactions 

Bs 
+ > 20 + 


Be 
O + O. > 30 


(The rate for the production of atomic nitrogen can similarly 
be estimated.) 
If x is the degree of dissociation, then for [22] 


dx 5.1 
— = 6.1087"? e~5-VkT — 4 — z)...{% 
6.10'87"/? + ) pin( x)... [23] 


This equation, together with the equation of state and the 
Rankine-Hugoniot relations, can be numerically integrated 
across the shock front. An approximate integration of [23] 
along with the corresponding equations for the production of 
atomic nitrogen and for the production of electrons [21] indi- 
cates that between altitudes of 100,000 and 200,000 ft the 
shock front thickness is of the order of centimeters,’ thus im- 
plying that the concept of an inviscid (or equilibrium) flow 
field above these altitudes may not be correct. 

The above discussion should serve as an introduction to the 
general problem of shock front structure and the distribution 
of electrons through the shock front. However, it is not in- 
tended that a detailed analysis should be presented in this 
paper. Such an analysis has been made for argon and some 


3 Recent experimental results obtained at Princeton University 
by V. Blackman and his collaborators confirm this theoretical 
result (private communication). 
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numerical results are given below. In obtaining the argon re- 


sults both the effect of an electron temperature different from ~~ 700 


the plasma temperature and the influence of various recom- 


bination processes were considered. 7 


The numerical results for argon (see Fig. 6) have received 
some rough experimental confirmation by LaPorte and his 
students at the University of Michigan and by Kantrowitz 
and his collaborators at Cornell. In addition, some measure- 
ments were made at the Los Alamos Scientific Laboratory 
which also seem to confirm the argon results. 

It has been difficult to obtain similar results for air because 
of the following reasons: Analytically the problem is compli- 
cated by the presence of additional degrees of freedom, the in- 
fluence of low-lying excited electronic states, the numerous 
types of chemical] species, and the lack of information on 
cross sections (or reaction rates). Experimentally, it has not 
been possible, until very recently, to produce sufficiently 
strong shocks in shock tubes to give the degree of ionization 
needed to make the same types of measurements as were 
made on the argon shock front, viz., measurements of electri- 
cal conductivity and of radiation. 


Electric Potential in Gaseous Shock Fronts 


With the above review of the flow structure about a hyper- 
sonic vehicle we are now in a position to discuss some of the 
plasma properties of the flow field. One such property is the 
positive electric potential that exists in strong gaseous shock 
fronts. Evidence of this potential was found about ten years 
ago by probe measurements made at the Los Alamos Scientific 
Laboratory. 

These positive potentials come about as follows: In the 
argon shock front, once the onset point has been reached, there 
is a sudden increase in pressure and number density (including 
electron density), and a decrease in temperature. The result- 
ing gradients exert a force on the electrons which is such as to 
cause a diffusion of electrons toward the front of the shock. 
The current may be symbolically written as 


j (total) = 7 (pressure) + 7 (density) — j (temperature) — 
j (Coulomb) 


The first term on the right represents the current resulting 
from the positive pressure gradient, and the second term, that 
resulting from the positive gradient in the number density of 
the electrons. The third term is negative because the tem- 
perature gradient is negative. The fourth term results from 
the electric or Coulomb force between the electrons and the 
ions. The first three terms are shown schematically in Fig. 7. 
It is evident that the current due to the gradient in the elec- 
tron density is by far the most important contribution to the 
total current. 

In obtaining the currents the gas was treated as a binary 
mixture consisting of electrons n; and partially charged posi- 
tive ions nm. where n = m, + m. The equation used was 

Op 


: —ne 
7 = —— Dis a 
Ne ox Ln np oX 
Di» 


Nop 


ST (me — natydX....... [24] 


The evaluation of the binary diffusion coefficient Dj. and 
the thermal diffusion ratio kr involves integrals containing the 
Coulomb force. Such integrals are usually divergent unless 
finite limits are assumed. 

In our case the limit is the Debye shielding distance which 
is simply the relaxation distance which an electron can be ex- 
pected to diffuse ina plasma. It is given by h where 


kT 
ht = — 


and is contained in the electrostatic potential which must 
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Fig. 7 Relative current densities due to electron diffusion in an — 
argon shock front - 
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satisfy Poisson’s equation, where the potential is given by 


o= (A = const) 


The concept of plasma shielding implies that the density of - 
free charges, i.e., ions and electrons, is sufficiently large so that — 
electron collisions are collective rather than individual. This — 
means that the electron sees the effective Coulomb field of a 
number of charges rather than that of each charge sepa- — 
rately. Such a condition exists for electron densities suf-— 
ficiently high to cause onset incidence. 

Equation [24] was numerically evaluated to obtain the 
charge density resulting from the diffusion of electrons toward 
the front of the shock. This diffusion process causes a dipole 
or double layer to form in the shock front. If the charge dis-_ 
tribution is integrated once the electric field is obtained, a 
second integration gives the electric potential. The numerical © 
results are shown in Fig. 8. 

The theoretical values given above for argon are in qualita- — 
tive agreement with the experimental values obtained at Los 
Alamos. The corresponding theoretical values for air have 
not been worked out but experimental values have been ob- 
tained. The measured potential for air is considerably smaller 
than it is for argon, in accordance with the theory of the 
shock front given in the previous section. 

Although the theory for air has been worked out, numerical 
values will be difficult to obtain until better values for some of 
the reaction cross sections are available. These cross sections 
are not as easy to compute for air as they are for argon be- 
cause of the presence of diatomic molecules including nitric 
oxide NO and because of the many low-lying electronic states 
that occur in high temperature air. 


Plasma Effects in Boundary Layers 


The air in the shock layer about a Mach 20 missile may be 
ionized as high as 1 per cent. In this case there may be a 
sufficient number of electrons to diffuse through the boundary 
layer causing a negative potential within the boundary layer. 
Equation [24] should then be applicable. The second term 
on the right, which involves the pressure gradient, is negligi- 
ble but the first and third terms both cause diffusion currents 
to flow toward the wall. (Note that in the shock front these 
two terms have opposite effects.) 

In general, a missile flying at velocities where electron dif- 
fusion might become important will also have a reacting wall 
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Fig. 8 Charge density, electric field and electric potential dis- 
tribution in an argon shock front 


which would submerge the effects of the electron diffusion. 
It is extremely difficult to predict what the wall effects may 
be and, although considerable work has been done on this 
problem in the study of probe potentials, the theory is far 
from being complete. 

Another problem which may involve plasma effects is that 
of transport properties. For example, the thermal conduc 
tivity of air may be expected to increase with temperature up 
to the point where plasma effects must be considered. For 
air, this temperature is about 5000 K and, for nitrogen, some- 
what higher. Theoretical and experimental values for the 
thermal conductivity of air and nitrogen are shown in Fig. 9 


(6). 


Boltzmann Equation Analysis 


It has been implied in a foregoing section that one may 
compute the quantity ov by simply averaging over a Maxwel- 
lian distribution. However, this is in general not true be- 
cause the velocity distribution functions are probably non- 
Maxwellian inside a shock front—especially for the electrons 
(5). In addition, the cross sections are not well known. 

The quantity ov also enters into the dielectric coefficient of 
the shock layer plasma. The dielectric coefficient, in turn, 
affects the transmission of microwaves through the shock 
layer. 

In both the computation of shock front reaction rates and 
microwave transmission the Boltzmann equation is used. 
(Of course, the hydrodynamic flow equations, with or with- 
out dissipation terms, may also be derived from the Boltz- 
mann equation.) In order to show this, a brief discussion 
and application to microwave transmission is given below. 


If f. (v, r,t) is the velocity distribution function for particle _ 


s, then the Boltzmann equation is 


— Of; ofs 


where m, F, is the force on a molecule of mass m, and the 
right-hand side is the collision term which denotes the rate 
at which f, is being altered by collisions. We have followed 
here the notation of Chapman and Cowling (7) so that 


FRONT OF SHOCK- Of 
I 


In case the distribution function is spatially independent, 
and in the presence of an electric field E(t), Equation [25] re- 


of _ 
ot Ol |e 


where the field is directed along the y-axis. If we now make 
the assumption that the plasma is weakly ionized so that the 
collision term contains only interactions between electrons and 
neutral particles, then the transport equation may be further 
reduced. This is done by using a suitable integrating multi- 
plier and applying the equation of change (7). The result is 
ot m 
In deriving this equation it has been assumed that the 
velocity distribution function of the neutral particles is Max- 
wellian. Hence, n is the number density of the neutrals. 
Also 
The total scattering cross section is given by o(v) and is 
equal to the momentum transfer cross section for low energy 
impact. 
In order to solve [27] it is necessary to know both f and 
a(v). In practice the distribution function is usually assumed 


to be Maxwellian. In addition, it is often assumed that 
either o(v) is constant or that ov is constant. Both cases 


have received considerable attention. (It is easy to show that, 
under certain assumptions, Equation [27] may be further re- 
duced to the standard wave equation for which many solu- 
tions are available.) 
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Fig.9 Thermal conductivity of nitrogen and of air as a function 
of temperature 
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Electron Scattering Cross Sections For example, one of the factors implicit in the assumption 


| the of Maxwellian distributions is that the external energy is 

nate In the above analysis it has been shown that the total cross negligible. If this is not the case, as for externally applied 

wed section o(v) plays an intimate role in the evaluation of reac- electric and magnetic fields, then the distribution function for 
tion rates. Also, the cross section for the scattering of elec- the electrons may well be highly non-Maxwellian. : 
trons from atomic or molecular systems is involved in the di- Let us assume we have solved the Boltzmann equation and — 
electric coefficient. In fact, the absorptivity of microwaves in the resulting distribution is a Druyvesteyn function (this is 
plasmas is proportional to some power of the cross section. In the distribution function resulting when the thermal energy — 

lent, ase the electron and plasma temperatures are the same, as of the gas is much less than the externally applied electric 

] re- seems likely for an equilibrium shock layer in air, then the energy). A simple analysis then shows that the electron 
cross section involved is the elastic scattering cross section for density may be several orders of magnitude less than that ob- _ 
electrons with the various atomic species. (This is probably tained from application of the Saha equation (with Maxwel-— 

[26] true for missile velocities less than about Mach 20, although lian distribution functions). 

it is possible that some inelastic processes compete.) 

In recent years there have been numerous attempts to com- ; , 

- pute or to measure the appropriate cross sections. Unfortu- Magnetohydrodynamics 

ad nately, the results of the various efforts are not in agreement. The above three sections have to do with differential plasma _ 

ther In addition, the numerous theoretical evaluations appear to properties. The term magnetohydrodynamics (or hydro-— 

be incomplete or oversimplified. Experimentally, the difficulty magnetics or magnetoaerodynamics) is primarily concerned _ 
+. is due to the lack of good environmental simulation. Theo- with integrated plasma properties. The application lies in 
" retically, the difficulty is due to the ane etad of exchange and the simple fact that the air flowing in the shock layer about 7 . 

m polarization effects at low energies. Exchange effects occur hypersonic missile may be in a plasma state due to thermal 

ia because of the Pauli principle; that is, the spins of each of the agitation or due to artificially induced ionization such as 
electrons involved in a collision can only have two orienta- sodium injection. If this plasma is then acted on by electro-_ 

the tions. Polarization effects occur because at the low energies magnetic fields the plasma properties may be altered in vari- 


AX- concerned the orbital electrons can be perturbed by the ous ways. Recently there have been numerous publications in _ 


foregoing analysis ie f(v) is Maxwellian. This is a good as- a Bond, J. W., Jr * Equilibrium Conditions Behind a Nor- 

ind sumption for the heavy particles, i.e., the atoms and ions, but the Averett Comp, 

’ MSD 1456, Dec. 1955. 

ed it is probably a poor assumption for the most important 2 Kolsky G.. “A for tha 

nat plasma constituent, the electrons, in the missile velocity range Transient Hydrodynamic Shock Problems in Two Space Dimen- _ 

ses with which we are concerned. This is because of the following sions,’’ Los Alamos Scientific Laboratory LA-1867, Sept. 1954. 

at, reason: At low temperatures the electron density is low and 3 Harlow, F. H., ‘A Machine Calculation Method for Hy-— 

re- atom-atom, atom-molecule, or molecule-molecule collisions drodynamic Problems,’ Los Alamos Scientific Laboratory, 

lu- predominate. In these collisions the masses involved are al- LAMS-1956, Sept. 1955. oa 
most equal and energy transfer is large, thus providing a _ 4 Harlow, F. H., “The Particle-in-Cell Method for Two- 
rapid means of temperature equilibration. At very high tem- me Se ns P — Los Alamos Scientific 
peratures the electron density is large and electron-electron 1 5. Bond, J. W., Jr., Physical Sinden, vol. 105, 1957, p. 1683. 
collisions predominate. Again, because of the equal masses 6 Finkelnburg, W., and Maecker, H., H. d. Physik, Bd. 
involved the temperature equilibration is rapid. However, for XXII, Springer-Verlag, Berlin, 1956. 
intermediate temperatures where collisions between electrons 7 Chapman, S., and Cowling, T. G., “The Mathematical 
and the various heavier particles are significant the tem- Theory of Non-Uniform Gases,’’ Cambridge, The University 
perature equilibration may be slow. Press, 1952. 

"] In most of the equation of state calculations that have been _ 8 Patrick, R. M., “Magnetohydrodynamics of Compressible 
applied to air in the temperature range of interest to us, a Fluids,” Ph.D. thesis, 1956, Cornell University, Ithaca, N. Y. 
modified Saha equation has been used to compute the equi- 9 Rossow, V. J., “On Flow of Electrical Conducting Fluids 
librium constant. This is simply the mass transfer law and overs Fiat Plate in the Presence of a Transverse Magnetic 

Field,” NACA TN 3971, May 1957. 
assumes that the distribution functions are Maxwellian. If 10. Pai, S. L., “Magnetohydrodynamics and Magnetogasdy- 
this is not the case, as seems likely for the electrons over a namics,” University of Maryland, TN BN-59, Sept. 1955. 

= wide energy range, then the electron densities could be in 11 Sears, W. R., “Recent Developments in Magnetoaerody- & 

error by a considerable factor. rye” namics,’’ IAS Summer Meeting, Los Angeles, June 1957. — 
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presence of the incoming electron. 


Velocity Distribution Functions 


Another differential plasma property is the velocity dis- 
It has been assumed in most of the 


this field, particularly in the Russian literature. Reference is — 
here made to a few of the articles available in the American — 


literature (8, 9, 10, 11). 
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Equations were developed relating solid propellant grain 
geometry to the important ballistic parameters of cross- 
sectional loading density, sliver or tail-off fraction, pro- 
gressivity ratio, initial surface and web. The equations 
were then solved by means of high speed computers and the 
results graphed in a readily usable form. If the loading 
fraction, sliver, progressivity ratio and web are given, 
grain designs which meet or most nearly meet given 
requirements can be found from the graphs. The graphs 
are useful in determining the effect of slight variations — 


in design parameters. Illustrations of the use of the _ 


graphs are given for the internal-burning star and the — 
internal-burning wagon wheel configurations. Calcula- 
tions have also been made on a modified form of the 
wagon wheel? which permits one, two or three levels of 
thrust during burning, 
Introduction 


NTERIOR ballistic design is at times severely hampered 

by limiting boundary conditions. Generally, design is 
based on a set of performance criteria such as payload, burnt 
velocity, total weight, burning time and motor dimensions. 
Within these criteria, the designer is free to specify pro- 
pellant, method of retaining the propellant charge in the 
motor, and charge configuration. Maximum performance 
depends on an optimum combination of these three vari- 
ables. 

An examination of the literature indicated that the oppor- 
tunity to optimize charge design has not always been used to 
advantage, primarily because the designer has had no prac- 
tical way to find the optimum geometry. Instead, charges 
have been designed empirically with compass and straight- 
edge until workable designs were obtained; in many cases 
designs were not optimum. Such empirical technique is 
not in agreement with the increased demands for more effec- 
tive ballistic systems. 

In recent years several articles ‘have been published on 
the mathematics of grain design (1, 2, 3). Indeed, the con- 
tributions to this subject in the open and classified literature 
have been many and varied. However, a gap still existed 
between theoretical studies and practical usage. In an 
effort to bridge this gap, the mathematics of three representa- 
tive charge designs have been reduced to usable form. The 
usefulness of this technique has been demonstrated by appli- 
‘ation to all internal-burning charges used by this Division. 
The simple, straightforward approach has been applied to 
the star, wagon wheel, and modified wagon wheel configu- 
rations. 


The Internal-Burning Star 


MAX W. STONE!’ 


A Practical Mathematical Approach to Grain Design 


Ts. 
Rohm and Haas Co., Redstone Arsenal Research Division, Huntsville, Ala. ie ned 


_—e of its versatility and prominence as a charge con- 
figuration, the internal-burning star was chosen first for anal- 
ysis. The star design shown in Fig. 1 may have desirable 
ballistic properties, but is not practical because thermal 
stresses will cause cracking of the grain at the points. Fig. 2 
shows a star design with small circular ares replacing the sharp 
points; this configuration will be considered here. 

-i 


Fig. 1 Internal-burning star 


The goal was a mathematical analysis that could easily Tig 4 


be used for the selection of grain designs. The techniques 
of geometry, algebra and trigonometry were chosen. 


Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif. June 10-13, 1957. 

1 Mathematician. 

2 Also referred to as H-R design. 


’ Numbers in parentheses indicate References at end of paper. oe Fig. 2 Internal-burning star with stress-relief fillets 
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Fig.3 Sector of internal-burning star-side burns i 


Assuming a constant cross section from one end of the pro- 
pellant grain to the other, and neglecting the effect of erosive 
burning, it is convenient to use a cross section for all calcula- 
tions. In addition, it is seen that there are 2N sectors sym- 
metric with the one indicated by the dotted radii in Fig. 2, 
where N is the number of star points. Therefore, the geo- 
inetric model, shown in Fig. 3, is the sector including one- 
half a star point with its attendant propellant mass. Be- 
cause the author’s experience has been with grains which were 
inhibited on the ends so that burning occurred only on the 
internal surfaces, this is the type burning considered here. 
The quantities to be used are defined in Fig. 3 and the follow- 
ing list. Final surface is understood to be the propellant 
surface existing at web burnout. 


N = number of star points 
D = outside diameter of the grain = inside diameter of the 
motor for case-bonded grain 
w = web = shortest burning distance 
r = radius of circular are at star point 
@ = angle which subtends half the original are of propellant 
at star point 
8 = angle which subtends half the are of propellant at web 
burnout 
B =7/N 
final propellant surface 
PR = progressivity ratio = 
initial propellant surface 
€ = loading fraction = 
total cross-sectional area of propellant 
cross-sectional area of motor 
SF = sliver or tail-off fraction = 


cross-sectional area of unburned propellant at web burnout 


cross-sectional area of motor 


During burning the surface will regress in parallel layers, 
and, depending upon the initial shape, the final surface may 
be in an are subtended by 6 (Fig. 3) or an are subtended by 
a plus a straight line segment (Fig. 4). Calculations con- 
cerning the configuration at burnout (SF and PR) will have 
to be treated separately for the two cases. The propellant 
present at web burnout may remain as sliver or it may burn 
at a reduced pressure to produce tail-off on the pressure-time 
trace. 

It should be noted that loading fraction and sliver both 
have the motor area as their point of reference. If it is de- 
sired to know sliver as a fraction of the propellant, SF may be 
divided by e. 

The existence of thermal stresses in star-perforated grains 
(and others as well) is well known, and has resulted in much 
experimental work to determine their cause and the correc- 
tions needed. Since it has been demonstrated that star 
points having no fillets lead to cracked grains during curing 
and/or temperature cycling, it is now customary to replace 
with ares all sharp indentions in the propellant surface. 
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Fig. 4 Section of internal-burning star-side does not burn out 


Thus, the stars used here have been equipped with circular 
arcs tangent to the sides of the point. All calculations made > 
in this study use 0.03-in. as the are radius for a 1-in.-diam 
grain. For larger grains the radius is linearly scalable. — 
Other methods for relieving the stresses could be employed 
as well, as far as computational technique is concerned. 

The equations relating grain geometry to the ballistic 
parameters of e, SF and PR were written on the basis of geo-— 
metrical considerations evident in Figs. 3 and 4. Thus, a | 
series of equations‘ were obtained which relate the grain de- _ 
sign parameters and the ballistic parameters of e, SF and 
PR. Since the angle a has no ballistic significance, it will 
be preferable to specify values of the other parameters and 
solve for a. Although an analytical solution is not obtain- 
able, the equations readily lend themselves to the iterative 


4 The equations are given in the Appendix. 
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Fig. 5 Design graph for internal-burning star; NV = 6 
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RADIANS 


€ 
6 Design graph for internal-burning star; NV = 


| 


procedures commonly used by large calculating machines. 

The equations were programmed by the Redstone Arsenal 
Computation Laboratory® for the IBM 650 and computed 
for a wide range of the parameters N; w and PR. The re- 
sults were graphed in the form shown in Figs. 5, 6 and 7, 
with a different graph for each value of N from 2 through 
16 (4). All calculations were based on D = 1 in. and r = 
0.03 in.; however, the results may be extended to any size 
motor by scaling on the diameter. All angles and ratios are 
unchanged while linear distances are multiplied by the 
diameter. 

There are two lines on the graphs which must be explained 
further. The theta curve represents the locus of points for 
which @ = a, the point at which the straight line part of the 
star point disappears. The value of 6 is read on the a scale. 
Designs chosen from the concave side of the theta curve will 
be such that the straight side of the star point will disappear 
before web burnout. Conversely, designs from the convex 
side of the curve will have the straight side remaining at web — 
burnout. Because the graph is crowded near the top, the — 
theta curve ends in an arrow, indicating that it does not cross 
the next web line as it continues upward. 

The other unexplained line on the graphs is a horizontal 
straight line. In the case N = 7 (Fig. 6), for example, it 
appears at a = 1.3992 radians. This line has to do with the 
shape of the surface-time trace for a particular design. For 
example, from the design graph, a design can be chosen which 
seems to be ideal, yet the surface-time history may be as in 
Fig. 8. This may happen in the case where the straight 
sides of the star point burn out before the web does, so that 
there are two distinct phases of the surface progression | with 


5 Now part of Army Ballistic Missile Agency. 
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Fig. 7 Design graph for internal-burning star; N = 8 
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Fig. 8 Undesirable surface-time trace 


each having its own separate progressivity ratio. Equations 
have been written to determine the value of @ required to 
give a certain PR, all other parameters being known, but the 
equations are too complex to be readily used because the 
second phase is nonlinear. A compromise solution was ob- 
tained by calculating the value of a required to raise the first 
portion of the surface-time trace to neutrality (PR = 1.0). 
This condition is simply a = tan (a — £), a constant for 
each N. This is the value represented by the horizontal line 
on each design graph. Now, for example, with that par- 
ticular value of a just mentioned and a PR = 1.1, the surface- 
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time trace will be as shown in Fig. 9. Dropping below the 
horizontal line increases the progressivity ratio for the first 
phase, going above the line decreases it. This information 
yields a qualitative picture of the pressure-time trace which 
can be checked by actually plotting the surface-time history. 

Surface as a function of x, the distance burned, can be 
plotted by means of the following equations 


S; = 2NDI x 


2 cos @ 
—(r+ 2) tana 


cos (a — B) 


+2) a + 
[1] 


where / is length of grain. Equation [1] is applicable until 
the second term becomes equal to the last (third) term, indi- 
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Fig. 9 Improved surface-time trace 
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Practical minimum value for angle a 
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cating that the side has disappeared. Then Equation [2] 


is valid. 
w) sing! 


r+2z 


= 2NDI (r + x) | + are sin 


[2] 
If only the shape of the trace is desired, the multiple 2VDI_ 


may be omitted. 
It should be noted that the pressure-time trace will not 


necessarily have the same progressivity ratio as the surface- — 


time trace because of the effect of the pressure exponent 
(6, 7). In addition, erosive burning, if present, will change | 
the burning pattern. 

There is a minimum value for a, shown in Fig. 10. This 
means that the inner projection (the base) of the star point is — 
the same distance from the motor wall as the star point itself. 
Since the result is 2N points with only alternate ones possess- 
ing the stress relief fillet, a more practical minimum is shown | 
in Fig. 11. This is outside the region of normal interest, as — 
seen from Figs. 5, 6 and 7, and is of no real concern. The 
maximum value for @ is indicated in Fig. 12. These limits 
were observed in the construction of the graphs. 


Example 1 
Required: 
PR—neutral, no dip or peak in surface-time trace 
approximately 75 per cent 
SF —low as possible 
-w —thin web, exact value not critical 


Solution: There is only one design for each value of N 
that gives a completely neutral trace; that design is found at 
the intersection of the horizontal line with the curve PR = 
1.0. From Figs. 5, 6 and 7, it is immediately apparent that 
only at N = 7 is the € requirement satisfied. At that point 
SF is 123 per cent and w = 0.144 in. 


Example 2 
Required: 


-PR—nearly neutral, slight progressivity acceptable 
—high, over 85 per cent 

—low as possible 

w —fairly thick web, in the neighborhood of 0.2D 


Solution: From the graphs the following design was 
chosen: N = 6, PR = 1.05, e« = 85.1 per cent, SF = 11.1 
per cent, w = 0.19 in. The surface-time trace is neutral for 
about 80 per cent of the burning distance and then becomes 
slightly progressive. 
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Fig. 15 Design graph for internal-burning wagon wheel; NV = 5 


The Internal-Burning Wagon Wheel Design 


Another basic configuration that has proved useful in sev-— 
eral slightly differing forms is the wagon wheel design, shown — 
in Fig. 13. The thickness of the spear of propellant is twice 
the web, that is, twice the shortest distance between inner | 
and outer surfaces of the grain. The two corners of the 
spokes (the port) are circular arcs tangent to the sides of 
the spears and to the short are between them, which is con- — 
centric with the grain perimeter. 

This configuration received the same treatment as the star. — 

Many of the same definitions and designations apply; the PR iool- 2.07% S05 | 
new ones may be seen in Fig. 14. It would be possible to 
calculate the length of line L if z were specified, but it seems 
more reasonable to include ZL in the equation for progressivity — 
ratio, specify a desired PR and solve for L. This was the — 
procedure followed in the calculations that have been per-— 
formed. 

The equations for the wagon wheel were programmed for 
the Redstone IBM 650 and solved for a wide range of values - 
of N,wand PR. As before, D and r were constant at 1 and 
0.03 in., respectively. 

The equation for the line L indicates that there is no wagon — 
wheel design giving a neutral surface-time trace (PR = 1.0) 
because the quantity (PR — 1) is used asa divisor. For all — 


practical purposes a neutral design can be obtained because — fi 

values of PR very close to 1.0, e.g., 1.01, were used success- 

The computer results were plotted in the form shown in” LN) 

Figs. 15 and 16 for values of N from 2 through 16 (5). Use _— 16 Design graph for internal-burning wagon wheel 
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of the graphs is the same as with the star, with the exception 
that L must be found from the second graph or hand com- 
puted from the equation after the design is chosen. The 
wagon wheel is a thin-web design, with its maximum being 
about 0.2D, while the star can go up to about 0.33D. In both 
designs sliver fraction increases with web; the star sliver frac- 
tion decreases with increasing N while the wagon wheel 
sliver fraction increases. Thus, of the two, only the wagon 
wheel provides small values of SF for small N. 


~ 
Fig. 17 Modified wagon wheel design 
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Fig. 18 Three-step surface-time trace 


Fig. 19 Sector of modified wagon wheel design 
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As one would expect, there is a maximum value of web that 
can be put intoa motor. This value can be found from 


+) sina =r 


Example 3 


Required: 
PR—nearly neutral, slight progressivity acceptable 
—-80 per cent 
SF —low as possible 

w —thin web 


Solution: Although only the graph for N = 5 is reproduced 
here, all traces for N = 6 and above will be regressive; for 
values of N less than 5, traces become increasingly progressive. 
Therefore, the only chance for a high loading fraction in a 
nearly neutral design is at N = 5. Thus from Fig. 15 is 
chosen the design NV = 5, € = 80.2 per cent, SF = 9.8 per 
cent, w = 0.112 in., PR = 1.04, and from Fig. 16 the line Z is 
found to be 0.15 in. 


Modified Wagon Wheel 


This useful design (Fig. 17) is a wagon wheel design with | 
alternately long and short spears of propellant. In this anal- 
ysis N will refer to the number of long spears (equal to the 
number of short spears) but not the total number of spears. 

A major innovation in the analysis of this design is the treat- 
ment of the three burning distances as not necessarily being 
equal. The three burning distances are determined by the 
thickness of the long and of the short spears and by the short- 
est distance from inner to outer propellant surfaces. This_ 
gives the possibility of a one-level (all burning distances _ 
equal), two-level (two burning distances equal) or three- | 
level (all burning distances unequal) thrust operation (Fig. | 
18). (The wagon wheel could be treated in this way also.) © 
The relative heights of the levels in the second and third cases _ 
‘an be adjusted somewhat by the choice of the long spear 
with its larger mass of propellant for the first step, or the — 
short spear with its smaller mass, or by choosing which of the © 
webs are to be equal. 

Some new lines and angles are added in this analysis; they — 
are indicated in Fig. 19. Because of the web treatment, the _ 
equations arrived at are more numerous and appear more _ 
formidable than previous ones. A compromise to practicality 
was made in defining the configuration so that the distance z _ 
(Fig. 19) applies between two spears as it does between the © 
long spear and the center lines. Thus, the spear with K as its — 
side is always the short one. 

In the same manner as before, equations were written to 
relate the grain geometry to e, SF and PR. These equations 
were then programmed for solution on the IBM 704 at the 
Redstone Arsenal Computation Laboratory, and solutions — 
obtained for a wide range of values of w;, we, w3, N and PR. 
This work has just been completed and the results will be 
graphed in a form similar to the other design graphs. 
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Equations for Grain Designs 


Internal-Burning Star 


The equations are derived directly from the geometry of 


: 
9 
9 : 
| 
1 78 
NS 
is SR 
S 
| / fk 
=~ 4 z - 4 


the design (Figs. 3 and 4) in terms of the definitions given in Actually, the final surface is S; + S,, but up until the in- 
the text of the article. For the case where the straight sides stant of burnout the spear of propellant remains. A separate 
of the star point disappear before web burnout, the following equation for line L cannot be written without also specifying 
equations are true z; it is more desirable to specify PR and solve for L in Equa- 
tion [10]. 
fa S; + S,+ — m — (S, + Se + m) PR 
2 cos 
+ r? (a — tana) : 
= cos (a — ) Expressions for the are lengths are easily obtained and wi! 
_[A-1] be omitted here. The area ratios are found by 
5 4w? 4 2 
tan B _ tan (8B — a) 
=D? .. [A-12 
N N 
4 
dD? (8 a) A(r 4+ w)? 6 tw) 
anc ‘ => = A-15 
A motor wD? 
N 
Modified Wagon Wheel Design 
> (r + w)@ ‘ 
PR = D 3 ... [A-2] Again the equations are found in a straightforward fashion 
- <oe + —.) sin B The angles appearing in Fig. 17 can be calculated by 
cos @ 
cos (a — B) aresin [A-14 
(? -r- v) sin B 7? Y = are sin [A-15} 
SF + w) 0+ (r+ w) v) sin 3] 6 = 90° +B —a@.............. [A-16] 
us 2 
[A-4] = 90° + ¥ [A-17] 
In the case where the side does not disappear during burn- — are cos [A-18] 
ing (Fig. 4), Equation [A-1], along with the following equa- r+ wy 
tions for progressivity ratio and sliver fraction apply 
D 
( sing = 
2 cos @ 
r(a — tana) + wla — tan (a — B)| + — ( B 
PR = ... [A-5] 
r 
r (a — tana) + 
cos (a — B) 
sp | tana) X [a tan (a 8) + 
xD? cos (a — B) st 


Internal-Burning Wagon Wheel Design 


From the geometry of the grain as seen in Fig. 14, the fol- 


lowing equations can be written t ‘ 
= — are cos [A-19] 
r+w r+ ws 


= — arcsin [A-7] 
= =P my The are lengths S; through Sg are calculated simply as the 
par e product of ungle and radius so those equations will not be 
ar The length of the line L will be determined from the PR 
oe oe [A-9] equations which follow, similar to the wagon wheel design. 
Since K is a function of L, it can be found by Equation [A-20]. 
w r+ wr Ws We r+ ws 
tan B tan Y tan (8 — a) 


. [A-20] 
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> in- In order to distinguish practical values of L, the following surface B = final surface = ~— ie : 
rate equations must be solved. The maximum value occurs when ws 
= L—m + — + (rt w3)0 + Sy + Ss... [A-32] 
ying tan B 
tan(8 —a) tanB 
and the minimum at k = 0 S; [ tan B 
—. tan8 tan (8 — a) tan Y tan B surface AB = surface A = surface B = 
will .. « (ee) bracketed term of Equation [A-33]. . [A-34] 
In a multistep thrust or pressure trace, it is obvious that surface BF = (r + w:) (0+ o) + Ss...... [A-35] 
FR as previously defined is inadequate. To be meaningful a . a 
a Separate PR must be calculated for each step. Thus, from 
lig. 16 the following definitions present themselves > > 
1 We — W3 W3 
A 4 = 
initial surface tan( 90° — 
surface at B 
‘13 surface at AB (r + ws) (0+ 1... {A-371 
tan B 
] surf: 
[A-25] 
surface at BF surface AB = K + —- 
sin 
The equation used to calculate PR, for any particular case sie : 
on vill depend upon the relationship of the webs to each other. Tay + (r + w;) (6+ o) + Sz... [A-38] 
The following equations give PR, for each case and also the tan (0° ’) 
intermediate surfaces from which PR, and PR, can be cal- 2 
14 culated by Equations [A-24 and A-25]. The denominator 
of the PR, equation is the initial surface (S,) and is the same surface B = K —n+ + (r+ (@+0-W4+S8s 
in every case. [A-39] 
PR, = PR, = PR, = _ Vir + (r + ws)? 
6} 
[A-26 
17] surface BF = (r + — + [A-41] 
wy = > ws final surface = Sy + S, + [A-42 
18] PR, K + + Ws Ws + > W3 > We 
Si a sin 8 B 
90° — l Ws 
—i an 
(r + ws) (0+ ¢) ...[A-27] 2 
nB (6+ ¢) +— 
sin B .# 
— wr tan{ 90° — 
surface AB = surface BF = K + ——— - 
sin B ws — ws 
7 surface AB = L+— sin B _ 
W3 
+ + we) O+ + Sr... [A-28] 
B 
surface B = final surface = 
+ ¢ ..[A-29] surface B = L — m + (r + ws) (0 o- 2) +8; 
= W3 > [A45] 
9) 1 where 
PRa= K+L—nt+ 6+ 9) + 
V(r + wa)? (r + ws)? 
be + Ss |... [A-30] : 
sin B (s0° surface BF = + ws) (0+ — Q) + [A-47] 
AB =s ino BF = L + 
)]. surface =surla sin B Although not necessary in the calculations, the value of 2 — 
. will aid in drawing the design. 
2) (0 Ss... [A-31 
+ (r + we) ( + $) + 8 [ ] (sin B) W3 
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The area ratios are 


= nB + + 2L w; — r? (0@+¢) + 


tan tan (8 — a) (; 
2 
SF = =. [> — (r+ w,)?A — (r + we)? tan (6 — A) + 
4 
D?(8@—a) (r+ wu)? (r + ws)? 
4 tan Y tan a) 
(r + ws)? tan (8 — » |. .[A-51] 
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Flight and orbital environmental conditions exert a 


Flight Environmental Conditions 


great influence on satellite design. As information on : hae? i 
these conditions changed and became firm the satellite The flight events that precede the ejection of the satellite th 
design went through numerous stages. Material selected from the third stage may be described as follows (Fig. 2). in 
for the satellite shell and structure is magnesium alloy > la 
FS-1. Design and operation of the structure and 7 of 
shell, antennas and separation mechanism are discussed. ; rape SATELLITE al 
Methods for control of satellite temperatures are thermal i as . 
isolation of the inner package, shell coatings and thermal oe) di 
switches. From a circular blank disk, a near-perfect ta! VANGUARD ” st 
hemisphere, only 0.026 in. thick, is formed to make the ti 
LAUNCHING VEHICLE SECOND 
satellite shell. Tubes are formed and welded to make up STAGE 15 
the “‘bird cage”’ internal structure. Fitting and assembly 
of the satellite involve jewelers’ precision on the part of +s = 
several groups. Satellite structural components are de- ; ad ll 
sign-tested individually at first. The assembly is then ae 7 en 
tested to the expected environmental conditions of ace th 
celeration, vibration, vacuum and temperature. ge 
STAGE 
N JULY 1955, President Eisenhower formally announced ba 
that the United States would undertake the launching of e : pa 
artificial satellites as a part of the research program for the th 
International Geophysical Year. One satellite has already th 
been launc hed by the U. S. Army. To place the Navy’s r Fig. 1. Vanguard launching vehicle wi 
small 20-in.-diam 211/.-lb satellite in an orbit around the »* be 
earth takes a three-stage rocket vehicle approximately 72 ft an 
long, weighing over 10 tons. Such a vehicle (Fig. 1) may INITIATE. 3RO STAGE SPIN THIRD STAGE ae 
be described as follows: The first stage is a liquid-fueled SEPARATE 2ND_ STAGE IGNITION SEPARATION 
rocket with a 27,000-lb-thrust rocket motor. The second gp OPTIMUM TRAJECTORY $ 
stage, also a liquid propellant rocket, attaches to the forward oy 
end of the first stage and houses the third stage and satellite . 
during its flight. The second stage contains the “brains” ; rs lin 
for the entire launching vehicle. The third stage is a solid Py hs > VELOCITY 25,000 FT/SEC fol 
propellant rocket. The satellite, mounted on the separation ‘sie ae - RANGE 1500 MILES Sif 
mechanism, is attached to the third stage rocket. = 
28 
we ta 
_ | Revised edition of paper presented at the ARS Spring Meet- rs 
ing, April 4—6, 1957, Washington, D. C. 1 
? Head, Satellite Structures Group, Project Vanguard. Mem. ‘ ; sti 
ARS. Fig. 2. Vehicle trajectory ny 
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(a) The: vehicle goes through a short and then 
starts a zero-lift trajectory. (b) Approximately midway 
through this zero-lift phase the first stage burns out, separa- 
tion of the second stage from the first stage occurs, and the 
second stage is ignited. (c) The second stage is controlled for 
optimum trajectory through part of its powered and all of its 
coasting flight. (d) The third stage is spun and the second 
stage is separated. (e) When the second stage is clear the 
third stage ignites and is spin-stabilized throughout its burn- 
ing period. (f) At burnout, approximately 10 min after 
launching, the satellite and the third stage will have a velocity 
of approximately 25,000 fps, an altitude of 200 to 400 miles, 
and a range of approximately 1500 miles. 

The satellite must survive vibrations imposed by three in- 
dividual rocket motors. It must survive the shocks of rocket 
starts. In addition, it must withstand steady state accelera- 
tions up to 35 g, aerodynamic heating to 300 F and rotation to 
150 rpm. 


Orbital Environmental Conditions 


The difficulty encountered in trying to design for orbital 
environmental conditions lies mostly in the uncertainty of 
these conditions. One of the major roles of the satellite is to 
get measurements of these conditions. Not much is known 
about ambient temperatures, meteors, space dust, cosmic 
rays, x-rays, ete., at altitudes from 200 to 1400 miles, the 
space in which the satellite must travel. Of the unknowns, 
temperature seems to stand out as the major problem. The 
batteries and transistors to be used can only operate in a com- 
paratively narrow temperature range (0 to 60 C). The orbit, 
the time of year of firing and the hour of firing contribute to 
the variations in temperatures expected. The satellite skin 
will fluctuate through a wide range of temperatures. It will 
become extremely cold when the earth is between the satellite 
and the sun. It will become hot when it is on the sunny side. 


Design 


me 


Design criteria had to be determined on the basis of pre- 
liminary estimates of the magnitude of the powered flight 
forces. Design factors of safety to be used in the de- 
sign calculations were then established. Before actual design 
could begin, the basic shape of the satellite had to be es- 
tablished. Some of the contributing factors that determined 
the size and shape were the maximum size that would fit into 
the vehicle with a minimum weight penalty for the second 
stage, minimum size needed to allow effective optical track- 
ng in the orbit, and minimum weight of shell and structure. 
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General 


Fig. 3 (a) Satellite with antennas in the folded position. 


}-in. OD by 0.024-in. wall aluminum tubing 6061ST6. They 


(b) Satellite with antennas in the 
radial position 


Material Selection 

The most important factor, other than reliability, that con- 
trols design is weight. For every pound the satellite weighs it 
takes approximately 1000 Ib of rocket vehicle to give the | 
necessary velocity to orbit. Every gram of satellite structural 
weight decreases the scientific payload the same amount. It 
was necessary, therefore, to make a very thorough study of 
materials for use in the satellite shell and structure. FS1 
magnesium alloy was chosen as a shell and basic structure 
material because of its light weight and also because it is one 
of the most commonly used and easily worked magnesium 
alloys. Calculations indicate that adequate strength could 


be obtained from FS1 in minimum thicknesses to be used for 


brics 
fabrication. 


The various loads incurred by the structure during the 
ascending trajectory may be analyzed as follows (Fig. 9): 
The major concentrated load which is the internal package is 
carried directly through the Kel-F support to the main hub at 
the base of the structure. This hub rests on the separation 
mechanism and the load is transmitted directly into the top 
of the third stage. The side loads incurred by vibration are 
transmitted from the internal package through the four Kel-F 
side rods. These rods transmit the load to the circumferential 
rim on the major axis of the satellite. Part of the side loads 
are taken up by the vertical bow-shaped tube, and part is 
transmitted around the ring into the shell of the satellite 
itself at the antenna supports. The major antenna load is 
carried by the structure but a portion of the load is trans- 
mitted through the shell. The load of the shell is picked up at 
the top ring of the structure, at the major or equatorial axis 
and also at the base. The highest stresses in the shell are 
located near the base of the satellite; hence the skin is thicker 
there. This section covers about 103 in. in diam at the base 
of the sphere and is 0.050 in. thick. Thinner skins (i.e., ap- 
proximately 0.025) show an oil-canning effect when being 
vibrated in the horizontal direction and fatigue quite rapidly. 
The frame is designed to take up to 63-in. maximum diam by 
134-in.-long internal package. This affords a considerable 
amount of flexibility for the internal package design. The 
method of base support and side support is such that it is 
relatively easy to incorporate different shapes and sizes of in- 
ternal packages. 


Structure and Shell 


Antennas 


The antennas are 29 in. long and constructed primarily of 


245 


a i 
1 
(a) 

= 

i 
IRE 


Fig. 4 Satellite in Coating Facility at ERDL, Fort Belvoir, Va. 


are operated by a spring which is a tension. The antennas 
are folded forward while the nose cone is on the second stage 
(Fig. 3 (a)). When the nose cone is jettisoned the antennas 
will rotate more than 90 deg to the radial position and lock 
on a tapered surface (Fig. 3 (b)). a ee 


Temperature Control 


The shell temperature can be controlled to a degree by its 
surface coating. The satellite coating must satisfy two funda- 
mental requirements. It must be highly reflective, in the 
visible, and it must emit absorbed solar energy to stabilize its 
temperature. Such a coating has been developed.* An evap- 
orative coating of aluminum provides high reflectance while 
an outer coating of silicon monoxide, transparent in the 
visible, provides the required emissivity in the infrared. Fig. 
4 shows the coating facility of ERDL, Fort Belvoir, Va. 

The internal package containing the electronics and bat- 
teries must be thermally isolated from the shell. Only two 
methods of heat transfer are involved—radiation and conduc- 
tion. Radiation transfer is minimized by using gold-plated 
surfaces (low thermal emissivity) and conduction is minimized 
by using low thermal conductance structural supports to con- 
nect the internal package to the shell. The isolation achieved 
in this manner will keep the internal package temperature 


3 Bradford, A. P., Cox, J. T., and Hass, G., “Coatings for In- 
frared Optics of the Satellite,’ U. S. Army Engineer Research 
and Development Laboratories, Fort Belvoir. Va. 
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Shells 


within 1 C of the mean orbital temperature. The shell tetn- 
perature will fluctuate +15 C about its mean orbital tempera- 
ture after equilibrium has been reached. 

A very thorough analytical and experimental program is 
being conducted on the thermal phase in order to get the best 
heat balance possible. If by chance estimates of orbital tem- 
perature conditions are extremely far off the result would |e 
that instead of the two-week designed life for the electronics 
we might get data for only a few days. These data woul, 
however, tell us where the estimates were in error and give 
the next unit a far greater chance for a full life. 


The Satellite Separation Mechanism 


Fig. 5 shows the separation mechanism which is a con - 
pletely self-contained unit that fits in a socket at the base 
of the satellite. It uses the characteristics of the third stage 
as a method of arming and operation. The mechanis: 
operates as follows: If approximately 12 g or more are applic 
for over 2 sec, the clock will be allowed to run for 10 see mor» 
to a stop on the g-weight arm; if 12 g or more is applied fi 
less than 2 sec, the unit will reset itself. Then after the g lev: 
is reduced below the 12-g level at third stage burnout th: 
timer will start to run and run approximately 26 sec. At the 
end of this period the timer arm closes circuits to the cater- 
pillar motors, which expand, pulling the hold-down pins i 
and rotating the spring release. The spring exerts a force 
which imparts a differential velocity of approximately 3 fp 
between the third stage and the satellite. The unit is d- 
signed with two completely independent circuits. Two RM-| 
batteries are used. One caterpillar motor is fired by each 
battery independently. Upon closing of the circuit by the 
timer arm, the batteries furnish enough voltage to the wire 
wound squibs in the caterpillar motors to cause them to ignite 
the charge and the expanding gases do the work of pulling the 
pins and rotating the spring release. Either of these circuits 
is capable of firing and releasing the satellite. The purpose of 
the duplication is to insure a higher degree of reliability. 


> 


4 


Fabrication‘ 


A circular blank disk 0.091 in. thick is cut out of sheet stock. 
This disk is then drawn with one die in one stroke of the press 
with the material heated to 650 F. After the flash is trimmed 
off, the drawn hemisphere is placed on a spinning lathe. The 
metal hemisphere is spun hot (600 F approximately) to give 
exact internal dimensions required. The inner segment of the 
pressure zone is spun from a flat circular blank of magnesium. 


4 “Manufacturing the Earth Satellite,” The Magazine of Mag- 
nesium, May 1957 edition. 
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one The desired section is trimmed out of the spinning. The two 


era- pressure connection fittings are welded in position and then 

the inner segment is welded onto its proper location in the 
n is interior of the hemisphere. Fig. 6 shows the resulting pressure 
best cavity. The unit is then stress-relieved and hydrostatically 
eli- pressure-tested for leaks. When the unit is found to be leak- 
1 be free it is placed on a mandrel and the outer surface is machined 
nics to give a skin thickness of 0.026 in. except at base section 
uli, where it is 0.050 in. thick. The hemispheres are then trimmed 
rive and all drilling, tapping and other machine work is done. 


The hemispheres are then ready for finishing and plating. 


The Structure 


yD - The tubular sections that make up the structure are first 
ase bent into the shape desired. These pieces are then machined 
ape as necessary in order to obtain the proper fits. The entire 
is}. assembly is fitted to a special jig and positioned so that it 
ie] will meet specifications and tolerances as to the concentricity, 
ore dimensions, etc. The tubular structure is then welded to- 
for gether while clamped in place on the jig (Fig. 7). With the 
vel siructure still on the jig the spherical surfaces at the top, 
th» bottom and where the antennas are located are then ma- 
the chined to fit the inside dimensions of the shell halves (Fig. 8). 

1 Assembly 
* After completion of the structure the shell halves are then 
... mounted onto the structure and secured in place by rivets and 
[-| screws. After the shells are assembled to the structure and 
ch inspected the unit is then ready for shipment to NRL where 


he further assembly of the scientific equipment and the coating 
of the shell will take place. 

The “Kel-F” support pieces are made at NRL and are on 
he hand when the unit arrives. The shells are taken off the 
structure and sent to Fort Belvoir for the silicon monoxide 


ts 
oF coating process. The internal structure is sent to the NRL 
shops for fitting of the internal package with the Kel-F sup- 
ports and the assembly of the antennas. Upon return of the 
shells from coating, the lower half of the shell is assembled to 
the structure. Gages, such as erosion gages, temperature 
gages, pressure gages, and so forth, are mounted on the lower 
hemisphere. The internal package is assembled; the an- 
kK. tenna connections are made; the wiring is done on the lower 
SS half of the sphere and the upper half of the sphere terminating 
ec| in the plugs previously mentioned. There is an extensive in- 
he spection and test procedure performed on each and every unit 
ve prior to, and after, the assembly procedure. 
he 
Fig. 7 Welding satellite structure 
g- Shells with Pressure Zones 


1 Hydrostatic. Each hemisphere is tested hydrostatically 
as it is received from the contractor. The zone is filled with 
distilled water and the pressure is increased gradually until a 
pressure of 25 psig is reached. This test is for the purpose of 
determining the structural adequacy of the zone and also to 
locate any leaks. 

2 Dry Gas. After completion of the hydrostatic tests the 
hemisphere is dried and filled with freon gas to approximately 
20 psig. A leak detector is employed to locate minute leaks. 
Leaks, if any, are marked and the zone is then flushed with 
nitrogen gas. 

3 Vacuum. The zone is connected to a vacuum system 
and the pressure is reduced to 0.6mm Hg vacuum. A sealer is 
| introduced into any small holes for sealing the leaks. The 

hemisphere is then oven dried again. Steps 2 and 3 are re- 
peated until no leaks can be detected. hes \.ehe 


Satellite Separation Mechanism 


In general the purpose of the tests on the satellite separation 
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mechanisms were to determine their suitability as separation 
units for separating the satellite after completion of third 
stage burning. The units must withstand the environmental 
conditions of the launching vehicle, provide sufficient delay 
and sufficient force to separate the satellite from the third 
stage after third stage burnout. 


Antennas 


The tests on the antennas were for two purposes. First, to 
check out the antenna design to insure that the antenna will 
survive the environment which exists during flight, such as 
aerodynamic heating, vibration, and acceleration. The second 
set of tests are checks to determine whether or not the anten- 
nas have been fabricated properly and have the strength as 
required, i.e., checking the heat-treatment of the antennas. 


; Assembly Design and Final Acceptance Tests 


Design acceptance tests and final acceptance tests are made 
independently under test conditions outlined by the Environ- 
mental Subcommittee of the Science Program Committee at 
NRL. The design acceptance tests consist of: (a) Random 
vibration up to 25 g parallel to thrust from 10 to 2000 cps and 
up to 20 g from 10 to 2000 cps traverse to thrust; (b) steady 
state acceleration to 75 g in the thrust direction; (c) extensive 
thermal-vacuum life tests between 0 and 60 C (internal 
package temperature) on the completely instrumental satel- 
lite. The final acceptance tests of flight units consist of the 
following tests: (a) Random vibration up to 10 g from 10 to 
2000 cps; (b) steady stage acceleration to 35 g in the thrust 
direction; and (c) thermal and vacuum tests as follows: simu- 
lation of initial heating of shell to 300 F, temperature cycling 
similar to design acceptance tests, and life test of one week 
under vacuum and thermal conditions. 

Two complete satellites, identical to flight units, are used for 
thermal testing. These units are instrumented with all 
gages, electronics and batteries. The assembly is placed in a 
vacuum chamber which is housed in a high altitude chamber. 
The environment in this chamber can be controlled. The 
temperature is cycled in a manner which simulates the ex- 
pected orbital temperature fluctuations. The pressure is re- 
duced to 0.1 mm of mercury for these tests. The temperature 


of the shell, the inner package and electronics are continually 
monitored. The equipment is operated as it will be in orbit 
and monitored constantly. Such a test as this takes about 
three weeks if no failures of equipment occur. = oe 
Summary 


After the results of the many and varied tests have been 
studied and analyzed and the design is considered satisfactory 
to meet all requirements, the resultant satellite (Fig. 9) may 
be described as follows: The highly polished silicon monoxide 
coated 20-in.-diam magnesium sphere has four antennis 
mounted 90 deg apart at the equator. These antennas are 
fastened by tubular rods to a tubular ring which is concentri:: 
with the cylindrical inner package. This ring is also supporte:|! 
by four bow-shaped vertical tubular members spaced 90 de. 
apart at angles of 45 deg to the antenna supports. Thes 
vertical members are fastened to the support ring of the ac- 
cess port at the top and to the main support column whic/ 
houses the separation mechanism at the bottom. Th 
cylindrical internal package is mounted on the main suppor 
column and is also secured to the concentric tubular ring b) 
four low thermal conductivity radial supports spaced 90 deg 
apart and angularly positioned between the antenna sup 
ports and the vertical members. The sphere is girdled b: 
two pressure zones or bands, one below the equator and on: 
above. On the shell are three microphones, a Lyman Alph: 
solar cell and two Lyman Alpha ion chambers. In addition to 
the microphones and the Lyman Alpha experiment, there arc 
experiments which consist of erosion gages and temperaturc: 
gages, which are located both at the north pole area anc 
around the equator. There is a pressure gage located on th: 
concentric tubular ring. Also included in the internal pack 
age are five 3-in.-thick modules; of these one module is 
for the Minitrack transmitter and associated electronics; onc 
module is for the Lyman Alpha electronics and batteries; one 
for the coding 48 channel telemetering system; one for the peak 
memory orbital switch unit; and one for the meteor counter. 
Below these five 2-in.-thick modules are the battery packages. 
One battery package is 14°; in. thick; the other battery pack- 
age is 2} in. thick. These seven modules are held to the top 
cover of the internal package through two }-in.-diam rods. On 
the top cover of the internal package are several connectors. 
One 36-pin connector; one 14-pin connector; one 50-ohm 
microdot connector and two standoffs for the Lyman alpha 
experiment. The connectors serve two purposes—not only 
do they connect the gages on the skin to the electronics, but 
they also connect the batteries with the electronics, thus act- 
ing as a turn-on switch. 


Conclusion 


When the Vanguard launching vehicle places the small 
IGY earth satellite sphere into orbit, it will have placed a re- 
liable laboratory into space. This laboratory together with 
the Sputniks and the Explorer will be instrumental in unlock- 
ing some of the hidden secrets of our planet and the space 
which surrounds us. Perhaps it will be an important instru- 
ment by which man can insure his survival when he first at- 
tempts to venture into the darkness beyond. 
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W. M. BROBECK,! R. E. CLEMENSEN? and W. 


An instrument for the measurement of flame tempera- 
tures from 1500 to 2700 C by the sodium-line reversal 
method is described in which the eye and hand of a trained 
operator are replaced by a closed-loop electromechanical 
system. The advantages are in continuous recording of 
varying flame temperatures, increased accuracy, and 
measurement on flames that would be dangerous to an 
observer. 


Introduction 


HE sodium-line reversal method of measuring high gas 

temperatures (1)* has been widely used and accepted. 
In this method a sodium salt is introduced into the hot gas, 
and the brightness of the characteristic D lines of the sodium 
spectrum is compared with the background formed by the 
continuous spectrum of a lamp filament observed through the 
same gas. With a fixed gas temperature, as the filament tem- 
perature is raised the sodium lines are first brighter and then 
darker than the background. At the “reversal point,” at 
which the lines and background have the same brightness, the 
temperature of a black body having the same brightness as 
the filament is equal to the gas temperature. By observing the 
filament with an optical or thermopile pyrometer, the equiva- 
lent black body temperature that equals the sodium-line re- 
versal (SLR) temperature of the gas can be obtained. 

The temperature reading is independent of the path length 
or concentration of sodium in the gas over a wide range. The 
effect of insufficient light is to make the location of the re- 
versal point less definite but not to move the reversal tem- 
perature, since the reversal point depends on the ratio of ex- 
cited to unexcited sodium atoms. — discussion of the 
method and its accuracy are given in (2, 3, and 5). 

In the usual application of this method to steady tem- 
peratures, the reversal point is observed by an operator look- 
ing through a spectroscope and adjusting the lamp filament 
manually. Methods have been devised for displaying the re- 
versal temperature on film (3) and of following the tempera- 
ture in the combustion chamber of an engine (4). 


2. 
Design 

In the application to be described it was desired to present 
the gas temperature on a conventional self-balancing potenti- 
ometer having a full scale travel in the order of 1 sec. 

The arrangement used is shown diagrammatically in Fig. 1 
Light from the filament B passes through condensing lens C 
and through the gas stream E. From the gas stream the light 


passes through the collimating lens G, slit H, prism I, and 


objective lens J of the spectroscope to the vibrating mirror K. 
From K the light, now consisting of the spectrum in the re- 
gion of the sodium D lines, passes through the exit slit M into 
the photomultiplier tube N. A block diagram of the electrical 
circuit is shown in Fig. 2. The output of the phototube is 
commutated by the chopper and fed through the amplifier 
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to the motor which ‘i of 
lamp filament to match the brightness of the sodium D lines 


from the flame. 
Operation 


The spectrum and electrical signals are shown in Fig. 3 
The two sodium lines are not resolved by the photocell. When 
the lamp is cooler than the gas the spectrum appears as in 
(a). As the mirror sweeps the line back and forth across the 
exit slit the phototube output appears as in (b). When 
chopped, the signal (c) is produced, having the necessary 60- 
cycle component to drive the two-phase motor in the increase 
direction. At (d), (e), and (f) are shown the corresponding 
signals for the case in which the lamp filament is too bright 
and the sodium lines appear dark against the brighter back- 
ground. 

The 60-cycle component in the bright light shown in (e) and 
(f) of Fig. 3 is due to nonlinearity in photocell response and 
difference in energy content of light of different colors. It is 
caused by scanning the spectrum in the area of the sodium 
lines as shown in (g). The 60-cycle component can be re- 
duced by reducing the mirror amplitude, and by use of a 60- 
cycle filter on the photocell output. The chopper phase is ad- 
justed to make any residual 60-cycle oe ‘ancel out as 
shown in (f). The signals in (b), (c), (e), (f) are obtained 
with no filter on the photocell output. W ith the filter in- 
stalled, input to the chopper is a smooth sine wave which de- 
creases to zero at balance, and then increases 180 deg out of 
phase when the balance point is passed. If the vibrating 
mirror were driven at 30 cycles instead of 60, the chopper 
would evidently not be required. However, this arrangement 
has not been tried. 

The temperature is observed by the radiation pyrometer 
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(A, Fig. 1) focused on the lamp filament and driving the re- _ 


cording potentiometer. It had been planned to record the lamp 
filament voltage but changes observed in the characteristics of 
the filament made the use of the radiation pyrometer prefer- 
able. The recorder is located several hundred feet from the 
optical system. Controls, not shown, are provided for adjust- 
ing the amplitude and phase of the mirror vibration. The 
effect of these adjustments is observed by a cathode-ray 
oscilloscope displaying the amplifier input. Ps 


LOW _LAMP_ INTENSITY 


BRIGHT LINE SPECTRUM 


PHOTOCELL 
(b) _Sl OUTPUT 
(c) | fl CHOPPER 
OUTPUT 
HIGH LAMP INTENSITY 


LINE SPECTRUM 


(e) PHOTOCELL 
OUTPUT 


CHOPPER 
(f) OUTPUT 


DARK SODIUM ABSORPTION LINES 


OUTPUT 
(g) SIGNAL 


INTENSITY 


[—] MIRROR SCANNING RANGE 


LIGHT WAVE LENGTH 


“Fig. 3. Optical and electrical signals 


190 
it = 
c 
| 
| | 140 
jf 
| 
1500 1680 1800 1880 1950 2010 2070 2115 2155 2190 


Fig. 4 Hysteresis of temperature measurement, propane-air 


The walls of the test section, in which the pressure is a few 
hundred psi, are protected by a graphite tube. The windows 
D and F (Fig. 1) can become obscured by carbon or metal 
vapor released in the test section. Deposits on the window F 
are important only as they affect the sharpness of the balance 
point, but changes in the transmission of D will affect the re- 
corded temperature. For this reason, the radiation pyrometer 
has been calibrated by observing the filament with an optical 
pyrometer through the glass of window D before and after 
several runs. In the present application, the change during a 
5-min run has not been appreciable. The correction for ab- 
sorption in the glass between the lamp and the flame is approxi- 
mately related by k = AT/T?. For a typical lens and view- 
port, k equals 4 X 10~* for the lens C, and 9 X 10~* for the 
viewpoint D. Using these correction factors, the pyrometer 
can be calibrated by an optical pyrometer looking directly at 
the lamp. 

Sodium vapor is introduced into the gas at a point of low 
velocity by a sodium “bubbler,” in which coiled nichrome wire 
filament is used to heat sodium chloride to its boiling point in 
aceramic cup. Sodium vapors from the cup are swept into the 
flame to provide the concentration of sodium necessary in the 
hot gas. 

The usable range of the instrument is limited to from 1500 to 
2700 C. The lower limit is reached when the signal strength is 
so low that the photocell cannot pick it up. The upper limit 
is determined by the maximum usable filament temperature. 
Although tungsten does not melt until it reaches 3370 C, ex- 
cessive vaporization causes obstruction of the light and rapid 
failure of the filament above 2700 C. 


Results 


Peformance of the instrument has been quite satisfactory. 
The balance point of the instrument agrees with that obtained 
by eye. Calibrations of repeated balance points agree to 
within the accuracy of the optical pyrometer, which is +5 C. 
The hysteresis as determined by the return to balance after a 
disturbance is +5 C, as shown in Figs.4and5. The observed 
flame temperatures in Figs. 4 and 5 agree quite well with tem- 
peratures reported in the literature of 1737 C for ethane in air 
and 2067 C for acetylene in air (5). Fluctuations in the base 
lines are due to variations in flame temperature. The deflec- 
tions up and down scale were deliberately introduced to check 
As can be 
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Fig. 5 Hysteresis of temperature measurement, acetylene-air 
flame 


JET PROPULSION 


- 
4: 
fa) 
~ > 
~ 
100 | i | | | 
| 
> 
: 


60~ ] 

MIRROR 
=~» 


(a) MIRROR POWER SUPPLY - 


fee «lew 


POWER AMPLIFIER 


VOLTAGE AMP 
r 


Spi 


Ly 


(b) DETECTOR CIRCUIT 
| 


Fig.6 Electronic circuits 


2 
seen, the speed of the motor drive limits use of this instrument balanced 120-cycle component. Phase of the mirror is adjusted 
to flames whose rate of temperature variations does not exceed to produce maximum sensitivity to 120 cycles and minimum 


50 C per sec. response to the 60-cycle component. 


Light leaving the exit slit enters the photomultiplier of the 
Most of the 60-cycle 
component, produced by the nonlinearity of the spectral sensi- 


Calibration detector circuit shown in Fig. 6(b). 


The recorder is calibrated by observing the lamp filament tivity of the photomultiplier and any high frequency inter- 
from the same position as the recording pyrometer with an ference, is shorted to ground in the filter, which is tuned to 
optical pyrometer. Several temperatures are checked, and pass only 120 cycles. 


the millivolt output of the recorder is correlated against the 


optical pyrometer readings corrected for the effect of lens C, 

viewport D, and the difference in tungsten emissivity be- P ee 

tween 6650 A, where the pyrometer is calibrated, and 5893 A Preliminary work was done using a conventional Bausch 


where the sodium D lines are located (7). The combined cor- and Lomb (6) spectroscope, with a vibrating mirror and 
rection is obtained from photocell system substituted for the eyepiece. Since this ar- 
Ze on rangement was not portable and was difficult to keep in align- 
ment, a compact instrument has been designed and is being 


AT = 0.0001177,,}-8 


constructed. 


Acknowledgments An assembly drawing of the proposed instrument is shown 
in Fig. 7. Light path and general operation have been de- 
scribed previously. The spectroscope design has been 
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beams regardless of prism position. Light is focused by lens C 
on the entrance slit D. Light leaving the entrance slit is 


The authors wish to acknowledge the valued assistance of 
R. A. Heckman, who initiated the project; also J. G. Bolger, 
R. N. Clark and A. C. Griffin for assistance with the mechani- 
-al design, the electrical construction and the mechanical con- 
struction, respectively. 

This work was sponsored by the U. 8. Atomic Energy Com- 
mission. 


Appendix A. Electronic Balancing Circuit ee 


The over-all operation of the electronic part of the in- 
strument has been described, and is shown in the block dia- 
gram of Fig. 2. Circuit details of the unique components in 
the system are shown in Fig. 6. 

The signal output is generated by an oscillating mirror 
which sweeps a small portion of the spectrum, including the 
sodium D lines at 5890 and 5896 A across the slit in front of the 
photomultiplier. The mirror power supply shown in Fig. 6(a) 
produces a half-wave 60-cycle component to the mirror mag- 
net. This causes the mirror to vibrate at a 60-cycle rate, 
producing a 120-cycle component of the sodium lines in the 
photomultiplier output. Amplitude of the mirror is ad- 
justed to sweep a range from (5887 to 5899 A to produce a 
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Fig. 7 Scanning spectroscope 
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made parallel by lens E. After passing through the prism B 
it is reflected from mirror A through lens F. Light from lens 
F is reflected from the vibrating mirror G and is focused on 
the exit slit H. Normally it passes through a hole in slide I 
into a photocell in box J. For spectrum and slit adjustment, 
a mirror on slide I can be pulled into position to reflect light 
into an eyepiece. 

Since construction of this improved instrument has not 
been completed, its operation has not been tested. 
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UMEROUS papers have been published in the past on 

the determination of equilibrium compositions for react- 
ing gases. A comprehensive but brief review on this subject 
can be found in (1).2 The governing equations, upon which 
the solutions of equilibrium compositions may be obtained are 
mainly of two kinds—the material balance equations and the 
equilibrium equations. Because of the nonlinear algebraic 
form of the latter equations, solutions of closed form are not 
obtained. 

In (1) the various approximate methods were classified into 
four categories: (a) Trial and error methods, (b) iterative 
methods, (c) graphical methods and use of published tables 
and (d) punched-card or machine methods. 

The iterative method presented here is somewhat different 
in nature from those previously published in that category in 
that the efficiency of iteration should be independent of 
whether the mixture is far from or close to the stoichiometric 
proportions. Furthermore, this method has the merit of 
providing explicit correction terms for each succeeding itera- 
tion; that is, the solutions of each iteration may be obtained 
by simple algebraic operations whose forms remain invariant 
for all orders of iterations. Convergence criterion for this 
method may be established mathematically for a simple 
problem of dissociation of diatomic gases which is employed 
mainly as a testing case for convergence. It indicates that 
the results for this case will converge even when the zeroth 
order of approximation involves an error of one hundred per 


cent in magnitude. = 


Analysis 


Introduction 


The equilibrium compositions of a mixture of perfect gases 
having chemical reactions are dictated by Dalton’s law and 
the conditions of equilibrium. For simplicity, we shall take 
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the case of stoichiometric reaction of water vapor as an il- 
lustration, which is represented by the following equations 


Dalton’s law 
pu,o + Pu, + Po. + Pou + Po + PH = P....... 
Conservation of atomic elements +s =>. = 


pu,0 + 2p0, + pou + po — B(2pu,0 + por + 
2pu, + pu) = 0.... [2] 


where @ is the ratio of (Vo/Nx). 


Equilibrium conditions 


Po(po,) 


pu, (po,)'/*(pu,o)~! = Ky,O.............. [5] 


In these equations, all the equilibrium constants are known 
functions of temperature. For a given set of values of pres- 
sure and temperature, this system is sufficient for the deter- 
mination of the six unknown partial pressures. However, be- 
cause of nonlinear algebraic equations involved, solutions of 
closed form have not been obtained. A method based on 
successive approximation has been given by R. Edse (2). 
Various other methods may be found in (1). 

In the following, an alternative method is suggested which 
has the merit of providing correction terms for the next 
iteration. Convergence criterion of this iteration process may 
be established. 

In the above algebraic system, Equations [1-6], the first 
two equations are linear while the other four are nonlinear. 
However, the last four equations are of such a form that their 
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logarithmic differentials are related by linear equations; 
that is 


(dpo/po) — (4dpo,/po,) = 0........... [7] 
(dpu/pu) — (4dpu,/pu,) = 0............. 
(dpu2/puz) + (4dpo,/po,) — (dpu,o/pu,o) = 0... 
(dpon/pou) + (4dpu,/pu,) — (dpu,o/pu,o) = 9. . [10] 


For given values of p and T, the values of po, and pu,o are 
first chosen and they are denoted by 


= 


Then the values of po, pu, Pou, pu, can be readily ob- 
tained from Equations [3, 4, 5, 6]. Now let us denote the 
first approximation to be given by 


= pO + [11] 


where the subscript 7 denotes any one of the components. 
V'rom Equations [1 and 2], we have 


puo® + 2po, + + po + dpu,o + 26po0, + 
dpon + bpo. ... {13] 


— B(2pu,.0 + + + pu + + 
Spon + 26pu, + = 0 
In terms of the difference form of zeroth approximation, 
Equations [7-10] may be written as 
(6po/po™) — (46po0,/po,) = 0.......-. [14] 
(6px /pu™) — (46pn,/pu,) = 0......... [15] 
(6px, + (46p02/po,) — = 0.. [16] 
(pon /pon™) + — = 0. . [17] 
Now, we can eliminate all 6p; except dpo, and épx,0 
from Equations [12, 13] by using [14, 15, 16, 17]. The 
resulting equations from [12 and 13] read 
( (pa.o + 4pon + + pu.) 
+ (po, — $pu, + ipon™ + — X 
(6p0,/po,) = p — Lip. [18] 
[(1 — + 4(1 — 28)pon — 26pu,° — 
+ [2po2 + — 28)pou + 
+ Bou, + Gpo,/po,) = + 
— (1 — — 2p0, 
— (1 — 28)pon — po)... . [19] 


These two equations may be expressed as 
+ = p — [20] 


+ DSpo, = [21] 

from which the quantities dpx,o and 6po, can be evaluated. 
In turn, we can obtain the first approximations pu,o‘ and 
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po,” from Equation [11]. Using Equations [3, 4, 5, 6], the 
first approximations for all partial pressures are determined. 

Succeeding iterations are to be performed in exactly the 
same fashion. The starting equations are 


+ = p — LUpi™....... [22] 
+ D™bpo,” = [23] 
where 


A” = + + + 
Bo = (po, + ipon™ + tpu™)/po,™ 


= [(1 — 28)pxo” + 40 — 28)pou™ — 
2Bpn, — 38pu 


D™ = [2po,™ + — 28)por™ + 4p0™ + 
Bpu,” + ]/po,™ 


ave 


B® = 28pm + — (1 


> 2po, — (1 — 28)pon™ — po™] 
The process is always to determine first 6px,0% and 6po,” 
from the known values of nth approximations of all the > 
partial pressures. Knowing and épo,™, the values of 
(n + 1)-th approximations for other partial pressures are then | 
obtained by using Equations [3, 4, 5, 6]. It is obvious that — 
iteration may be terminated as soon as the right-hand sides of 
Equations [22 and 23] become insignificantly small. 

Generalization of this method to other chemically reacting — 
systems may be easily made in a straightforward fashion as 
soon as the number of independent chemical reactions in the 
system can be ascertained and the governing equations of 
equilibrium are formulated. 

In the following, we will study the convergence criterion of 
the above iteration scheme. The method is applied to the 
simple case of dissociation of diatomic gases whose exact 
solution of equilibrium composition can easily be determined. 
The reason for choosing this particular example is the simplic- 
ity in establishing the convergence criterion for this case. 
Similar criteria may be established for more complex cases. 
However, they will require tedious algebraic manipulation. 

The governing equilibrium equations for the dissociation of 
diatomic gases, (1/2)Y2 =— Y, are 


In applying the present iteration method to this case, we 

write 

+ Spy,” = p — py — py2......... [26] 


— = 0. . [27] 


First we choose a value of py,; then the value of py is ob- 
tained from Equation [25], i.e. 


Substituting this into [26 and 27], we have 


Spy + = p — K(py,)'? — py,...... [29] 
Spy — 4K(py,)'Spy2 =0........... [30] 
253 
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Eliminating dpy from [29 and 30], we obtain 
[1 + 4K(py,%) = p — K(py,)'* — py,®.. [31] 


Since py, = py, + dpy,, we can rewrite Equation [31] 


as 
= 1+ py, — KV + 

1K V py2)]... . [32] 
By induction, the following general relation may be derived 


= 1+ [33] 


py,'?/py.) 


where 
K 


In Equation [32], we see that if py2 is estimated to be 
larger than the correct answer, the second term on the right- 
hand side of that equation will be negative; thus py,“ will be 
smaller than py2. On the other hand, if py, is estimated 
to be smaller than the correct answer, the second term on the 
right-hand side of that equation will be positive; thus py,“ 
will be larger than py,. The same characteristic is also 
exhibited by the general equation [33]. Therefore, qualita- 
tively we see the trend of this iteration scheme proceeds in the 
right direction. Furthermore, it is seen that the iteration will 
be automatically terminated as soon as the second term on 
the right-hand side of Equation [33] approaches zero. 

It follows from this discussion that the convergence of this 
iteration method can be proved if the following two conditions 
are shown to be valid 


To prove condition [34] is equivalent to requiring 
| 
| (n) (n) (n +1) (n +1) 
+ VPM Pp — py2 


or 
| py," (p py,” KV 
| (py. (p KV py,"*»)| 


By using Equation [33], the term on the left-hand side of the 
above inequality, after rearranging, may be writtenas 


K 
x (n) (x += + 


Similarly, the term on the right-hand side of the inequality 
may be written as = 


| 
| py. {( ps,» + Qm(1 + 
i 


Examining these expressions closely, we see that the above 
inequality holds provided that 


To prove condition [35] is equivalent to showing 


because the denominator of Q“*» is always positive. Now 
we will express the above quantity in terms of Q™ as follows 


KV + QM, 


The last term on the right-hand side can be expanded by 
binomial series when |Q™| < 1 


py2"t) K V = el + + 
KV pve [1 + Q™ — 


As a consequence of [34], the absolute value of Q™ will 
successively decrease as n increases, higher order terms in- 
volving square of Q™ or higher may be dropped from the 
above expression when n approaches large integer V. That is 


+ K V > py2™ (1 +Q™) + 
KV py (1 + 4Q) =p 


Therefore, [35] is verified under the proviso of [36]. 
Consequently, the condition |Q™|< 1 may be taken as the 

criterion of convergence of this iterative method for the case 

of dissociation of diatomic gases. In terms of the initial 


choice, zeroth order approximation, this condition reads 


@ 


lp — py, — KV py, - 
K 
| py,” V py, 


which is not at all stringent. For example, the zeroth order 
estimation involving +100 per cent error will still satisfy this 
condition. 

It should be emphasized here that the study for the case 
dissociation of diatomic gases is mainly for the purpose of 
illustration since exact solution of the equilibrium composition 
for this case can be explicitly obtained. Furthermore, the 
rate of convergence is not established. Numerical computa- 
tions, however, indicate fast convergence; first or second order 
approximations will provide accurate results in most cases. 

For many cases where the mixture is far from the stoichio- 
metric proportions, judicious selection of the zeroth order 
approximations may be made along the line suggested in (2) 
which will shorten the number of required iterations. 

In view of the invariance of the form of the iteration 
formula for a given reacting system, this method should be 
particularly well adaptable to automatic machine com- 


putation, 
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Technical Notes 


Ignition of Electrolytic Monopropellants 
by Submerged Electrical Discharge’ 


MARJORIE W. EVANS,? FRANK I. GIVEN? and G. M. 
MULLER* 


Stanford Research Institute, Menlo Park, Calif. 


LECTROLYTIC liquid monopropellants can be ignited 
in the liquid phase by submerged high voltage electrical 
discharge. The authors have studied the nature of high volt- 
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age discharge in electrolytic liquids, the factors controlling — 
the ignition or failure to ignite of monopropellants, and the | 
propagation rate and shape of the reaction front following 
ignition. Experimental techniques included measurements of — 
current through, and potential drop across, a submerged gap _ 
as a function of time, with a precision of +10 per cent. Simul-— 
taneously, a modified Beckman and Whitley Model 189 Fram- — 
ing Camera was used to obtain time sequence pictures at a— 
rate of 0.25 frame/ysec. 

Upon sudden (rise time = 5 usec) application of 1.5 to 4 _ 
kilovolts from a storage condenser to a localized electrode 
submerged in the monopropellant, a corona of finite lifetime 
is formed in the neighborhood of the electrode. This corona 
grows and then fades, disappearing at 30 to 40 usec. If the 
anode-cathode geometry is such that the corona dies before a 
continuous path is formed through the corona from anode 
to cathode, no arcing occurs, and the remaining portion of the 
energy on the condenser may be discharged through the elec- 
trolytic liquid, the current being carried by the ions in elec- 
trolytic conduction. If, however, a continuous plasma path 


Fig. 1 Corona discharge and ignition in hydrazine, hydrazine nitrate monopropellant—voltage, 3500 v; condenser capaci- 
tance, 55 uf; framing rate, 0.25 frame/usec; exposure time, 1.3 usec 


Eprror’s Note: The Technical Notes and Technical Comments sections of JET PROPULSION are open to short manuscripts describing new 
developments or offering comments on papers previously published. Such manuscripts are published without editorial review, usually 
within two months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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Fig. 2. Cross section reaction vessel; width of chamber = 5.1 
cm; all dimensions are in centimeters 


is effected during the corona lifetime, arcing occurs and the 
condenser rapidly discharges through the low resistance path 
provided by the are. 

For proper conditions of geometry and voltage, the 
monopropellants can be ignited either by are or corona dis- 
charge, and the initial reaction occurs within the region af- 
fected by the are or by the corona. A study of corona ignition 
for a variety of anode-cathode geometries and voltages 
showed that ignition is a function of energy density within 
the corona. By suitably controlling the geometry and volt- 
age, ignition can occur upon delivery to a corona of as little 
as 1 to 10 joules. Criteria for effecting ignition with the low- 
est total expenditure of energy from the storage condenser 
were tentatively established. 

Upon successful ignition, a self-propagating reaction zone 
moves through the liquid at the rate of 70 to 100 meters per 
sec for as long a time as the vessel walls remain intact and 
hold the liquid together. For a hemispherically growing reac- 
tion zone, this propagation velocity means that, within 100 
usec of the initial application of voltage to the electrode, reac- 
tion begins in approximately 2 ce of liquid; and, within 250 
usec, 8 ce of liquid have undergone at least partial reaction. 
Assuming very fast chemical reaction, the 65-30-5 (per cent 
weight) hydrazine, hydrazine nitrate, water propellant used in 
these experiments will produce in 250 usec as much as 0.5 
mole of gaseous products. 

A series of framing camera pictures of monopropellant ig- 
nited by a corona discharge is shown in Fig. 1. A cross section 
of the reaction chamber is given in Fig. 2. Plexiglas windows 
1 in. thick bolted to a heavy metal vessel formed two sides of 
the chamber. The protruding finger from the top of the ves- 
sel in Fig. 1 was floating at electrolyte potential and had no 
significance in the experiment. The vessel was completely 
filled with monopropellant. The first seven frames show the 
growth and decay of the corona. Ignition occurred within the 
corona during that time (28 ywsec). The remaining frames 
show the self-propagation of the reaction zone into the liquid 
for 36 wsec more. According to other data taken simultane- 
ously, the reaction zone propagated for a total of 230 usec, at 
which time the Plexiglas windows shattered under the pres- 
sure of the products of reaction. Ignition for this particular 
geometry (not by any means the optimum geometry for 
economy of energy) required delivery to the corona-electrolyte 
system of 115 joules, of which 35 were delivered to the corona 
itself and were instrumental in igniting the monopropellant. 


function played by anode-cathode geometry and condenser 
voltage in ignitibility of electrolytic monopropellants is to be 
presented at the Seventh Combustion Symposium, Oxford, 


_A Note on Sonic Sudden Enlargements 


R. H. PAGE! 


Stevens Institute of Technology, Hoboken, N. J. 


Introduction 


IGURE 1 has been found quite useful in estimating pres- 

sure ratios across choked abrupt area changes followed 
by short tubes. Since the permissible range of pressure ratios 
across a sonic sudden enlargement can be ascertained from 
this chart, problem solutions may often be expedited by its 
use. 


Use of Fig. 1 


The geometry considered in this paper and a typical pres- 
sure profile are shown in the figure. How the acoustic state at 
section 1 is approached (e.g., Fanno line or isentropic ex- 
pansion) does not influence these results. For a given area 
ratio A»/A;, the maximum value (curve A), probable mini- 
mum value (curve B), and absolute minimum value (curve C), 
of the static pressure ratio P,/P, are plotted. Curve B is 
called the probable minimum because it is more likely to occur 
than curve C. The curve C solution is somewhat analogous 
to the strong shock solution for angle shock since it must be 
forced on the flow to apply. (However, curve B represents a 
larger entropy change across the sudden expansion than curve 


The small charts at the top of the figure present the stagna- 
tion pressure ratio, Po,/Po, (maximum and minimum values), 
and the length of the “recovery zone” L/D» (where Dz is the 
downstream tube diameter). Since the length of the recovery 
zone is dependent on a number of factors, the dashed curve 
only represents an approximation. 

The maximum values, curve A, are independent of the 
geometry of the expansion and may be used for both outward- 
expanding and inward-expanding axially symmetric and two- 
dimensional sonic sudden enlargement. The minimum values, 
curves B and C, are, strictly speaking, only applicable to out- 
ward expanding axially symmetric configurations (the most 
frequently encountered geometry). However, minimum 
values calculated for two-dimensional sudden enlargements 
and sudden enlargements around abrupt central bodies in 
circular conduits fall within 10 per cent of these curves. Thus, 
Fig. 1 may be used, with some reservations, for configurations 
other than axially symmetric ones. 

The top chart emphasizes that this figure is restricted to ex- 
pansions followed by a short tube or to the short portion of a 
longer tube. The over-all pressure ratio for a sonic sudden en- 
largement followed by a rather long tube includes the down- 
stream pressure ratio P;/P2 which must be calculated from the 
downstream flow path (e.g., Fanno or Rayleigh line). For 
long tubes acoustically choked at section 1, the over-all pres- 
sure ratio P;/P, (where P; is the pressure at tube exit) may 
become rather low, limited by section 3 becoming sonic. For 
that case: P;/P; = Po,/Po, = A:/Az. If the tube is very long, 
a double throat problem may exist, and it will be necessary to 
check (1)? to see if acoustic flow will ever be reached at section 
1. For, if in the transient starting operation, section 3 (the 
end of the pipe) chokes before section 1, the sudden enlarge- 
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A detailed description of the behavior of electrical dis- Brofesor of Mechanical 


charges in electrolytic monopropellants and an analysis of the 
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Fig. 1 Pressure ratios for axially symmetric flow 


ment will remain subsonic. (Subsonic sudden enlargements 
are solved by well-known procedures (2, 3) for which con- 
venient charts are available (4).) 


Origin of Curves 


A modified one-dimensional adiabatic analysis in which the 
shear stress at the wall is neglected has been used for the cal- 
culations. In applying the momentum principle, the force 
acting on the annular area or “base area’? (A»-A;) is 
puted from the pressure acting on the base area. This “base 
pressure” is treated two ways: 

1 For curve A, the base pressure was set equal to the sonic _ 
pressure P; and thus, in the manner of previous investigators — 
(2, 3), a maximum pressure regain curve was calculated. 
This approaches unity asymptotically as the area ratio A2/A;_ 
is increased. Of course large values of A»/A, invalidate the — 
one-dimensional assumption and, therefore, the curve is not _ 
continued for areas ratios larger than three. 

2 For curves B and C, the minimum pressure ratios were 
desired. Low pressure ratios P2/P, are obtained (2, 5) when 
the pressure acting on the annular area (base pressure) be- 
comes less than the sonic pressure P;. Indeed, the lowest 
pressure ratios, P2/P,, occur when the base pressure is at a 
minimum. Lowest limiting base pressure values are available 
(6) and were used in the analysis to calculate curves B and C. 
In this way the rather complex axially symmetric analysis (6) 
was merged into a one-dimensional treatment. 

The recovery zone length is not subject to a simple defini- 
tion and depends upon parameters not included in this note. 
This provisional curve was constructed as a compendium of 
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experimental observations (2,7,8). It is not accurate, but is 
intended to present the order of magnitude of the recovery 
length. 


Further Details 


For most engineering purposes the maximum and minimum 
values of the pressure ratio P2/P, will provide sufficient design 
information. The exact determination of the pressure ratio 
P./P,, rather than the determination of maximum and mini- 
mum values, is much too involved to be considered in this 
note. However, detailed studies have been conducted (6, 9). 
W. L. Chow in a private communication indicates that 
further investigations are being carried out and that a com- 
prehensive report may be forthcoming. 


References 


1 Shapiro, A. H., “The Dynamics and Thermodynamics of 
Compressible Fluid Flow,’’ vol. 1, Ronald Press, New York, 
1953, pp. 143-203. 

2 Cole, B. N., and Mills, B., “The Theory of Sudden Enlarge- 
ments Applied to the Poppet Exhaust-Valve, with Special 
Reference to Exhaust-Pulse Scavenging,’ Proceedings of the In- 
stitution of Mechanical Engineers, vol. 1B, no. 8, 1952-1953, pp. 
364-377. 

3 Hall, W. B., and Orme, E. M., ‘Flow of a Compressible 
Fluid Through a Sudden Enlargement in a Pipe,’ Proceedings of 
the Institution of Mechanical Engineers, vol. 169, no. 49, 1955, pp. 
1007-1015. 

4 Unterberg, W., “Communication on ‘Flow of a Compressible 
Fluid Through a Sudden Enlargement in a Pipe’ by Hall, W. B., 
and Orme, E. M.,”’ Proceedings of the Institution of Mechanical 
Engineers, vol. 169, no. 49, 1955, p. 1018. 

5 Wick, R.S., “The Effect of Boundary Layer on Sonic Flow 
Through an Abrupt Cross-Sectional Area Change,” Journal of 
the Aeronautical Sciences, vol. 20, no. 10, 1953, pp. 675-682. 
“Further Comments on the Effect of Boundary Layer on Sonic 
Flow Through an Abrupt Cross-Sectional Area Change,” Journal 
of the Aeronautical Sciences, vol. 22, 1955, pp. 135-137. 

6 Korst, H. H., Page, R. H., and Childs, M. E., ““A Theory for 
Base Pressures in Transonic and Supersonic Flew,” University of 
Illinois, Engineering Experiment Station, ME-TN-392-2, OSR- 
TN-55-89, 1955. 

7 Hall, W. B., and Orme, E. M., ‘Authors’ Reply—Com- 
munication on ‘Flow of a Compressible Fluid Through a Sudden 
Enlargement ina Pipe’,”’ Proceedings of the institution of Mechan:- 
cal Engineers, vol. 169, no. 49, 1955, p. 1020. 

8 Keochendorfer, F. D., and Rousso, M. D., ‘Performance 
Characteristics of Aircraft Cooling Ejecters Having Short Cylin- 
drical Shrouds,”” NACA RM E51E01, 1951. 

9 Korst, H. H., Chow, W. L., Zumwalt, G. W., and Smith, 
S. J., ‘Some Investigations of the Effect of Approaching Flow 
Boundary Layer on the Base Pressure Solution for Axi-Sym- 
metric Sudden Expansion,’’ University of Illinois, Internal Re- 
port, Mechanical Engineering Department, 1956. 


> 
> 


Reflections in Two-Dimensional 
_ Diffusers From Hydraulic Analogy 
Experiments 


E. V. LAITONE! and J. E. STOUT? 


The separation distance required to prevent the forma- 
tion of a Mach reflection between two wedges, representing 
either a two-dimensional diffuser or a Busemann biplane, 
is obtained by assuming that the normal shock wave could 


not have been created if the first expansion wave from the 
shoulder of the wedge would intersect it. This hypothesis 
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is verified experimentally for the special case corresponding 
to a specific heat ratio of y = 2 by means of the shallow 
water hydraulic analogy. 


IGURE 1 gives an indication of the range of validity of the 
shallow water flow hydraulic analogy (see (1, 2 or 3)*) to 
two-dimensional gas dynamics for the particular cases occur- 
ring when either determining the minimum Mach number 
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Fig. 1 Minimum free stream Mach number for a given stream 


tube deflection, N = 1, or for a reflection from a plane wall, 
N=2 
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Fig. = 1.68, H/C = 2.41 


that must be attained in order to prevent a detached norma! 
shock wave from forming in front of a given stream-tube de- 
flection of 6.(N = 1), or when determining the Mach num- 
ber that must be exceeded in order to prevent the formation 
of a Mach reflection (V = 2). The usefulness of the direct 
analogy to two-dimensional compressible flow is limited to 
the regime wherein the oblique shock wave relations are 
approximated to the desired degree of accuracy by the isen- 
tropic compression values. As pointed out in (1 and 2), no 
hydraulic analogy exists for the mixed problem of transonic 
flow (having both subsonic and supersonic domains inter- 
acting), or whenever isentropic domains interact with the 
shock wave, unless the entropy increase through the shock 
wave can be neglected. That is, the direct analogy fails 
whenever the exact relations depicted in Fig. 1 are not suffi- 
ciently accurate approximations to one another. 

It is obvious from Fig. 1 that the model laws determined 
by the transonic similarity approximation, which are gener- 
ally used to convert hydraulic analogy experimental results 
(3) from a specific heat ratio equivalent to y = 2 to the usual 
gas dynamic case of y = 7/5, cannot have as large a range of 
validity as the hydraulic analogy itself possesses in checking 
the fundamental equations of motion with y = 2. Conse- 
quently, as pointed out in (1 and 2), the hydraulic analogy is 
more limited when used as a direct model test than it is when 
used as a check on any theoretical derivations or hypotheses 
concerning the behavior of the fundamental! nonlinear equa- 
tions of motion. Even when the hydraulic analogy to two- 
dimensional compressible flow completely breaks down, the 
shallow water tests are still useful for checking various meth- 
ods used to develop approximate mathematical solutions of 
the boundary value problems, similar to the type occurring 
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Fig. 6 Minimum separation distance which prevents Mach 
reflection between two wedge profiles 


n gas dynamics, arising from the nonlinear equations de- 
fining the hydraulic jump and the potential shallow water 
theory. The only restrictions then are that the vertical 
iccelerations outside the hydraulic jumps be negligible, and 
that all pressures be approximately hydrostatic. As shown in 
(1 and 2), these restrictions are satisfied, at least at transonic 
conditions, if the water depth everywhere remains less than 
1/,-in. and the model size is suitably large. 

As an example of the application of shallow water experi- 
ments in formulating a theoretical solution valid for the gas 
dynamics case, even though the hydraulic analogy itself would 
fail as a direct model test, we shall determine a criterion for 
the necessary separation distance between two wedge profiles 
(simulating either a Busemann biplane or a two-dimensional 
diffuser) that would prevent the formation of a normal shock 
wave in between them; that is, the separation distance that 
would eliminate a Mach reflection. 

Figs. 2 and 3 show typical cases of a normal shock wave 
producing a choking in the internal flow between two wedges. 
The photographs were taken by towing two 10-in. chord, 
10-deg nose angle wedges at constant velocity through still 
water at a depth of 0.2 in. in a tank 20 ft long and 4 ft wide 
as described in (1). Fig. 2 shows the normal shock wave in 
front of both wedges since the free stream or towing Mach 
number is 1.50 which is less than the minimum value of 1.51 
corresponding to a detached hydraulic jump for a 10 deg nose 
angle in Fig. 1 (V = 1). 

In Figs. 3, 4 and 5, the free stream or towing Mach number 
is 1.68, which is greater than that defined by N = 1 in Fig. 1, 
but is still less than that given by N = 2. Therefore, a Mach 
reflection occurs whenever the attached nose shock waves 
intersect before they are met by the expansion waves arising 
from the wedge shoulders. That is, as the separation dis- 
tance is increased, the expansion waves from the shoulder 
must eventually cancel the normal wave as in Fig. 5. Since 
the flow is everywhere supersonic, the nose shock waves form 
‘straight lines; consequently the minimum separation distance 
Hi, in order to have the first expansion wave arising at the 
‘shoulder meet the shock wave intersection point as shown in 
Fig. 6, is easily evaluated as 


H |= + are sin 1/M.) — tan | 
= tan 0 


Cc tan (0, + are sin 1/M.2) — tan 0, 
ein 
(= 


C 
[cos — 0) — — 1sin (6, — 


Fig. 6 presents Equation [1] for both the shallow water 
and the ideal gas flow having y = 7/5. The experi- 
mental data given in Fig. 6 were obtained from the hydraulic 


-analogy tests, the minimum separation distance which would 
-eliminate the normal wave or hydraulic jump being deter- 


mined from photographs similar to Figs. 4 and 5. The data 
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indicate that a separation distance only slightly greater than 
that given by Equation [1] for shallow water theory will pre- 
vent the occurrence of a normal wave. It would now be of 
great interest to compare Equation [1] with corresponding 
experimental data in gas dynamics. Such tests would deter- 
mine the range of validity of Equation [1], and, at the same 
time, should verify the similarity in the physical behavior 
between a shallow water flow and a two-dimensional com- 
pressible gas flow, even though the hydraulic analogy itself 
is not directly applicable as a model test for the entire range 
of the experimental data shown in Fig. 6. 
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Friction and Heat Transfer in a Rough 
Tube at Varying Prandtl Numbers 


R. C. HASTRUP,! R. H. SABERSKY,? D. R. BARTZ and 


M. B. NOEL4 


- California Institute of Technology, Pasadena, Calif. — 


A set of experiments was carried out with the purpose of 
extending the available data on heat transfer in rough 
tubes as well as the information concerning the relation- 
ship between heat transfer and friction. The results are 
presented in terms of the heat transfer coefficient C, and 
the friction coefficient C;. Particular attention was paid 
to the effect of Prandtl number. The experiments were 
carried out in an electrically heated tube with water as a 
working fluid. By using water at various temperatures ex- 
periments could be performed over a continuous range of 
Prandtl numbers from about 1 to7. For the tube and the 
type of roughness used in this experiment the results con- 
firm that roughness increases the heat transfer as well 
as the friction coefficient, but that the ratio C,/C, de- 
creases with increasing roughness. The latter effect 
seemed slightly more pronounced at the lower Prandtl 


numbers. 
Nomenclature 
c = specific heat 
Cy, = friction coefficient 


C;, = heat transfer coefficient 


d = tube diameter 

de = equivalent diameter for rough tubes 
g = gravitational constant 

k = conductivity of tube wall 

L = length of tube section 

p = pressure 

q = heat transfer rate per unit area 
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Fig. 1 Schematic diagram of test installation 
Re = Reynolds number = Vd/y 
s = thickness of tube wall 

t = temperature 

t, = bulk temperature of fluid 


ty; = temperature at inside of tube wall Seopa 
two = temperature at outside of tube wall oo 
u = local velocity 
V = average fluid velocity 
ey = turbulent exchange coefficient for friction jr 
€, = turbulent exchange coefficient for heat transfer a 
v = kinematic viscosity pay 
p = density ; 
o = Prandtl number = cu/k 
7 = shear per unit area a 

Introduction 


HE effect of wall roughness on friction in circular pipes 

has been well known fora long time. The most extensive 
set of experiments is perhaps that due to Nikuradse (1).5 In 
these experiments Nikuradse prepared pipes of different de- 
grees of roughness by means of sand grains of various sizes. 
For each pipe the friction drop was measured and the results 
were presented in a graph of the friction factor vs. Reynolds 
number. The experiments covered a range of relative rough- 
ness up to 1/;5, where the relative roughness was defined as 
the ratio of the sand grain diameter to the pipe diameter. 

A similar set of experimental data for the heat transfer 
coefficient in tubes of different roughness would be equally 
useful. It would be of theoretical interest as well as of prac- 
tical importance, for example, in the design of heat exchanger 
surfaces. In contrast to the rather complete data available 
on the friction coefficient, the information on the heat transfer 
coefficient in rough tubes is still relatively scarce. Of course, 
some work has been done in this field—one of the most valu- 
able studies perhaps being the work recently published by 
Nunner (2). 


* Numbers in parentheses indicate References at end of paper. 
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The lack of information becomes particularly apparent 
when it is realized that much more data are required for a 
proper representation of the heat transfer coefficient than for 
the representation of the friction coefficient. This results 
from the fact that the friction coefficient is determined ade- 
quately by two parameters: the relative roughness and the 
Reynolds number. At least one additional parameter, the 
Prandtl number, is, however, required for the proper deter- 
mination of the heat transfer coefficient. For a complete 
investigation of heat transfer in rough tubes it is, therefore, 
not sufficient to represent the heat transfer coefficient as a 
function of Reynolds number for various relative roughnesses, 
in analogy with Nikuradse’s graph. A graph of this type 
will have to be obtained for various Prandtl numbers. It was 
the purpose of these experiments to contribute to the infor- 
mation necessary for a complete survey of the heat transfer 
in rough tubes. The phase which was selected for this in- 
vestigation was the effect of Prandtl number on the heat 
transfer coefficient in a tube of given roughness. 


Test Equipment and Procedure 


In order to conduct the desired experiments it was decided 
to use water as the working fluid. Water has the character- 
istics convenient for this experiment; that is, its Prandtl 
number is markedly affected by temperature, and its density 
is essentially constant. At room temperature the Prandtl 
number is approximately 7, and it can be made to decrease to 
unity by raising the temperature to about 300 F. By chang- 
ing the temperature, therefore, a continuous variation of the 
Prandtl number could be effected without having to change 
fluids. 

The experimental installation is shown schematically in 
Fig. 1. The water is contained in a reservoir and heated 
under pressure to the desired temperature; it is then forced 
out of the reservoir by pressure from a nitrogen supply tank. 
The water then flows through the test section into a receiver 
tank. The pressurized system produces a very steady flow 
and allows accurate flow and pressure measurements. The 
tank capacity was sufficiently large for run durations of the 
order of 5 min. This time interval was adequate for estab- 
lishing equilibrium operating conditions. 

The test section itself consisted of a stainless steel tube 
(Type 321) of */s in. diam. The wall thickness of the tube 
was 0.010 in. The length of the complete tube was 39 in. 
Of this, only 18 in. at the downstream end were heated, and 
this portion represented the actual test section. The up- 
stream portion of the tube, corresponding to about 50 diam- 
eters, served as entrance length to establish an equilibrium 
velocity profile. The entire test section together with about 
15 diameters of the entrance length were artificially rough- 
ened by knurling from the outside. During the knurling 
operation the inside of the tube was filled with a low melting 
point alloy (‘Cerro bend’’) which was melted out after com- 
pletion of the operation. The knurl pattern used was dia- 
mond shaped and relatively coarse. The diagonals of the 
diamond were approximately 0.10 and 0.05 in., and the depth 
of impressions produced on the inside of the tube was approxi- 
mately 0.008 in. Immediately upstream and downstream of 
the 18-in. test section, a small length of tube was left smooth. 
These smooth regions were provided so that static pressure 
taps could be located at points where there would be no local 
disturbances due to surface roughness. 

In addition to the rough tube a smooth test section with a 
smooth entrance section was also procured. The smooth test 
section was used because it seemed desirable to conduct ex- 
periments with such a tube first, so that the results could be 
compared with those of other experimenters as well as with 
the results of the rough tube. 

Heat transfer to the 18-in. section of both the rough and the 
smooth tube was provided by electrical means. Alternating 
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current was passed through the stainless steel tube, the thin 
wall of the tube itself serving as a resistance heater. Electric 
heating has proved to be quite convenient in experiments of 
this type, as it allows rapid adjustment and accurate deter- 
mination of heat flow rates. 

The principal measurements that were made consisted of 
(a) the temperature of the water entering the test section, (b) 
the surface temperature at two places on the outside of the 
test section, 9 and 16 in. from the entrance, (c) the pressure 
drop across the test section and (d) the electric power dissipa- 
tion in the test section. In addition the flow rate was deter- 
mined by means of a Venturi meter. The dimensions of the 
test section were, of course, also required. In the case of 
the rough tube an effective tube diameter was defined on the 
basis of the volume of the tube as 


d. = 


where vol is the volume of the tube and L the corresponding 
length. The volume was determined by measuring the 
amount of water required to fill the tube. 

The quantities which we sought in these experiments 
were the friction coefficient C; and the heat transfer coeffi- 


cient C,. These coefficients are given by Sw 
dp d Bt. 


— th) 


where ¢,.; is the temperature on the inside surface of the tube 
wall and ¢, is the bulk temperature of the fluid. Both C; and 
C, may be computed from the aforementioned measure- 
ments. In computing C;, the pressure gradient dp/dz was 
assumed to be constant along the test section. It was there- 
fore set equal to Ap/L where Ap is the pressure drop across 
the test section and L is its length. The velocity V was ob- 
tained from the volume flow rate and the cross-sectional area 
of the tube. In the case of the rough tube the diameter was 
taken to be equal to the effective diameter and the cross- 
sectional area was also based on this diameter. The heat 
transfer rate per unit area in Equation [3] was obtained as- 
suming uniform power consumption along the tube. This 
assumption was justified as the specific electrical resistance 
of the tube was practically constant for the temperature 
changes which occurred in the tests. The inside wall tem- 
perature ¢,; was calculated from the temperature measure- 
ment at the outside of the wall t..., correcting for the temper- 
ature drop through the wall by the expression 


qs {1 8 
tes = — — — 4 
k (; 


The bulk temperature at any point along the tube was com- 
puted by adding to the temperature of the fluid entering the 
test section, the temperature increase caused by the total 
heat transfer to the fluid from the tube wall up to the sec- 
tion in question. With these data and knowing the fluid 
properties, the coefficients C; and C, were computed. It 
may be pointed out that the coefficients thus obtained repre- 
sent local values applying to a particular section of the tube 
and that they are not averages over the length of the tube. 

Great care was taken in making all the measurements. 
Particular attention was given to the problem of temperature 
measurements. In order to obtain satisfactory accuracy it 
was found necessary to calibrate each thermocouple in place. 
With precautions of this type it was possible to obtain repro- 
ducibility of the experimental coefficients C; and C, within 
1 and 2 per cent, respectively. This reproducibility was also 
achieved in cases where a given Reynolds number was ob- 
tained in different ways, i.e., by means of different sets of 
values of velocity and viscosity. 
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Fig. 2. Friction coefficients for the rough and the smooth tube 
as a function of Reynolds number 
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Fig. 3 Heat transfer coefficient for the smooth tube as a func- 
tion of Reynolds number 
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Fig.4 Heat transfer coefficient for the rough tube as a function 


The largest uncertainty in determining the absolute accu- 
racy of the measured coefficients was probably introduced by 
the value for the conductivity of the stainless steel tube. This 
quantity was needed in determining the temperature tui 
as shown by Equation [4]. Values ranging from k = 0.245 
10-* Btu/in.-sec-°F at 100 F to 0.28 10-3 Btu/in.-sec-°F 
at 300 F were used for the present computation. Taking 


10* 


into account possible deviations from this value as well as 
other inaccuracies in the experimental results, it is estimated 
that the friction coefficient C; is accurate within 2 per cent 
and that the heat transfer coefficient C;, is accurate within 


5 per cent. 
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Fig. 5 The ratio 2C;,/Cy as a function of Prandtl number for — - 


Re = 1.5 X 10°; the results corresponding to the Prandtl ana- 


logy, to the von Karman analogy, and to the empirical relation 
aan of Colburn are also shown for comparison 

pee 


Experimental Results 


The experimental results are presented in the form of sev- 
eral graphs. The first graph (Fig. 2) shows the variation of 
the friction coefficient with Reynolds number for the rough as 
well as for the smooth tube. The smooth tube experiments 
are in close agreement with published data. The curve 
C, = 0.046/Re®-2, which fits the published data for smooth 
tubes well, is shown in Fig. 2 for comparison. The rough tube 
shows a marked increase of the friction coefficient over that 
of the smooth tube. It is seen that the friction coefficient of 
the rough tube is still dependent on the Reynolds number. 
This shows that, at least for the range which was investi- 
gated, the tube cannot be considered ‘perfectly rough.” 
For the type of roughness used here, it appears that the Reyn- 
olds number dependence extends to higher values of Reyn- 
olds numbers than for a sand-grain-surface of equal relative 
roughness. In Figs. 3 and 4 the heat transfer coefficient C;, is 
shown as a function of Reynolds number for various constant 
values of Prandtl number. Fig. 3 applies to the smooth 


tube and Fig. 4 to the rough tube. For comparison, the 
equation 


has also been plotted in Fig. 3. This equation is one of sev- 
eral similar ones frequently recommended for computing heat 
transfer in a smooth tube. Agreement between the empirical 
equation and the measured data is considered adequate to 
give confidence in the experimental techniques of the present 
series of tests. The data presented in Figs. 3 and 4 show 
that the roughness improves the heat transfer coefficient 
somewhat over that of the smooth tube. In each case, 
however, the increase in the friction coefficient caused by the 
roughness exceeds the corresponding improvement in the 
heat transfer coefficient. 

It should be pointed out that the coefficient C, in Figs. 3 
and 4 is the isothermal value of C,. The isothermal value 
was obtained experimentally as follows: At each Reynolds 
and Prandtl number the heat transfer coefficient was deter- 
mined for several heat transfer rates g. The coefficient C, 
was then plotted against g, and the curve was extrapolated to 
a zero heat transfer rate. The value of C, at ¢ = 0 corre- 
sponds to a zero temperature difference between the wall and 
the bulk fluid. This is the value of C;, which has been used 
in this presentation. The isothermal value was selected for 
this investigation because it was felt desirable not to veil the 
effect of roughness by the effects which are due to temper- 
ature and property variations across the boundary layer. 
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The effects of roughness on friction and heat transfer can 
perhaps be compared best by forming the ratio 2C,/C,;. This 
ratio has been plotted in Fig. 5 as a function of Prandtl 
number for a given Reynolds number. Two experimental 
curves are shown, one for the smooth tube and one for the 
rough tube. For comparison three theoretical relationships 
are also shown. The first is based on the empirical relation 
for C;, given by Equation [5] and on 


which very closely fits the data for the smooth tube in Fig. 2 
The other two curves are based on the analogies betwee 
heat transfer and friction as analyzed by Prandtl and voi 
K4rmén, respectively. Fig. 5 has been prepared for a par- 
ticular Reynolds number. Graphs for other Reynolds num 
bers within the range of this investigation are quite similar. 


Discussion of Results “6 «ip? 


The principal results of this set of experiments are that 
for the type of roughness tested, (a) the heat transfer coeffi- 
cient C,, increases with roughness and (b) the increase in the 
friction coefficient C; exceeds that of the heat transfer coeffi- 
cient. As a consequence, the ratio 2C;,/C; decreases with 
roughness. These results are in agreement with the findings 
of previous experimenters (see for example (2 and 3)). For 
the rough tube tested, the ratio 2C),/C; is below that for the 
smooth tube for the full range of Prandtl numbers. The 
difference seems to be more pronounced at the lower end of 
the range. In the neighborhood of ¢ = 1 the difference in the 
ratios 2C),/C; for the two tubes is about 15 per cent and it 
decreases to about 10 per cent at o = 7. 

The results may be analyzed further in the light of Reyn- 


olds analogy. Following the usual approach, one may 
write for the shear near the wall of a pipe 
@: T du 
- = (vy + €&) 3 [7] 
+. p y 
The heat transfer rate per unit area may be expressed as 
v dT 
pc o dy 


The quantities ¢, and e&, are the turbulent exchange coeffi- 
cients for friction and heat transfer, respectively. According 
to Reynolds analogy e, is proportional to ¢. In most analyti- 
cal work the proportionality constant is set equal to unity for 
lack of sufficient experimental data. In that case it can be 
shown from the above two equations that for ¢ = 1 the ratio 


= 1 


It may be noted here that the experimental value of 2C,/C; 
obtained for the smooth tube does not equal 1.0 but rather 
0.9 as seen from Fig. 5. The data are of course subject to 
the limitation in accuracy mentioned before. Nevertheless 
one may take this result as an indication that the propor- 
tionality constant between ¢, and ¢ differs from unity. In 
this connection it might also be mentioned that experiments 
with air in smooth tubes (see (2 and 3)) yield values of about 
1.18 for the ratio of 2C,/C;. The Prandtl number for air at 
the temperature in question is approximately 0.70. If the 
coefficient C, is now assumed to be proportional to «~*-§ 
(see Equation [5]), the data for air may be extrapolated to a 
Prandtl number of unity. This extrapolation also leads to a 
value of 2C;,/C, slightly below unity. 

Let us now further examine Fig. 5 in which the ratio 2C,/C;, 
is plotted as a function of Prandtl number for the rough as 
well as for the smooth tube. As is seen, both curves de- 
crease rapidly with Prandtl number. The fact that 2C,/C, 
for the rough tube decreases with increasing Prandtl number 
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is somewhat disappointing. One might have imagined that 
the effect of roughness would be to increase the turbulent ex- 
dt] change coefficients to such an extent that vy < ¢& and v/a < 


tal e, right up to the wall. In that case it would follow from 
‘he Equations [7 and 8] that 

ips 2C, 

on —_ = 


Cy 


If it is further assumed that ¢,/¢; is independent of rough- 
6 ness, one would conclude that the ratio 2C,,/C; for rough tubes 
is approximately equal to unity and independent of Prandtl 
2 number. Unfortunately, the actual experimental evidence 
sccumulated by the present authors as well as by previous 
01 experimenters (see for example (2 and 3)) does not support 
r the above hypotheses, 

n For the types of roughness tested so far, it appears that 
friction is more strongly affected by roughness than is heat 
transfer and that there is no effective increase of the turbulent 
exchange coefficients within the laminar boundary layer. 
Nunner (2) suggests the explanation that each roughness 
t element, in addition to exerting surface drag, exerts a form 
drag. Vortices may be shed from these elements and flow 
disturbances will result. However, these disturbances are 
ipparently too far from the surface to effectively increase 
h the exchange coefficient €, near the surface. 

is As was stated in the introduction, it is hoped that the data 
presented in this paper will contribute to the large body of 


e experimental information which is still needed to clarify the 
e problem of heat transfer in a rough tube. Additional experi- 
f ments will be required and several aspects of the problem may 
e have to be studied. In view of the suggested role of form 
t drag, an investigation of the effect of the shape of the indi- 

vidual roughnesses seems particularly challenging. One 
is might even speculate on the possibility of developing rough- 
, ness shapes with low form drag, which should be especially 

suited to heat transfer applications. If successful, the hy- 

pothesis regarding the equality between C, and C,/2 for rough 
) tubes might still be realized in some cases. In spite of the 


fact that investigations involving a certain variety of rough- 
ness shapes have not so far shown any such favorable effects 
(2), further study of this aspect may still be justifiable. 
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Generalized Trajectories for Free- 
Falling Bodies of High Drag' 
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Nomenclature 


A = 
Fz = 
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area of body consistent with definition of Cp, ft? 
drag force, lb, 


g acceleration due to gravity, fps? 


go = gravitational constant = 32.2 lb,,ft/lb, sec? 
h = actual altitude, ft 
h* = reference altitude, ft 
h; = initial actual altitude, ft 
m = mass of body, lbm 
s = oblique distance, ft ; 
t = time, sec 
V = velocity, fps 
V; = initial velocity, fps 
x = vertical distance measured downward, ft 
= h* — hy ft 
a = constant, ft 
p = local air density, lb,/ft® =| 
po = standard sea level air density, lb,,/ft® 
6 = angle of descent between vertical and direction of flight, 
degrees. 

Dimensionless Numbers 
Cy = drag coefficient 
h = h/a 
= hi/a 
R= 3.85 x 10-2 = 

m cos 0 

V/V ag 
V; = V;/ V ag 
x = z/a 
&; = 


Introduction 


, ¥ PRESENT considerable interest exists in the utiliza- 
tion of high drag devices such as parachutes for the re- 
covery of missile equipment and as brakes for missiles them- 
selves. This has led to the increased need for trajectory in- 
formation on free-falling bodies over an extended range of 
drag-to-weight ratios. One specific application of this trajec- 
tory information involves the prediction of the temperature 
history of a high drag body under the influence of aerody- 
namic heating. 

The equation of motion for a free-falling body in a gravita- 
tional field may easily be established and solved assuming an 
isothermal atmosphere and constant drag coefficient. The 
resulting solution appears in series form. For bodies of low 
drag-to-weight ratio, such as projectiles, this series usually 
converges rapidly and easily yields the approximate trajec- 
tory; however, as the drag-to-weight ratio becomes greater, 
the series converges too slowly for convenient use of the series 
solution. However, a judicious selection of the reference sys- 
tem for the equation of motion will allow the graphical presen- 
tation on a single coordinate diagram of the trajectories of a 
body of high drag for a wide range of weight-to-drag ratios, 
initial altitudes and initial velocities. Such a diagram 
eliminates the necessity for repeated evaluation of the series 
for each new problem. This method and its application form 
the content of this paper. 


The Equation of Motion and Its Solution 


Vertically Free-Falling Body 

The equation of motion is formed by equating the time rate 
of change of momentum of the body to the sum of all forces 
acting, which for the case of a free-falling body are the gravi- 


tational force mg/go and the drag force Fa. ai bs 
d(mv) mg 
godt Jo 
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In the conventional choice of coordinates for this equation, 
distance is measured upward from a coordinate system fixed 
at the surface of the earth. The method presented in this 
paper deviates from common practice and fixes the coordinate 
system at some point in space h* and measures the distance 
traveled downward, z, from this reference altitude (see Fig. 1a). 
The true altitude h of the body is then 


This choice of coordinate system makes possible the general- 
ized graphical presentation of trajectory data as will be dis- 
cussed. 

The drag force Fg may be expressed in terms of the dimen- 
sionless oe coefficient Cp by 


dinsiiaies Equation [3] into Equation [1] and assuming 
the mass m to be constant, there results 


CpApV? 


Changing the independent variable of time ¢ to distance 
measured downward z by use of the relation 


qv _dVdz_ av 
ze 


_yd@V _ 1 
‘dt dz (2 


dt dt dx 
[4] now becomes 


dV? cater 
a = 
x 


This differential equation cannot as yet be solved since the 
density p is a function of the distance x. However, if an iso- 
thermal atmosphere is assumed, the density then varies with 
altitude A in the following way 


= pe/* = ** (6] 
Introducing this expression into Equation [5] gives 
dx m 


This differential equation may now be expressed in dimen- 
sionless form if the following basic dimensionless numbers are 
defined and the quantities g and @ are assumed to be con- 
stants 


Equation [7] now assumes the dimensionless form i 


+ eS w= 2............ [8] 
dk m 


This equation is to be solved subject to the boundary con- 
dition 

= V;when# =0.................[9] 

Inspection of Equation [8] and its boundary conditions [9] 


reference oltitude _feference oltitude 


(a) (b) 


Fig. 1 Reference systems for equation of motion. (a) Verti- 
cally free-falling body. (b) Obliquely free-falling body 
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reveals that the solution will be a function of the dimension. 


less groups 


If the coefficient of drag Cp is assumed ed constant, the | the di- 
mensionless parameter 


CpApoa 
m 


becomes independent of the distance z, and Equation [8] may 
be solved by standard methods. Its solution appears in th« 
following infinite series form 


V2 — = 2e-Ke*® X 


This series solution may be easily evaluated for bodies of 
low drag and large mass, such as projectiles, since the series 
then converges rapidly. However, for bodies of high drag, 
such as parachutes, the infinite series converges too slowly to 
allow for convenient use of this equation. In such instances 
a graphical presentation of the general solution would be use- 
ful, since having once been prepared it could be repeatedly 
used. In attempting to accomplish this, one difficulty imme- 
diately arises from the fact that Equation [10] contains four 
variables. If this number could be reduced to three, a simple 
graphical presentation of the solution would become possible. 
Consideration of the four dimensionless numbers appearing in 
Equation [10] reveals that no limitation is imposed on the 
generalized nature of the solution if the dimensionless group 


CpApoa 
m 


is held constant. This may appear at first to fix the drag-to- 
weight ratio CpA/m but such is not the case. By virtue of the 
choice of reference system (see Fig. 1a) the reference altitude 
h* is allowed to vary and may take on any value greater than 
or at least equal to the true initial altitude of the body; there- 
fore, holding the above dimensionless numbers constant 
simply fixes the relation that must exist between the drag-to- 
weight ratio and the reference altitude h* as given by 


CoA poe | nt/a _ 


where the terms pp and @ are constants. The number K may 
be somewhat arbitrarily chosen. In the present analysis a 
value of K = 3.85 X 10~? was selected in consideration of 
certain scaling needs. With the afore-mentioned dimension- 
less group held fixed, the solution to the equation of motion 
now contains only three variables 


This makes possible a graphical presentation of the solutions 
as a family of curves on a single coordinate diagram as shown 
in Fig. 2. These curves are actually the trajectories of a free- 
falling body for a wide range of initial altitudes, initial veloci- 
ties and weight-to-drag ratios. 


Obliquely Free-Falling Body 


The equation of motion and solution developed for the 
special case of a vertically falling body may easily be extended 
to the more general case of an obliquely falling body (see Fig. 
lb). For such a condition the equation of motion may be 
written as 


d(mV) mg cos 6 
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Fig. 2 Generalized trajectory diagram 


In terms of the drag coefficient, and assuming constant 
mass, Equation [12] becomes 


Changing the independent variable of time ¢ to distance 
measured vertically downward (x) by use of the relations 


dV dVds_ dV dx ds V cos 6 dV? 
di ds dc dsdt dt 2 de 
Equation [12] now becomes 
CpApV? 
(14] 
m cos 6 


By assuming the same variation of air density with altitude as 


previously given by Equation [6] and by using the same di- 
mensionless numbers as before, Equation [14] becomes 


dV? CoApre 


e = 2 
d& m cos 6 


This equation and its solution are identical to those pre- 
viously given for the vertically free-falling body with the ex- 
ception of the presence of the cos @ term in the coefficient. 
The graphical solutions given in Fig. 2 are still applicable. 
However, the relation between the drag-to-weight ratio Cp 
and the reference altitude h* is now given by 


6 
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This represents the more general situation of which the 
vertically falling body is a special case for cos @ = 1. Graphi- 
cal presentation of this equation appears in Fig. 3, which pre- 
sents CpAa/m as a function of h*/a for various angles of 
descent. 

The next section discusses in more detail the characteristics 
of the generalized trajectory plot given in Fig. 2 and its ap- 
plication in solving for trajectories. 


Generalized Trajectory Diagram for Bodies of High Drag _ 


The graphical presentation of the solution to the trajectory 
equation is given in Fig. 2 and applies over an pees 
range in drag-to-weight ratio CpA/m from 4 to 4 X 10° 
ft?/lb,,, in initial altitude from sea level to 250,000 ft, and in | 
initial velocity from zero to 10,000 fps. This figure presents — 
the dimensionless velocity 7 vs. the dimensionless distance = _ 
for various dimensionless initial velocities with the dimension- | 
less number K held constant at a value of 3.85 X 10-2. The | 
resulting curves represent the solution to the equation of | 
motion given in Equation [10]. 

To simplify use of these curves, the basic dimensionless 
numbers on the coordinates have been evaluated for the usual 
‘ase of g = 32.2 fps? and a = 25,000 ft with this scale also — 
presented in addition to the dimensionless coordinates. It 
should be noted that the resulting dimensionless scales are 
slightly restricted by the choice of a = 25,000 ft. However, 
if the altitude range of interest is more limited, the best a 
ble approximation may be a slightly different value of a. 

In Fig. 2 the ordinate scale begins at the top with ¢ = 0 and | 
all distances are measured downward from this point. In 
order to apply significance to such a measurement it is neces- 
sary to fix the origin = 0 with respect to the surface of the 
‘arth. This is accomplished through use of the reference al-— 
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titude h* which appears in the dimensionless number K 


CpA po 


m cos 


For. convenience in calculation, Equation [16] has been in- 
terpreted graphically in Fig. 3 and ,* may be read directly 
from this figure once the drag-to-weight ratio and angle of 
descent are specified. The reference altitude in this way 
established the origin of the ordinate scale with respect to the 
surface of the earth. In general, the reference altitude will 
not be equal to the initial altitude of the body, the only re- 
quirement being that h* must not be less than the initial al- 
titude. Once this ordinate scale has been fixed in the above 
manner for a body of specified drag-to-weight ratio, the prob- 
lem of applying these curves to situations involving various 
initial altitudes and velocities has to be considered. In this 
regard it should be noted that a characteristic of the trajec- 
tory of a free-falling body makes it possible to enter the curve 
representing the trajectory once it has started, provided one 
enters with a body of the same drag-to-weight ratio at an 
initial velocity corresponding to the local velocity for that al- 
titude. Consideration of the equation of motion shows this 
to be true since the forces acting on the body in both cases will 
be the same. With this in mind it is now possible to enter Fig. 
2 at any initial altitude less than the reference altitude and at 
any initial velocity up to 10,000 fps, provided the ordinate 
scale has been fixed by the weight-to-drag ratio as previously 
mentioned. 

In the following section several sample problems will be 
worked to illustrate the use of the trajectory diagram. 


1 Vertically Free-Falling Body of High Drag ae 


Sample Problems 


Given: 


Drag-to-weight ratio CpA/m = 1 ft?/Ib,, 
Initial velocity V; = 5000 fps 

Initial altitude h; = 200,000 ft 

Angle of descent 6 = 0 deg 


Find: Velocity vs. altitude of the body. 

Step 1: Assuming a value of a = 25,000 ft enter Fig. 3 at 
CpAa/m = 25,000 for 6 = 0 and read off h* = 10.84. The 
point = 0 on the trajectory diagram now corresponds to an 
actual altitude of 


Step 2: Establish the relationship between the true altitude 
h and the distance traveled downward zx through use of the 
general equation 
h=h*—-x or h=h* 

which for this special problem becomes 
271,000 

a 


or x = 271,000 —h 


= 


Step 3: Enter the trajectory diagram on Fig. 2 by locating 
the points 


5000 
25,000 X 32.2 
4 Ordinate: X; = 10.84 onde 
ag rdinate: Aj = 
25,000 


or ad ‘etal V; = 5000 fps 
Ordinate: X; = 271,000 — 200,000 = 71,000 ft 


ll 


Step 4: Select the trajectory curve passing through the 
point as located in Step 3. In most instances it will be neces- 
sary to interpolate between existing curves. Read off the 
value for the dimensionless velocity 7 for selected values of < 
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which are then made dimensional through multiplying by 
ag and a, respectively. The distance x is converted to 
true altitude h, as before, through the equation 


h=h*—-2 h = 271,000 — x 


Results: See Table 1. a 
Table 1 Velocity vs. altitude data for parachute in sample 
problem 1 

ft z,! ft 3 V,! fps 
200,000 71,000 2.84 5.56 5000 
175,000 96 , 000 3.84 3.36 3000 
150,000 121,000 4.84 1.01 900 
125,000 146,000 5.84 0.38 345 
100,000 171,000 6.84 0.25 220 
75,000 196,000 7.84 0.14 125 
50,000 221,000 8.84 
25,000 246, 000 9.84 
Sea level 271,000 10.84 
1 For a=25,000 ft and g = 32.2 fps?. 


2 Obliquely Free-Falling Body of High Drag 
Given: 


Drag-to-weight ratio CpA/m = 1 ft?/lb,, 
Initial altitude h; = 200,000 ft 
Initial velocity V; = 5000 fps 
Angle of descent 8 = 30 deg - 


Find: Velocity vs. altitude of the body. 

Step 1: Assuming a value of a = 25,000 ft, enter Fig. 3 at 
CypAa/m = 25,000 for @ = 30 deg, and read off p* = 11.00. 
The point ¢ = 0 on the trajectory diagram now corresponds 
to an actual altitude of 


h = h* = h*qa = 275,000 ft 


Step 2: As in the case of a vertically falling body establish 
the relationship between the true altitude h and the distance 
traveled downward z through use of the general equaton 


or h=h* 


which for this problem becomes 


h = h* = ah* = 271,000 ft ~~ - ‘+ h= “lie or ¥ = 275,000 — x 


Step 3: Enter the trajectory diagram on Fig. 2 by locating 
the point 


x = 2x; = 275,000 — 200,000 = 75,000 ft 


Step 4: Select the trajectory curve passing through the 
point as located in Step 3. In most instances it will be neces- 
sary to interpolate between existing curves. Read off the 
value for the dimensionless velocity 7 for selected values of z 
which are then made dimensional through multiplying by 
Vag and a, respectively. The distance x is converted to true 
altitude h as before through the equation 

h=h*—-—2z or 
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Kodak reports on: 


running up the flag of arsenic-cross-linked selenium. ..a use for Biphenyl in 


aerodynamics research... big deal over 0.009 roentgen 


Out to 


We find that mechanically and op- 
tically decent inorganic polymers 
can be made out of sulfur’s lesser 
known sister, who was named after 
the moon. Acquaintance with ar- 
senic-cross-linked selenium has rip- 
ened into intimacy. The issue is a 
kind of glass with refractive index 
2.45 or higher. Black as pitch by 
reflected visible light, it has a trans- 
mittance that is optically useful out 
to a prodigious 25y, practically 
long enough to modulate with your 
favorite TV situation comedy. 

This 92% Se - 8% As glass is 
quite a different animal from arse- 
nic trisulfide, which begins falling 
off at 9u. Except for an absorption 
band at about 12.6u, most of its 
transmission loss arises from Fres- 
nel reflection due to the high refrac- 
tive index, and there are said to be 
ways of fixing that. The high index, 
on the other hand, is a blessing 
when infrared detectors are em- 
bedded in the glass in order to 
boost their geometrical radiation- 
collecting effectiveness. 

Kodak Selenium Glass retains its 
rigidity to about 70 C and does not 
crystallize after long exposure to 
this temperature. 

Though the production scale is 
still modest, we would like to try 
supplying it in molded shapes to 
manufacturers determined to outdo 
us on optical systems employing it. 

Potential competitors and other 
friends are most cordially invited to 
submit inquiries about Kodak Selenium 
Glass to Eastman Kodak Company, 
Apparatus and Optical Division, Roch- 
ester 4, N. Y. 


Wind writing 


“Photography” means “light writ- 
ing.” To get light to write meaning- 
ful records of air movement, the 
wind tunnels of the world are fes- 
tooned with complex and delicate 
optical arrangements. 

Why not use the wind itself to 
write with? the men of the Ames 
Aeronautical Laboratory of the Na- 
tional Advisory Committee for Aer- 
onautics asked themselves. Like the 
wind patterns in snow around a 
fence. It worked fine. 


At (a) is a wind-tunnel model before 
and at (b) after spraying with a thin 
coating of Biphenyl; at (c), after 
the wind has etched into the sub- 
limating stuff a pattern that tells of 
laminar, turbulent, and separated 
flows and locates the shock waves. 

Our connection with this is that we 
supplied the Biphenyl. Not that Bi- 
pheny! is big business for us. (Color and 
black-and-white film for schlieren and 
interferometric studies in aerodynamics 
is much bigger business!) The point we 
make is that somebody had an idea re- 
quiring a certain ofganic compound, 
looked into our catalog of some 3600 of 
them, and found it there. You, too, can 
get a copy by asking Eastman Organic 
Chemicals Department, Distillation 
Products Industries, Rochester 3, N. Y. 
(Division of Eastman Kodak Company). 


A modest contribution 


Man comes home from work and 
his wife asks him what kind of a 
day he’s had. “Same old yackety 
yak,” he reports. ‘“Heap-big all- 
afternoon powwow over the fact 


that whereas it used to take 0.019 
roentgens of radium gamma radia- 
tion to produce a density of .05, 
they have now made x-ray film that 
takes 0.010 roentgens of gamma ra- 
diation to produce a density of .05. 
Cauliflower tonight again?” 
e 

Eastman Kodak Company is 
pleased to announce that significant 
progress has been made in its con- 
tinuing program of improving the 
sensitivity of film to ionizing radia- 
tion. One result of this program is 
that physicians are enabled to ob- 
tain required diagnostic informa- 
tion from Kodak No-Screen Medi- 
cal X-ray Film with fully 50% less 
exposure of the patient to x-rays. 

Another result, of interest to 
those engaged in industrial radiog- 
raphy and in the atomic energy 
field, is that the film badges worn 
to check on their body exposure to 
radiation can now be made capable 
of at least 50% more sensitivity 
than previously. 

While it is realized that the end 
of the effort to reduce human expo- 
sure to radiation is nowhere in sight, 
it is felt that even modest techno- 
logical advances in this direction 
serve the interests of general wel- 
fare. 

For cumulative monitoring of 
occupational exposure to x-rays 
and gamma rays, the enhanced 
sensitivity is now available in Kodak 
Personal Monitoring Film, Type 1. 
Maximum sensitivity is attained for 
radiation quality of approximately 
50 kvcp. Approximately 30 times 
the same roentgen value is required 
to produce equivalent densities with 
radium gamma rays. 

In Kodak Personal Monitoring 
Film, Type 2, with a low-sensitivity 
emulsion on one side of the film 
and a high-sensitivity emulsion on 
the other side, the range of assess- 
ment is extended to emergency 
levels of exposure. Other modifica- 
tions are available for the measure- 
ment of neutron exposure. 


Full information about these monitor- 
ing films is obtainable from Eastman 
Kodak Company, Special Sensitized 
Products Division, Rochester 4, N. Y. 


This is another advertisement where Eastman Kodak Company 


probes at random for mutual interests and occasionally a little 
revenue from those whose work has something to do with science 
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New Patents 


Inwardly opening breech closure for rocket 
launcher (2,820,398). W.E. Good, River- 
ton, N. J., assignor to the U.S. Army. 

Launcher tube with a breech housing 
secured to a flanged rearend. The breech 
block has a conical surface to engage the 
inner part of a conical seat with a tight 
fit, providing a gas-tight seal. 

propellent support (2,820,410). 
Tarr, Sierra Madre, Calif., assignor 
to the War Dept. 

A member is Fostted around a reinforc- 
ing base and sealed against the walls of the 
tube. A propellant grain fits within the 
tube and is bonded to the sealing member. 
Axial flow compressors (2,820,588). A. 
J. Penn and L. A. Nevard, Liverpool, 
— assignors to D. Napier & Son, 

uti 

Disengageable coupling means between 

first and second rotors. When the second 
stage rotor attains a predetermined speed, 
windmilling of the first compression part 
is stopped and becomes driven with the 
later stages. 
Dual unit jet propulsion plant (2,820,599). 
J. Ackeret and C. Keller, Kusnachi, 
Switzerland, assignors to Aktiengesell- 
schaft fuer Technische Studien. 

Combination of a turbojet power unit 

in which the exhaust is directed as a pro- 
pulsion jet. A further jet propulsion unit 
ae ises a turbine and compressor where- 
by the — air is directed as a pro- 
pulsion jet. 
Automatic jettisoning system for aircraft 
seats (2,820,602). KE. V. Foster, Man- 
hattan Beach, Calif., assignor to Northrop 
Aircraft, Inc. 

Personnel survival system consisting of 
ejectable seats catapulted upon release of a 
pressurized fluid initiated by a gas genera- 
tor. 


A 


Combustion chamber construction in a gas 
turbine engine (2,821,067). H. C. Hill, 
Seattle, Wash., assignor to Boeing Air- 
plane Co. 

Each combustion chamber comprises an 
outer wall opening at one end to the col- 
lector chamber. An inner shell within the 
wall opens at one end to the nozzle box. 
A burner is admitted to each inner shell 
at the end away from the box. 
power plant (2,819,765). E. 

Chatterton, London, England, assignor 
to D. Napier & Son, Ltd. 

Master control connected with an ad- 
justable governor to regulate the speed of 
a reciprocating engine. The control also 


regulates the transmission ratio between 
the engine crankshaft and the compressor 
and turbine rotors. 

Exhaust temperature regulator for gas 
turbine power plant (2,821,065). N. E. 
Starkey, Schenectady, N. Y., assignor to 
General Electric Co. 

Means of controlling the rate of fuel con- 
sumption in the hot gas generating area 
as a function of load output desired from 
the turbine rotor. 


2,821,351 : 


Airplane wing structure embodying jet 
engine (2,821,351). V. V. Utgoff, Alex- 
andria, Va. 

Pulsejet units side-by-side inside the 
wing. The top surface skin provides an 
intake. An exhaust communicates from 
the jets to the wing exterior spanwise of 
the trailing edge. The major portion of 
exhaust flow is deflected below a trailing 
edge flap. 

Firing system for jet type perforating gun 
(2,821,136). J. H. Castel, Houston, Tex., 
assignor to P. G. A. C. Development Co. 

Shaped charges and a cable for lowering 
a perforating gun into a borehole. Explo- 
sive charges are sequentially fired by trans- 
mitting direct and alternating currents to 
the gun. 


Military rocket structure (2,821,137). 
G. A. Lyon, Detroit, Mich. 
Bowed spring strips with helically 


angled edges movable laterally relative to 
the projectile body. Edges of the strips 
cooperate with rifling in the projectile 
barrel. 

Dive control system for pilotless aircraft 
(2,821,349). H. Sohn, Essex, Md., as- 
signor to The Martin Co. 

Torque producing means for precessing 
a gyroscope at a predetermined variable 
rate. Torque is modified by accelerations 
due to wind lift. Aerodynamic means 
responsive to the gyroscope position con- 
trol the attack angle of the aircraft. 

Jet airplane construction (2,821,350). J. 
Smurik, Fort Lee, N. J. 

Fins on each side of the body longitudi- 
nally dividing the space into a central 
channel with a channel on either side. 
Combustion products from the central 
burner are conducted to the side burners, 
raising their temperature to facilitate 
starting. 

Fin stabilized projectile (2,821,924). L. J. 
Hansen and P. Rosenberg, Chicago, IIl., 
assignors to the U. S. Army. 

Shell comprising a boat-tail fixed to the 
rear and carrying nested extensible fins 
imbedded in the propelling charge. Fins 
extend outward upon movement of pistons 


+t George F. McLaughlin, Contributor 


actuated in response to differential pres- 


sure between the chamber and the atmos- 
phere. 

Improvements in a device for starting gas 
turbine plants (2,822,665). C. R. Nicolin, 
Finspond, Sweden. 

Separate compressor in front of the air 
intake of a turbine for delivering com- 
pressed air. A slidably mounted sur- 
rounding cover is placed into engagement 
to form a closed air passage between tlie 
compressors during the starting period. 
Propellant feed system (2,822,667). 
Drexel, Los Angeles, Calif., assignor to tlie 
Garrett Corp. 

Container with a movable wall dividing 
it into propellant and gas chambers. \ 
combustion chamber in the gas chamber 
permits direct transfer of heat energy to 
cause thermal expansion of the gas. 
Aircraft with jet propulsion engine and gun 
firing means (2,822,732). T. L. Metcalfe, 
Ashbourne, England, assignor to Rolls- 
Royce Ltd. 

Gun firing means which opens a valve 
to bleed off from the delivery of the com- 
pressor air compressed when the gun is 
actuated. 

Thermal arrangement for rocket missiles 
(2,822,756). N. G. Kuller and K. T 
Thorildsson, Bofors, Sweden, assignors to 
Aktiebolaget Bofors Corp. 

Electric terminals defining a groove fo: 
connecting the detonating device with « 
distant power supply. Terminals have 
outwardly tapering bushings with portions 
protruding into the groove. 


Antenna mounting for a guided missile 


(2,821,706). H.J. Rels and W. V. Foley, 
Middle River, Md., assignor to The 
Martin Co. 


Stub pattern antennas for transmitting 

radiant energy. A disk-shaped fairing 
forms a ground plane whereby radiation 
extends above a horizontal plane passing 
through the antennas. 
Cursor generator with a height finding 
adaptation (2,823,374). L. S. Michels, 
Inglewood, Calif., assignor to Gilfillon 
Bros., Inc. 

Production of pulses representing a 
cursor line for use as a height finder. Each 
pulse corresponds to the intersection in 
time between a first signal representing the 
cursor angle and a second variable angle 
signal representing a scan angle position. 
Target seeker head for guided missiles 
(2,823,612). A. Cox and C. Ledda, New 

‘ork, N. Y., assignors to the U. S. Air 
Force. 

Pyramidal-shaped body with flat glass 
plate side windows separated by narrow 
stringers at the base. Light rays focus on 
a plane through the central opening of an 
annular concave mirror for the use of the 
control mechanism. 


Eprtrors Note: 


Patents listed above were selected from the Official Gazette of the U.S. Patent Office. 


Printed copies of patents may 


be obtained from the Commissioner of Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 


268 


JET PROPULSION 


: 
a 
. 
‘ 
Dk 
= 
i 
/ 


y 


© Armament 
© Ballistics 
2 Radar Antennas 


© Guided Missile 

Support Equipment 

@ Auxiliary Power Supplies 
7 © Control Systems 


Aprit 1958 


Throughout the Pattern of Modern Defense... —_ 


AMF has experience you can use 


e Most defense problems fit into AMF’s “big picture”...a picture drawn from AMF’s 
tremendous backlog of experience in defense work. Wherever you are, there’s probably 7 

an AMF defense component or integrated system...a product of AMF experience... | 

operating near you. e Guided missiles ride on AMF equipment to AMF launching sites 

...to be loaded, fueled and fired by AMF-built equipment. Under the sea, AMF-built = 
weapons wait for unfriendly submarines on the prowl. Along our borders, AMF radarscopes 
search the sky for “stranger” aircraft. @ There is little room for failure where the job is the | 
nation’s defense. And the nation looks to companies like AMF to design, test and produce 
a variety of defense products. With a wide range of experience in the most sensitive fields 
of defense work, AMF may well be the answer to your problems...in defense or industry. 
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Book Reviews 


Principles of the Properties of Materials, 
by J. P. Frankel, McGraw-Hill, New 
York, 1957, 228 + xiv pp. $6. 

Reviewed by Dona.p 8. CLarK 
California Institute of Technology 


This is an interesting book based upon 
a philosophy to which all educators will 
not subscribe. Some teachers of en- 
gineering materials will agree with the 
author that one should concentrate on 
principles rather than practice or tech- 
niques. This is the more modern way of 
looking at engineering training. Such an 
approach becomes necessary because of the 
more scientific character and stringent re- 
quirements of engineering. While much 
of the practice of engineering is empirical, 
more of it can be approached from an 
analytical point of view. The student 
should have a good grounding in principles 
which will permit him to understand pres- 
ent practices and by which he may be 
able to obtain better solutions to complex 
problems. This is the basis on which this 
book has been written. In some two 
hundred pages, Professor Frankel has 
attempted to present the principles of 
knowledge concerning the properties of 
materials. 

The entire approach to the subject is one 
of apparent sophistication. This is an 
attitude that should not be laid on so 
thick. The author has probably over- 
emphasized the basic concepts without 
sufficient illustration to clarify the sig- 
nificance of the approach. While the 
material presented is interesting, it has its 
deadly aspect as far as the student is con- 
cenred. Much of the material could be 
expanded, particularly as to application, 
and yet an instructor would be necessary 
to assist in interpretation. 

The author states in the preface that 
the book is a text, not a treatise. Cer- 
tainly he is correct that it is not a treatise. 
However, one may also question if it is a 
text. A text is a complete story for the 
purpose of educating the reader in a partic- 
ular subject. It may or may not be a 
treatise. This reviewer is of the opinion 
that the book is too much of a survey text 
to be of real value for the education of 
sophomores in engineering. 

The book lacks clarity in many in- 
stances, particularly in definition. For ex- 
ample, reference is made to phases, but 
there is not an adequate definition of a 
phase. Some reference is made to age 
hardening, but the student is not made 
familiar with the mechanics of age hard- 
ening in preparation for the discussion of 
such a property. There seems to be an 
atmosphere of fear of mentioning specifics 
in spite of the existence of some excellent 
examples of practical applications. The 
sophisticated approach seems to be of 
more importance than a clear understand- 
ing of application of principles of physics 
to materials. 

The subject is covered in two parts, the 
first on physicochemical foundations and 
the second on properties. The first part 
is a good discussion of fundamental con- 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


cepts starting with certain phases of 
atomic physics, thermodynamic considera- 
tions, crystal structure, cohesion and the 
metallic state. These subjects are covered 
very nicely as a survey, but possibly with 
a “holier than thou’’ attitude, and with- 
out sufficient depth. 

The second part of the book on prop- 
erties is fine, but it does not show the stu- 
dent the real significance of these prop- 
erties in engineering. For example, an 
understanding of phase relations is not 
easy to grasp by the beginner. The prin- 
ciples laid by the author are not sufficient 
to permit the student to secure a basic 
understanding of phase diagrams, etc., 
by merely reading another book. The 
subjects covered include conductivity, 
heating, diffusion, electrochemistry, elas- 
ticity, plasticity, fracture, viscosity and 
general subjects. 

The author is to be complimented on 
bringing this type of material to the sopho- 
more level. The criticism given in this 
review is not intended to discourage the 
use of the author’s approach, but to point 
out certain basic deficiencies. Possibly 
in the next edition, the material can be ex- 
panded so as to clarify many of the points. 
This review recommends that every teacher 
of the subject of engineering materials 
not only read this book, but make use of 
much of the basic material in his course. 
Such use will mean decreasing the amount 
of material given on specific reduction 
practices and composition of specific alloys 
and their properties. A reduction of this 
type of material is possible without 
seriously interfering with a proper educa- 
tion of our future engineers. The follow- 
ing statement at the close of the book can 
well be a closure for this review, because 
many teachers are not sufficiently cogni- 
zant of its importance: ‘‘Not only will the 
student be called on soon to solve prob- 
lems his teachers could not solve, but also 
he must be able to solve problems his 
teachers could not even imagine.”’ 


Exploring the Atmosphere’s First Mile. 
Vol. I: Instrumentation and Data 
Evaluation; Vol. II: Site Description 
and Data Tabulation, by Heinz H. Let- 
tau and Ben Davidson, Pergamon Press, 
New York, 1957, 376 pp. (Vol. I), 578 
pp. (Vol. II). 

Reviewed by 8. F. Sincer 
University of Maryland 


The volumes report the proceedings of 
the Great Plains Turbulence Field Pro- 
gram carried out in Nebraska. These were 
elaborate meteorological experiments per- 
formed in open prairie country with the 
primary objective to study the profiles of 
wind, temperature and other meteorologi- 
cal quantitities in the lowest part of the 
atmosphere. The subject matter properly 
falls into the fields of micrometeorology. 
Volume I contains an introductory article, 
giving history, summary and general con- 
clusions, and 65 scientific papers which 
describe in detail the multitude of instru- 
ment types used in the field program. 


Volume II presents the data of the seven 
observational periods. 


The Science of Engineering, edited by J. 
E. Goldman, John Wiley & Sons, New 
York, 1957, 528 pp., $12. 

Reviewed by J. P. FRANKEL 
University of California 


This book contains a series of lectures 
demonstrating the application of know!- 
edge in chemistry and physics to the 
understanding of the properties of mate- 
rials. The first four lectures deal with the 
applications of solid state science in engi- 
neering and basic information as to the 
structure of the material. Lectures 5 to 10 
are concerned with metals and alloys 
lecture 11 with surfaces; the remaining 
seven lectures are divided almost entirel, 
among magnetic properties, properties of 
semiconductors and dialectrics, and non- 
crystalline materials. 

The treatment is that of a review o! 
basic results in the field. Few or no der- 
ivations are enclosed. Generally, the 
results are stated in the “it can be shown’ 
style. For complete understanding of the 
subject material, the reader must bring 
to the book an understanding of the basic 
concepts of quantum mechanics and sta- 
tistical thermodynamics. The treatment 
in lecture 2 of the foundations of quantum 
mechanics is inadequate for the reader with 
no background and superficial to those who 
have already studied quantum mechanics. 
There is no review of statistical thermo- 
dynamics corresponding to the review of 
quantum mechanics. 

Not withstanding these defects (which 
may perhaps be a result of the different 
authorship of different chapters), the book 
is an excellent summary of the methods 
used to investigate and rationalize the be- 
havior of the materials. With respect to 
the mechanical properties of solids, the 
treatment is more complete than that 
usually found in text books on solid state 
physics. 

In summery, the book is an excellent 
reference for people in other fields who 
would like a quick view of the basic pro- 
cedures of solid state physics, and there- 
fore can be considered a “must’’ for teach- 
ers of courses in the properties of materials 
in engineering colleges. 


Solar Eclipses and the Ionosphere, by W. 
J. G. Beynon and G. M. Brown, Perga- 
mon Press, New York, 330 pp. $21. 

Reviewed by S. F. SINGER 
University of Maryland 


This volume forms No. 6 of a series of 
supplemental volumes to the Journal of 
Atmospheric and Terrestrial Physics. The 
first one was “Rocket Exploration of the 
Upper Atmosphere’’; another in the series 
has been “Vistas in Astronomy,” a two- 
volume edition. 

The present work is the result of a sym- 
posium held in London in August 1955 to 
discuss in the fullest possible way the 
effects of solar eclipses on the ionosphere. 
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1,500 mile 


ofa 


thrust. 


The very heavens shake with the thunder of 
some 135,000 horses as Thor streaks across 
the sky at 10 times the speed of sound. De- 
signed by Douglas Aircraft to deliver total 
destruction to targets as far away as 1,500 
Thor represents the Air Force’s 


miles, 


Tmk(s) 


Preserver of Peace ... 


striking arm where ground objectives are 
concerned. For this surface-to-surface inter- 
mediate range ballistic missile, with its 
inertial guidance system, RCA has developed 
and is supplying electronic units to help Thor 
declare its mission: the prevention of war. 
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physicals, or smoothness 


a beyond the ordinary 

you should be interested 

in the unusual 
foundry methods of 

Morris Bean & Company 

Yellow Springs 6, Ohio se 


Such observations are of great value in the 
interpretation of many geophysical phe- 
nomena. This book gives a full account 
of the various influences of solar eclipses 
on the earth’s atmosphere. 

The effects, of course, relate mainly to 
the sudden removal of solar ultraviolet and 
x-rays which normally cause the iono- 
sphere. There are therefore some theo- 
retical considerations on these points fol- 
lowed by a large section giving results of 
ionosopheric radio observations during 
the eclipses. Then follows a series of 
theoretical papers which discuss the radi- 
ation from the sun and its production of 
the ionosphere. Next, there are four 
papers dealing with the recombination 
process in the ionosphere which removes 
ions and electrons. 

Because of the fact that the ionosphere 
carries electric currents and contributes to 
the earth’s magnetic field, eclipses also give 
geomagnetic effects. These are discussed 
in a series of four papers. An interesting 
set of papers deals with solar radio noise 
during the eclipse. 

An extensive bibliography which goes 
back as far as 1912 lists the literature of 
solar eclipses and the ionosphere. 


Books Received 


Dynamics of Machinery, by A. R. 
Holowenko, Wiley, New York, 1955, 
461 pp. $7.50. 

Glass Reinforced Plastics, by Phillip 
Morgan, Philosophical Library, New York, 
1957, 276 pp. $15. 

Basic Mathematics for Radio and Elec- 
tronics, by F. M. Colebrook and J. S. 
Head, Philosophical Library, New York, 
1957, 359 pp. $6. 


Radio Aids To Air Navigation, by J. H. 
H. Grover, Philosophical Library, New 
York, 1957, 138 pp. $6. 

The Expanding Case for the UFO, by 
M. K. Jessup, Citadel Press, New York, 
1957, 253 pp. $3.50. 

Arizona’s Meteorite Crater, by H. H. 
Nininger, World Press, Denver, 1956, 
232 pp. $3.75. 


Industrial Electronics Circuits, by R. 
Kretzmann, Philosophical Library, New 
York, 1957, 194 pp. $10. 

Case Problems in Transportation Man- 
agement, by George P. Baker and Gayton 
IE. Germane, McGraw-Hill, New York, 
1957, 685 pp. $8.50. 


Radio Telemetry, by M. H. Nichols and 
L. L. Rauch, 1956, 461 pp. $12. 


Control-system Dynamics, by Walter 
R. Evans, McGraw-Hill, New York, 1954, 
282 pp. 


The Chemistry of Organometallic Com- 
pounds, by Rochow, Hurd and Lewis, 
John Wiley & Sons, New York, 1957, 344 
pp. $8.50. 

Advances in Geophysics, Vol. II, by 
H. E. Landsberg, Academic Press, New 
York, 1957, 378 pp. $8.80. 

Digital Differential Analyzers, by 
George F. Forbes, published by G. F. 
Forbes, Pacoima, Calif., 1957, 154 pp. 


Modern Chemistry for the Engineer and 
Scientist, by G. R. Robertson, McGraw 

Hill, New York, 1957, 442 pp. $9.50 
This volume is a compilation of the lec- 
tures given at U.C.L.A. by the following 
outstanding scientists: Kenneth 

Pitzer, Paul H. Emmett, W. Albert Noyes 

Jr., Norman Davidson, Clifford S. Garner 

Laszlo Zechmeister, Saul Winstein, Her 

man F. Mark, Nelson W. Taylor, Eugen 

G. Rochow, Theodore W. Evans, Georg: 

C. Kennedy, John E. Dorn, Joseph A 

Pask, Donald K. Tressler, Robert L. Met 

calf, Gordon A. Alles, Theodore A. Geis- 
man, Linus Pauling. 


Elements of Heat Transfer, by Max 
Jakob and George A. Hawkins, John Wile, 
& Sons, New York, 1957, 317 pp. $6.75 
This constitutes the third edition of the 
well-known elementary textbook. It in- 
cludes new data, and places more emphasis 
on analysis of heat transfer problems. 


Advances in Geophysics, by H. E. 
Landsberg, Academic Press, New York 
1956, 378 pp. $8.80. This book shoul 
be of interest to all geophysicists, particu- 
larly those involved in the IGY. It con- 
tains the following contributions: 

A. P. Crary—Arctic Ice Island Research 

Zdenek Sekera—Recent Developments 

in the Study of the Polarization of 
Sky Light 

Perry Byerly—Subcontinental Structure 

in the Light of Seismological Evidence 
E. C. Bullard, A. E. Maxwell and R. 
Revelle—Heat Flow through the 
Deep Sea Floor 
J. A. Jacobs—The 
Earth 

P. H. Jones, H. E. Skibitzke—Subsur- 
face Geophysical Methods in Ground- 
Water Hydrology 

S. F. Singer—Geophysical Research 

with Artificial Earth Satellites 


Interior of the 


Internationale Tagung Uber Stau- 
strahlen und Raketen, edited by Eugen 
Sanger, Irene Sanger-Bredt, Ernst von 
Olnhausen, Germany, 1956, 358 pp. This 
volume constitutes the proceedings of the 
International Conference on Jets and 
Rockets held in Freudenstadt, Feb. 6-8, 
1956. It includes the addresses and 
technical contributions of the following: 
Franz Ludwig Neher, B. Eckert, H. C. 
Seebohm, A. Seifriz, H. Veit, E. Sanger, 
Maurice Roy, J. Dupin, I. Sanger-Bredt, 
R. W. Porter, J. Corbeau, W. Michely, 
KE. Sharer, A. Gerber, W. Pilz, L. J. 
Sedov, G. G. Czerny, H. Maecker, Th. 
Peters, H. J. Kaeppeler, F. Winterberg. 

This compilation is a contribution of the 
Forschungsinstitut fiir Physik der Strahl- 
antriebe. 


Higher Oxo Alcohols, by L. F. Hatch, 
John Wiley & Sons, New York, 1957, 120 
pp. $4.50. 


Principles of Electrical Measurements, 
by H. Buckingham and E. M. Price, 
Philosophical Library, New York, 1957, 
600 pp. $15. 

Automation: Its Purpose & Future, 


by Magnus Pyke, Philosophical Library, 
New York, 1957, 191 pp. $10. 


Practical Astronomy, by W. Schroeder, 
Philosophical Library, New York, 1957, 
206 pp. $6. 
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P-E optical range instrumentation 
provides vital missile data 


checks behavior from production through flight 


Engineers must see what is happen- 
ing to a missile at every stage of its 
life. And Perkin-Elmer optical range 
instruments help them see .. . on the 
production line where stable plat- 
forms are aligned by P-E theodolites 
... at the launching site where P-E 
instruments first lay the missile for 
flight and then record vital data about 
its ascent ... downrange where P-E 
tracking instruments let observers 
watch — and learn from — its meteoric 
journey across the skies. 

Data grows more complex as the 
state of the art advances and the de- 
mands upon optical range instrumen- 
tation multiply. That calls for a high 
order of creative engineering. The 


number of P-E instruments in the — 


missile program today demonstrates 
the emphasis that Perkin- Elmer places 


on this singular talent. 


ROTI (Recording Optical Tracking Instrument) 
Mark | — Twin telescopes equipped with 
cameras track and record flight data on 
missiles and airborne objects. 

ROTI Mark Ii -— Single telescope and cam- 
era with velocity memory circuits for 
tracking objects that pass behind clouds 
or other obstructions. 

TPR (Telescopic Photographic Recorder) — A 
mobile unit similar in function to ROTI, 
but designed for easy mobility to any site 
accessible to prime mover. 

Kth57 Cine-Theodolite—Intermediate range 
instrument for evaluation of AA fire, 
monitoring bombing runs, etc. 

Azimuth Alignment Auto-Theodolite — Short, 
intermediate and long range models align 
inertial guidance systems at assembly 


Azimuth Alignment 
Theodolite 


Satellite Tracking 
Camera 


Sighting Telescope 


and at the missile launching site. - 
Sighting Telescope — Combines wide field 
of view and high magnification for quick 
target acquisition. 

SATRACK (Satellite Tracking Camera) — Em- 
ploys P-E aspheric optical systems. High 
light-gathering power, wide field of view 
will enable camera to photograph IGY 
satellites. 

Auto-Zoom* Lens— Extends versatility and 
range of standard closed circuit TV cam- 
eras. Permits remotely controlled wide- 
angle or closeup observation of action, 
target or instruments. 


An interesting booklet, “Optical Track- 
ing Instruments,” describes these and 
other P-E instruments for the Space Age 
more fully. Write for it. *T.M. 


ENGINEERING AND OPTICAL DIVISION 


NORWALK, 


Perkin-Elmer 
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The Visicorder has charted the orbit of Sputnik | 


A Model 906 Honeywell Visicorder Oscillo- type radio receivers. The beat oscillators gen- 
graph wrote this record of the signals from erated audio output signals, a semi-conductor _ 
Sputnik I for the Department of Electrical Engi- bridge curcuit rectified them, and the d-c output, 
neering at the University of Illinois at Urbana. filtered by an R-C network with atime constant —© 
The marginal notes are those of Edgar Hayden, of about .003 seconds, was used to drive the 
the research associate who took the record. Visicorder galvanometers directly. . 

Interferometer-type antenna systems (2 dipole The Visicorder, teamed with the interferome- — 
elements 4 wave length above ground spaced sev- ter antenna, quickly established a record of the © 
eral wavelengths along a north-south baseline) orbit of Sputnik I. 


is a record of Sputnik I | 


The HONEYWELL VISICORDER is the first high- | 
frequency, high-sensitivity direct recording oscil- 
lograph. In laboratories and in the field every- 
where, instantly-readable Visicorder records are 
pointing the way to new advances in product 
design, rocketry, computing, control, nucleonics 
...in any field where high speed variables are — 
under study. 

To record high frequency variables—and mon- 
itor them as they are recorded—use the Visi- _ 
corder Oscillograph. Call your nearest Minne- _ 
apolis Honeywell Industrial Sales Office for a _ 
demonstration. 
Reference Data: Write for Visicorder Bulletin 
Minneapolis Honeywell Regulator Co., 
Industrial Products Group, Heiland Division 


5200 E. Evans Ave., Denver 22, Colo. hae 
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Specialists 
wanted for 


Nuclear 
Powerplants 


1. Senior Engineer Research or 
Engineering: Nuclear physicist, 
PhD. (or equivalent experience) 
for basic studies in the application 
of nuclear power to air vehicles. 


2. Senior Engineer Research: 
Nuclear Engineer, M.S. degree or 
equivalent, with a background in 
unconventional powerplant design. 


3. Research Engineer: B. S. degree 
in nuclear engineering. College 
work in general powerplant design. 
4. Engineering Specialist or Sen- 
ior Engineer Research: M.S. or 
equivalent, with a strong back- 
ground in hydrodynamics and 
fluid mechanics to study basic re- 
lationships of propulsion equip- 
ment for air vehicles. 

5. Senior Engineer Research or 
Engineering Specialist: PhD. de- 
sired with a background in nuclear 
energy, servo mechanism, controls 
or instrumentation. Should have 
preferably 4 years’ additional 
experience. 

6. Senior Engineer Research: 
PhD. or M.S. degree with a back- 
ground in heat transfer relating 
to gas or liquid cooled nuclear 
power plants, preferably of air- 
craft type. 

7. Senior Structural Engineer or 
Stress Analyst: Several years’ ex- 
perience in mechanical stress and 
structural analysis. M.S. or B.S. 
degree or equivalent. 


Write: Mr. A.D. Jamieson, 
Rocketdyne Engineering 
Personnel Dept. 

6633 Canoga Ave., Canoga Park, 
California 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 


BUILDERS OF POWER FOR OUTER SPACE 
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Technical Literature Digest___ 


and Rocket 
Propulsion Engines 


The Design and Testing of an Axial 
Compressor Having a Mean Stage Tem- 
perature Rise of 30 deg. C, by A. D. 5S. 
Carter, S. J. Andrews and E. A. Fielder, 
Gt. Brit., Aeron. Res. Council, Rep. & 
Mem. 3985 (formerly A.R.C. Tech. Rep. 
16775; Nat. Gas Turbine Estab., Rep. 
R148), 1957, 14 pp. 

A Survey of the Alternating Pressures 
Exciting High Frequency Vibrations in Gas 
Turbines, by J. R. Forshaw, Gt. Brit., 
Aeron. Res. Council, Rep. & Mem. 2989 
(formerly A.R.C. Tech. Rep. 17133; 
Nat. Gas Turbine Estab., Rep. R148), 
1957, 20 pp. 

Design and Performance of an Experi- 
mental Axial Discharge Mixed Flow Com- 
pressor, III: Over All Performance of 
Impeller and Supersonic Diffuser Com- 
bination, by Ward W. Wilcox and William 
H. Robbins, NACA RM E51A02, April 
1951, 26 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
119, p. 16, 8/16/57.) 

Performance Investigation of Can Type 
Combustor from British Turbojet Engine 
and Comparison of This Combustor with 
Can Type Combustor from U. S. Turbojet 
Engine, by William P. Cook and Richard 
G. Koch, NACA RM E9H29, Feb. 1950, 
25pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 119, 
p. 16, 8/16/57.) 

Approximation for the Effect of Twist on 
the Vibration of a Turbine Blade, by 
A. I. Martin, Aeron. Quart., vol. 8, Aug. 
1957, pp. 291-308. 

Aero-engine Survey, 


Aeroplane, vol. 


| 93, Sept. 6, 1957, pp. 376-380 


Ignition Systems Must Keep up with 
Engine Performance, by Kurt R. Stehling, 
Aviation Age, vol. 28, Sept. 1957, pp. 
56-61. 

Liquid-Rocket Engine Control, by Ru- 
dolph H. Reichel, Missiles and Rockets, 
vol. 2, Sept. 1957, pp. 101-102, 104, 106. 

Propulsion Systems Evaluation, by 
George P. Sutton, Missiles and Rockets, 
vol. 2, Sept. 1957, pp. 123-124, 126. 


Aerodynamics of Jet 
Propelled Vehicles 


Aspects of High Energy Ballistics, by 
Louis Gold, J. Franklin Inst., vol. 264, 
Oct. 1957, pp. 301-311. 

Criteria for the Generalized Dynamic 
Stability of a Rolling Symmetric Missile, 
by Charles H. Murp' tg , J. Aeron. Sci., 
vol. 24, Oct. 1957, pp. od 774. 

A Comparative Analysis of the Per- 
formance of Long-range Hyper-velocity 
Vehicles, by Alfred J. Eggers Jr., H. 
Julian Allen and Stanford E. Neice, 
NACA TN 4046, Oct. 1957, 66 pp. 

A Study of the Motion and Aerodynamic 
Heating of Missiles Entering the Earth’s 
Atmosphere at High Supersonic Speeds, 


Ep1tor’s Note: Contributions from E. R. 
G. Eckert, J. P. Hartnett, T. F. Irvine Jr. and 
P J. Schneider of the Heat Transfer Labora- 
tory, University of Minnesota, are gratefully 
acknowledged. 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


by H. Julian Allen and A. J. Eggers Jr., 
NACA TN 4047, Oct. 1957, 61 pp. 
(supersedes RM A53D28). 

Motion of a Ballistic Missile Angu 
larly Misalined wih the Flight Path fines on 
Entering the Atmosphere and Its Effect 
Upon Aerodynamic Heating, Aerodynamic 
Loads, and Miss Distance, by H. Julian 
Allen, NACA TN 4048, Oct. 1957, 66 pp. 
(supersedes RM A56F15). 

On the Nature of the Sonic Line for 
Supersonic and Hypersonic Flow over 
Blunt Bodies, by Ronald F. Probstein, 
Wright Air Dev. Center, TN 57-349 
(ASTI A-A D 142037), Sept. 1957, 17 pp. 


Heat Transfer and 
Fluid Flow 


Simplified Method for Calculating the 
Heating Surfaces at Low Combustion 
Product Temperatures, by Jiri Possner, 
Strojirenstvi, vol. 7, Feb. 1957, pp. 107-111 
(in Czechoslovakian ). 

Effect of Agitation on the Critical Tem- 
perature Difference for a Boiling Liquid, 
by F. S. Pramuk and J. W. Westwater, 
Chem. Engng. Progr., Symposium Series, 
vol. 52, no. 18, 1956, pp. 79-83. 

Ebullition from Solid Surfaces in the 
Absence of a Pre-existing Gaseous Phase, 
Trans. ASME, vol. 79, May 1957, pp. 
735-740. 

Heat Transfer to Flowing Gas-solids 
Mixtures in a Circular Tube, by Leonard 
Farbar, and M. J. Morley, Ind. Engng. 
Chem., vol. 49, July 1957, pp. 1143-1150. 

Exchange Processes in Steam Contain- 


ing Air, by F. von Bosnjakovic, Kyllek- 
nisk Tidskrift, no. 3, 1957, pp. 142-146 
(in Swedish ). 


Calculating Three-phase Flow in Hori- 
zontal Pipe Lines, by D. P. Sobocinski 
and R. L. Huntington, Pipe Line Ind., 
vol. 6, no. 1, Jan. 1957, pp. 43-45. 

Plane Couette Flow at Low Mach 
Number According to the Kinetic Theory 
of Gases, by Hsun-Tias Yang and Lester 
Lees, GALCIT Hydrodynamics Lab., Rep. 
M-36, Feb. 1, 1957, 37 pp. 

Heat Transfer in Pipe Flow at High 
Speeds, by J. E. Bialokoz and O. A. 
Saunders, Proc. Inst. Mech. Engng, 
vol. 170, no. 12, 1956, pp°389=406. 

Study of Heat Transfer in a Spiral 
Canal, by V. G. Fastovskii and A. EF. 
Rovinskii, Jeploenergetika, vol. 4, Jan. 
1957, pp. 39-41 (in Russian). 

On the Pressure Losses for Turbulent 
Flow in Smooth Pipe Bends, by H. Ito, 
Tohoku University, Inst. High Speed 
Mech., Rep., vol. 7, 1956, pp. 149-174. 

Capillary Two-phase Flow, by H. A. 
Whitesel, Refriger. Engng., vol. 65, April 
1957, pp. '42- 44, 98-99. 

One-dimensional Transient Flow in 

Pipe with Two Gases, by R. S. Benson, 
Engng., vol. 202, Nov. 16, 1956, pp. 687- 
691. 
Pipe Friction Loss at High Pressures, 
by J. G. Slater, J. R. Villemonte, and 
H. J. Day, Proc. Am. Soc. Civil Engrs., 
vol. 83 (J. Hydraulics Div.), Feb. 1957 
(Paper 1163), 21 pp. 

Some Plane Problems in the Theory of 
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The new Du Pont Lino-Writ Rapid Process- Two big extras from Du Pont: quicker proc- 
ing Chemicals Kit enables you to process _ essing than ever before... without the disad- 
Lino-Writ 4 faster at temperatures up to vantage of offensive odor. Ask your dealer 
135 degrees. And we’ve done away with about the new Du Pont Lino-Writ Rapid 
the disagreeable odor previously associated Processing Chemicals Kit. He’s stocking 
with high-temperature processing. Now you __it now. 

can develop up to 20 feet per minute, and 


still get the same sharp, clear results that E. I. du Pont de Nemours & Co. (Inc.), 
Lino-Writ is famous for. That means that on —_ Photo Products Department, Wilmington 98, 
some runs you can now process your paper Delaware. In Canada: Du Pont Company of 
twice as fast as you previously could. oe Canada (1956) Limited, Toronto. _ 
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APL OFFERS 
GREATER 

FREEDOM OF 
ACTIVITY 


APL has responsibility for 
the technical direction of 
much of the guided missile 


__-: program of the Navy Bureau 


of Ordnance. As a result 
staff members participate in 
assignments of challenging 
scope that range from basic 
research to prototype test- 
ing of weapons and weapons 
systems. 


A high degree of freedom 


members to give free rein 
to their talents and ideas. 
Thus, professional advance- 
ment and opportunities to 
accept program responsi- 
bility come rapidly. Promo- 
_ tion is rapid, too, because 

of our policy of placing 
professional technical men 
at all levels of supervision. 


APL’s past accomplish- 
ments include: the first ram- 
jet engine, the Aerobee high 
altitude rocket, the super- © 
sonic Terrier, Tartar, and 
Talos missiles. Presently the 
Laboratory is engaged in 
solving complex and ad- 
vanced problems leading to 
future weapons and weap- 


ons systems vital to the ~~ 
y 


National security. Interested 
engineers and physicists are 
? invited to address inquiries 


Professional Staff 
Appointments 


The Johns Hopkins University 
Applied Physics Laboratory 


8617 Georgia Avenue, Silver Spring, Md. 


of action enables APL staff | 


Thermal Waves, by B. G. Korenev, 
Doklady Akademii Nauk SSSR, vol. 112, 
Jan.-Feb. 1957, pp. 29-32 (in Russian). 
Thermodynamical Theory of Thermal 
Conduction of Dielectrics under Electric 
Fields, by S. Mascarenhas, Nuovo Cimento, 
vol. 5, no. 4, May 1957, pp. 1118-1121. 


Analog Computers Calculate Heat Trans- 
fer, by Robert 8. Schechter, Petroleum 
Refiner, vol. 36, Feb. 1957, pp. 112-114. 


Difference Methods for the Calculation 
of Temperature Processes in One-dimen- 
sional Heat Flow in Bodies, by Herbert 
Schuh, VDI Forschungsheft, no. 459, 
1957, 43 pp. 

The Solutions of Boundary Value Prob- 
lems on a AC Analog Computer, 
by M. Yanowitch, Computers and Auto- 
mation, vol. 6, Feb. 1957, pp. 26-29. 


Half Values as Units of Measurement 
in Steady State Exchange, by F. von 
Bosnjakovic, Chemic-Ingenieur-T echnik, 
vol. 29, no. 3, 1957, pp. 187-197. 

Experimental Study of the Velocity and 
Temperature Distribution in a High- 
velocity Vortex-type Flow, by J. P. Hart- 
nett and E. R. G. Eckert, Trans. ASME, 
vol. 79, May 1957, pp. 751-758. 

On the Formulation and Analysis of 
Numerical Methods for Time-Dependent 
Transport Equations, by Herbert B. Keller 
and Burton Wendroff, New York Univ., 
Inst. Math. Sci., (U. S. Atomic Energy 
Comm.), NYO-7694, Dec. 1956, 38 pp. 

Turbulent Mixing of Plane Hot Air 
Flows, by W. Szablewski, /ng. Arch., 
No. 1, 1957, pp. 10-25 (in German). 

Aerothermal Seminar of the Faculty of 
Sciences of Paris, France, Ministere de 
UV Air, Note Tech. 63, 1957, 110 pp. (in 
French). 

Thermal Transpiration: Application of 
Liang’s Equation, by M. J. Bennett and 
F. C. Tompkins, Trans. Faraday Soc., 
vol. 53, Feb. 1957, pp. 185-192. 

Determination of Thermal Conductivi- 
ties of Metals by Measuring Transient 
Temperatures in Semi-infinite Solids, 
by S. T. Hsu, Trans. ASME, vol. 79, July 
1957, pp. 1197-1202. 

Analysis of Heat-driven Oscillations of 
Gas Flows, by H. J. Merk, Appl. Sci. 
Res., Sec. A., vol. 6, no. 4, 1957, pp. 317- 
336. 

What Are Heat Transfer Characteristics 
of Wire and Tube Condensers? by O. W. 
Witzell and W. E. Fontaine, Refrig. 
Engng., vol. 65, March 1957, pp. 33-37, 
127. 


Heat Transfer Characteristics of Liquid 
Metal in Closed Thermosyphon, by F. J. 
Bayley and N. Bell, Engng. vol. 183, 
March 8, 1957, pp. 300-302. 

The Use of Sodium and of Sodium- 
Potassium Alloy as a Heat-transfer Me- 
dium, by W. B. Hall and T. I. M. Crofts, 
Proc., Inst. Mech. Engrs., vol. 170, no. 10, 
1956, pp. 281-305. . 

Temperature and Velocity Distribution 
and Transfer of Heat in a Liquid Metal, 
by H. E. Brown, B. H. Amstead, and B. E. 
Short, Trans., ASME, vol. 79, Feb. 1957, 
pp. 279-285. 

Utilization of Solar Furnaces in High- 
temperature Research, by Pol Duwez, 
Trans. ASME, vol. 79, July 1957, pp. 
1014-1018. 

Research on Application of Cooling to 
Gas Turbines, by J. B. Esgar, J. N. B. 
Livingood, and R. O. Hickel, Trans. 
ASME, vol. 79, April 1957, pp. 645-651. 

Analysis of Large Aperture Parabolic 
Mirrors for Solar Furnaces, by Joseph 
Farber and Burton I. Davis, /. Optical 


i 


Soc. America, vol. 47, March 1957, pp. 
216-220. 

The Cooling Process in Periodic Equip- 
ment. Heat Exchangers, Economic Size 
of the System. Operating Costs, by Ugo 
Fasoli, Chimica e l’Industria, vol. 39, Feb. 
1957, pp. 83-92 (in Italian). 

Heat Transfer to Vibrating Air Columns, 
by W. B. Harrison, W. C. Boteler, T. W. 
Jackson, A. Lowi, Jr., and F. A. 
Thomas, Jr., NACA (unpublished), (N- 
94857 ), 1955, 38 pp. 

Measurements on the Velocity of Sound 
in Nitrogen under High Pressure, by A. 
van Itterbeek, W. de Rop and G. Forrez, 
Applied Sci. Res., vol. 6, Sect. A, no. 5-6, 
1957, pp. 421-432. 

On the Use of Matrices for Solvin: 
Periodic Heat Flow Problems, by G. W. 7 
White, Applied Sci. Res., vol. 6, Sect. A 
no. 5-6, 1957, pp. 433-444. 


On the Existence of Periodic Waves 
Near Critical Speed, by Walter Littman 
Comm. Pure Appl. Math., vol. 10, Ma: 
1957, pp. 241-269. 


Deposition of Suspended Particles fron 
Turbulent Gas Streams, by S. K. Fried- 
lander and H. F. Johnstone, Ind. Engng 
Chem., vol. 49, July 1957, pp. 1151-1156. 


Electrical Conductivity of Thermally 
Ionized Air Produced in a Shock Tube, b) 
Lawrence Lamb and Shao-Chi Lin, J 
Appl. Phys., vol. 28, July 1957, pp. 754- 
759. 


Shock Hugoniot of Air, by W. E. Deal, 
J. Appl. Phys., vol. 28, July 1957, pp. 
782-784. 

Quantum-mechanical Effects in Gas 
Dynamics, by E. V. Laitone, J. Chem. 
Phys., vol. 26, June 1957, pp. 1560-1566. 

Thermal Conductivity of Gas Mixtures 
in Chemical Equilibrium, by James N. 
Butler and Richard 8. Brokaw, J. Chem. 
Phys., vol. 26, June 1957, pp. 1636-1642. 


An Electronic Analogue for Supersonic 
Flow, by Leslie S. G. Kovasznay, J. 
Fluid Mech., vol. 2, June 1957, pp. 383- 
396. 

On the Stability of Plane Shock Waves, 
by N. C. Freeman, J. Fluid Mech., vol. 2, 
June 1957, pp. 397-411. 

Simplified Equations for Calculating 
Local and Total Heat Flux to Noniso- 
thermal Surfaces, by E. R. G. Eckert, 
J. P. Hartnett, and Roland Birkebak, 
J. Aeron. Sci., vol. 24, July 1957, pp. 549- 
551. 

The Relation between the Flow of Non- 
Newtonian Fluids and Turbulent New- 
tonian Fluids, by R.S. Rivlin, Quarterly 
Appl. Math., vol. 15, July 1957, pp. 212- 
215. 

Intermolecular Free Length in Liquids 
in Relation to Sound Velocity, by Otohilo 
Nomoto, J. Physical Soc. Japan, vol. 12, 
June 1957, p. 739. 

The Direct Measurement of Local Skin 
Friction on Aerobee Hi Rockets in Free 
Flight, by Felix W. Fenter and W. C. 
Lyons, Jr., Texas Uni., Defense Res. 
Lab. CM-877 (DRL 391), May 1957, 
47 pp., 27 fig. 

The Evaporation of Fuel Sprays, I. 
Theorectical Treatment, by G. W. Benson, 
Canada, Nat. Aeron. Estab., Lab. Rep. 
LR-181, (ASTIA AD _ 121803), Nov. 
1956, 38 pp. 

States of Partial Equilibrium behind 
Shock Waves in O., CO, COs, and 
by Wayland C. Griffith and Anne Kenny, 
Princeton Univ., Dept. Phys., Tech. Rep. 
II-23, March 1957, 16 pp., 6 tab., 6 fig. 

Measurements of the Thrust Produced 


by Convergent-Divergent Nozzles at Pres- 
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The burst of extra energy that lofted the 
Army's Jupiter-C was generated by a new liquid 
fue! mixture based on UDMH. This unique 
first-stage fuel, which was developed by the 
Rocketdyne Division of North American Aviation, 


Inc., was formulated and supplied by Westvaco. 


We welcome the opportunity to discuss DIMAZINE 


and its properties with your staff. 


U.S. Army Photograph 
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sure Ratios up to 20, by P. F. Ashwood, 
G. W. Crosse and Jean E. Goddard, Gt. 
Brit., Aeron. Res. Council, Current Paper 
326, 1957, 15 pp., 10 fig. (formerly ARC 
Tech. Rep. 18,771; Gt. Brit., Natl. Gas 
Turbine Estab. Mem., M.288). 

Heat Transfer Research Holds Key to 
Better Rocket Cooling, by Kurt R. Steh- 
ling, Aviation Age, vol. 28, Aug. 1957, 
pp. 36-43. 

Mechanical Design Analysis of Several 
Noncritical Air-Cooled Turbine Disks and 
a Corrugated-Insert Air-Cooled Turbine 
Rotor Blade, by Merland L. Moseson, 
Morton H. Krasner and Robert R. Zie- 
mer, NACA RM E53E21, July 1953, 
49 pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 119, 
p. 19, 8/16/57.) 

Investigation of Water-spray Cooling 
of Turbine Blades in a Turbo-jet engine, 
by John C. Freche and William J. Stelp- 
flugy,; NACA RM E53A23, March 1953, 
45 pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 119, p. 
19, 8/16/57.) 

Analytical Comparison of Turbine-blade 
Cooling Systems Designed for a Turbojet 
Engine Operating at Supersonic Speed and 
High Altitude, I: Liquid-cooling Systems, 
by Wilson B. Schramm, Alfred J. Nachti- 
gall and Vernon L. Arne, NACA RM 
E52J29, Jan. 1953, 52 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abs. 119, p. 18, 8/16/57.) 

Analytical Comparison of Turbine-blade 
Cooling Systems Designed for a Turbojet 
Engine Operating at Supersonic Speed 
and High Altitude, II: Air-cooling Systems, 
by Wilson B. Schramm, Vernon L. Arne 
and Alfred J. Nachtigall, NACA RM 
E52J30, Jan. 1953, 50 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abs. 119, p. 18, 8/16/57.) 

An Analytical Method for Evaluating 
Factors Affecting Application of Transpira- 
tion Cooling to Gas Turbine Blades, by 
Jack B. Esgar, NACA RM _ E52G01, 
Sept. 1952, 68 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
119, p. 17, 8/16/57.) 

An Experimental Cascade Study of the 
Effects of a Solidity Reduction on the Two- 
Dimensional Aerodynamic Characteristics 
of a Turbine-rotor Blade Suitable for 
Air Cooling, by Henry W. Plohr and Wil- 
liam J. Nusbaum, NACA RM E52B27, 
May 1952,17 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abs. 119, p. 17, 8/16/57.) 


Cold-air Investigation of a Turbine with 
Nontwisted Rotor Blades Suitable for 
Air Cooling, by Thomas R. Heaton, Wil- 
liam R. Slivka and Leonard F. Westra, 
NACA RM E52A25, March 1952, 34 
pp. (Declassified from Confidential by 
authority of NACA Res. Abs. 119, p. 17, 
8/16/57.) 

Aerodynamic Studies in the Shock Tube, 
by Josef Rabinowicz, Calif. Inst. Tech., 
Hypersonic Res. Project, Mem. 38, June 
10, 1957, 111 pp., 38 ref. 

Pressure Recovery, Drag, and Subcriti- 
cal Stability Characteristics of Three Coni- 
cal Supersonic Diffusers at Stream Mach 
Numbers from 1.7 to 2.0, by Theodore J. 
Nussdorfer, Leonard J. Obery and Gerald 
W. Englert, NACA RM E51H27, Feb. 
1952, 31 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
119, p. 17, 8/16/57.) 

Pressure Recovery, Drag, and Subcriti- 
cal Stability Characteristics of Conical 
Supersonic Diffusers with Boundary- 
layer Removal, by Leonard J. Oberey, 
Gerald W. Englert and Theodore J. 
Nussdorfer, NA 


CA RM E51H29, Feb. — 


5 
1952, 25 pp. (Declassified from Con- 


fidential by authority of NACA Res. 
Abs. 119, p. 17, 8/16/57.) 

Development and Application of a 
Technique for Steady State Aerodynamic 
Heat Transfer Measurements, by Frederic 
W. Hartwig, Calif. Inst. Tech., Hypersonic 
Res. Project, Mem. 37, June 1, 1957, 48 
pp., 13 ref. 

Formulas for Thermal Conductivity of 
Ternary Gas Mixtures, by B. N. Srivas- 
tava and §. C. Saxena, J. Chem. Phys., 
vol. 27, Aug. 1957, pp. 583-584. 


The Laminar and Turbulent Mixing of 
ets of Compressible Fluid, Part I. Flow 
‘ar from the Orifice, by L. J. Crane and 

D. C. Pack, J. Fluid Mech., vol. 2, July 
1957, pp. 449-455. 


The Secondary Flow in a Cascade of 
Turbine Blades, by W. D. Armstrong, 
Gt. Brit., Aeron. Res. Council, Rep. & 
Mem. 2979 (formerly A.R.C. Tech. Rep. 
17442), 1957, 14 pp. 

Radiant Heat Exchange in a Gas Filled 
Enclosure, by H. C. Hottel and E. &. 
Cohen, Project Squid Tech. Rep., MIT-14- 
P (ASTIA AD 138820), 1957 (available 
only on microcards). 

Results of Systematic Investigations 
on Secondary Flow Losses in Cascades, 
Part I: Secondary Flow Losses in Com- 
pressor Cascades of Profile NACA 8410, 
by K. Gersten, Braunschweig, Tech. Hoch- 
schule, Inst. fiir Stromungsmechanik, Rep. 
54/29a, June 1955, 32 pp. 

Equations, Tables, and Figures for Use 
in the Analysis for Helium Flow at Super- 
sonic and Hypersonic Speeds, by James N. 
Mueller, NACA TN 4063, Sept. 1957, 
178 pp. 

Mollier Chart for Air in Dissociated 
Equilibrium at Temperatures of 2000°K 
to 15000°K, by I. Korobkin and 8. M. 
Hastings, NAVORD Rep. 4446 (Aero- 
ballistics Res. Rep. 368), May 1957, 7 pp. 


Analysis of Shock Motion in Ducts dur- 
ing Disturbances in Downstream Pres- 
sure, by Herbert G. Hurrell, NACA TN 
4090, Sept. 1957, 11 pp. 

An Analysis of the Interaction of Shock 
Waves with Unseparated Turbulent 
Boundary Layers, by Andrew G. Hammitt, 
Princeton Univ., Dept. Aeron. Engng., 
Rep. 371 (AFOSR-TN-57-108; ASTIA 
AD 120460), Feb. 1957, 14 pp. 


A Basic ee’ to Shock Front 
Analysis, by Ph. J. Theodorides, Univ. 
Maryland, Inst. Fluid Dynamics Appl. 
Math., TN BN-93 (AFOSR-TN-57-59; 
ASTI A-AD 115099), Jan. 1957, 31 pp. 


The Structure of a Shock Wave in a Gas 
Having a Long Relaxation Time, by L. 
Talbot, Univ. Calif., Inst. Engng. Res., 
Tech. Rep. HE-150-145, March 1957, 23 
pp. 

Some Notes on Shock Wave Boundary 
Layer Interactions and on the Effect of 
Suction on the Separation of Laminar 
Boundary Layers, by K. N. C. Bray, 
Gt. Brit., Aeron. Res. Council, Current 
Paper 332 (formerly A.R.C. Tech. Rep. 
18631; Handley Page Ltd., Aero. Res. 
Rep. 24), 1957, 7 pp. 

A Phenomenal Theory of the Free 
Turbulent Flow of a Non-Reacting Binary 
Gas Mixture, by Luigi G. Napolitano, 
Brooklyn Polytech. Inst., Dept. Aeron. 
Engng. Appl. Mech., Rep. 317 (AFOSR- 
TN-57-39; ASTIA AD 115077), Dec. 
1956, 46 pp. 

Turbulent Mixing of Two Non-Reacting 
Gases, by Luigi G. Napolitano, Brooklyn 
Polytech. Inst., Dept. Aeron. E: . Appl. 
Mech., Rep. 323 (APOSR-TN-57-41: 


_ ASTIA AD 115079), Jan. 1957, 29 pp. 
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... cross new 
frontiers in system 
electronics at THE 
GARRETT CORPORATION 


Increased activity in the 
design and production of sys- 
tem electronics has created 

openings for engineers in the 

following areas: 


ELECTRONIC AND AIR DATA Bee 
SYSTEMS Required are men of proj- 
ect engineering capabilities. Also 
required are development and 
design engineers with specialized 
experience in servo-mechanisms, 
circuit and analog computer design 
utilizing vacuum tubes, transistors, 
and magnetic amplifiers. 


SERVO-MECHANISMS 

AND ELECTRO-MAGNETICS Complete 
working knowledge of electro-magnet- 
ic theory and familiarity with materi- 
als and methods employed in the design 
of magnetic amplifiers is required. 


FLIGHT INSTRUMENTS AND 
TRANSDUCER DEVELOPMENT 
Requires engineers capable of analyz- 
ing performance during preliminary 
design and able to prepare proposals 
and reports. 


FLIGHT INSTRUMENTS 

DESIGN Requires engineers skilled 
with the drafting and design of light 
mechanisms for production in which 
low friction, freedom from vibration 
effects and compensation of thermo 
expansion are important. 


HIGH FREQUENCY MOTORS, 
GENERATORS, CONTROLS Requires 
electrical design engineers with 
BSEE or equivalent interested in 
high frequency motors, generators 
and associated controls. 


Send resume of education 
and experience today to: 
Mr. G. D. Bradley 


9851 S. Sepulveda Blvd. 
Los Angeles 45, Calif. 
DIVISIONS: 
AiResearch Manufacturing 
Los Angeles 
AiResearch Manufacturing 
Phoenix 
AiResearch Industrial 
Rex—Aero Engineering 
Airsupply — Air Cruisers 
AiResearch Aviation 
Service 
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cortre 


Ratio pressure reg $ 
16 


18 types 


lute pressure regulator 


ba 


AIRCRAFT 


Shutoff volves 
349 typ 


Differential pressure 
regulators 
414 types 


AiResearch is the largest 
designer and manufacturer of 
pneumatic controls for the air- 
craft and associated industries. 
During the past 10 years more 
than 300,000 units have been 
produced and are in service. 


Systems, and AIRCRAFT, MISSILE, ELECTRONIC, NUCLEAR AND INDUSTRIAL APPLICATIONS 


Apri 1958 


Compartment regulators 


mperatyre control 


Pneumatic switches 
10 types 


Temperatures of the fluids 
(including gas and liquids) 
range from —400°F to 
+2000°F at pressures to 
+6000 psig. The units operate 
at any ambient pressure at 
ambient temperatures from 


Flow regulators 
28 types 


Mixing vaives 
25 types 


Pneumatic engine 


6 types 


computer-controllers 


—300°F to over +1000°F. 
Line diameters range from % 
inch to 15 inches. 

This equipment is developed 
and tested in the finest pneu- 
matic facilities in the world. 
Your inquiries are invited. 


ENGINEERING REPRESENTATIVES: AIRSUPPLY AND AERO ENGINEERING, OFFICES ‘MAJOR 


AiResearch Manufacturing Divisions 


Los Angeles 45, California « Phoenix, Arizona 
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ia More than 20 models available 


covered by patents 


HIRTEEN BRISTOL 
HIGH-SPEED RELAYS 
IN THIS CONVERTER! 


Twelve-and-a-half microvolt resolution at 20 readings per 
second! That’s the outstanding feature of the analogue-to- 
digital converter, developed by Non-Linear Systems, Inc., 
Del Mar, California, to “digitalize” the output of low-volt- 
age transducers in either ground or airborne service. 

It’s significant that Non-Linear Systems engineers se- 
lected thirteen miniature Bristol Syncroverter* high-speed 
relays (inset, top) for use in the converter scanning cir- 
cuits. This versatile, high-speed, polarized relay has earned 
an enviable reputation for reliability, long life and immu- 
nity to shock and vibration in just such critical low-level, 
dry-circuit applications. 


Are dry circuits your problem? 


If so, we believe we have the answer. Dry-circuit reliability 
and long life are outstanding features of the Syncroverter 
high-speed relay. It’s unaffected during severe shock and 
vibration. It has fast pull-in and drop-out and negligible 
contact resistance, and it operates reliably over a wide tem- 
perature range. 


You can specify Bristol Syncroverter high-speed relays 


in an extremely wide variety of operating characteristics 
and in various case and mounting arrangements. Ask us 
for complete details. Write: The Bristol Company, 175 
Bristol Road, Waterbury 20, Conn. 8.34 
*T. M. Reg. U. S. Pat. Of. 


= BRISTO FINE PRECISION INSTRUMENTS 
FOR OVER 68 YEARS Fos 


Approximation Methods in Compres- 
sible Fluid Dynamics, Part I, Chapters 
I-IV; Part II, Chapters V—VII, by Isao 
Imai, Univ. Maryland, Inst. Fluid Dy- 
namics Appl. Math., TN BN-95 (AST/A 
AD 128411), March 1957, 184 pp. 

Further Investigations of the Laminar 
Turbulent Transition in a Free Jet 
(Annular Nozzle), AFOSR-T N-57-31 (AS- 
TIA AD 126494), Dec. 1956, 22 pp. 


Statistical Studies on the Turbulent 
State of Motion, Part I: Kinematics, 
by Hans Gebelein and Erik Svenson, 
AFOSR-TN-57-278 (ASTI A AD 12653°), 
1957, 42 pp. 

A Hydrodynamic Study of Small Di- 
ameter Orifices in Thin Plates, by Henry 
L. Herman and Milton V. Ashley, Pica- 
tinny Arsenal, Samuel Feltman Ammuni- 
tion Labs., Tech. Rep. 2450, Aug. 1957, 
19 pp. 


Combustion, Fuels and 
Propellants 


Combustion Performance Evaluation of 
Magnesium Hydrocarbon Slurry Blends 
in a Simulated Tail Pipe Burner, b 
Leonard K. Tower and J. Robert Bra: 
stetter, NACA RM E51C26, May 1951 
53 pp. (Declassified from Confidentia 
by authority of NACA Res. Abs. 116, | 
14, 7/5/57.) 

Effect of Particle Size and Stabilizin; 
Additives on the Combustion Propertie 
of Magnesium Slurry, by Albert M. Lor 
and Vernida E. Evans, NACA RAJ 
E52K12, Jan. 1953, 28 pp. (Declass 
fied from Confidential by authority « 
NACA Res. Abs. 116, p. 15, 7/5/57.) 

An Experimental Investigation of the 
Combustion Properties of a Hydrocarbon 
Fuel and Several Magnesium and Boror 
Slurries, by Albert M. Lord, NAC. 
RM E52B01, April 1952, 30 pp. (De 
classified from Confidential by authority 
of NACA Res. Abs. 116, p. 14, 7/5/57.) 

Experimental Investigation of Physical 
and Combustion Properties of Several 
Residual Fuel Oils and Magnesium-fuel- 
oil Slurries in a Ram-jet-type Combustor, 
by Preston N. Cook, Jr., Vernida E. Evans 
and Erwin A. Lezberg, NACA RM 
E53D30, June 1953, 23 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 116, p. 15, 7/5/57.) 

Magnesium-Slurry Combustion Per- 
formance in 6/5-inch-diameter Ram-jet 
Engine Mounted in Connected-pipe Fa- 
cility, by J. Robert Branstetter, James B. 
Gibbs and Warner B. Kaufman, NACA 
RM E53E27, Aug. 1953, 63 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 116, p. 15, 7/5/57.) 

The Internal Ballistics of a Recoilless 
High-Low Pressure Gun, by J. N. Kapur, 
Appl. Sci. Res., vol. 6, Sect. A, no. 5-6, 
1957, pp. 445-466. 

Boron Fuel Production to Start, by J. S. 
Butz, Jr., Aviation Week, vol. 67, July 
15, 1957, p. 27. 

Decomposition of Alcohols over a 
Fischer-Tropsch Iron Catalyst, by S. R. 
Srinivasan and N. G. Basak, Fuel, vol 
36, July 1957, pp. 277-285. 

Supersaturation in Hydrocarbon Sys- 
tems, n-Pentane in the Liquid Phase, 
by W. B. Nichols, L. T. Carmichael, and 
B. H. Sage, Ind. Engng. Chem., vol. 49, 
July 1957, pp. 1165-1172. 

Decomposition of Hydrogen Peroxide 
Vapor on Relatively Inert Surfaces, by 
C. N. Satterfield and T. W. Stein, Ind. 
Engng. Chem., vol. 49, July 1957, pp. 
1173-1180. 
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GIANT CAPACITOR SYSTEM POWERS ADVANCED 


Lockheed experiments on electrical power 
sources for space craft are using electrical 
energy stored in a giant capacitor system 
one of the largest in the aircraft industry. 
This propulsion research —based on 


Missile Systems—where expanding electrohydrodynamics — is an example of the 
advanced work being conducted in Division 
missile programs insure 1: laboratories near Stanford University in 
=: the electrical, chemical and nuclear fields. 
-§ more promising careers 


The giant capacitor bank is capable of 

a 5 million kilowatt jolt of electricity, 

used to fire Lockheed’s new “hotshot” wind 
tunnel to 18,000 degrees Fahrenheit 

with speeds up to 15,000 miles an hour. 


Research in such sophisticated areas by 
Division scientists has earned Lockheed 
leadership in missiles. Nose cone reentry 
studies, for example, have helped speed 
delivery dates on Polaris, Lockheed’s 
submarine fleet ballistic missile. 


With expansion on our IRBM and other 
advanced programs, new positions have 
opened for engineers and scientists... 
positions that call for the utmost in 
creative ability. 


These areas include: Propulsion, _ 
Aerodynamics-Thermodynamics, | 
Guidance, Electronics, Flight 
Controls, Structures, Design. 7 


If you are qualified, you are invited to write 
M.W. Peterson, Research and Development 
Staff, Sunnyvale 20, California for | 
further information. 


MISSILE SYSTEMS 


A DIVISION OF LOCKHEED 
AIRCRAFT CORPORATION 


SUNNYVALE PALO ALTO VAN NUYS 


SANTA CRUZ* CALIFORNIA 


CAPE CANAVERAL, FLORIDA 


ALAMOGORDO, NEW MEXICO 


Discussing characteristics of the capacitor 
system are research scientist Dr. Edward Fisher 
of the Propulsion Staff, right, and Dr. Terry E 
Turner, Hot Shot Tunnel Group Leader. 

The two groups are engaged in a joint effort 

on advanced propulsion experiments. 
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_ Flame Propagation Rates, Chemical 
Nature of Attachment Surface, by Leon 
Lapidus, J. B. Rosen and R. H. Wilhelm, 
Ind. Engng. Chem., vol. 49, July 1957, pp. 
1181-1186. 

Derivation and Precision of a New Vapor 
Pressure Correlation for Petroleum Hy- 
drocarbons, by J. B. Maxwell and L. 8. 
Bonnell, Ind. Engng. Chem., vol. 49, 
July 1957, pp. 1187-1196. 

Theory of Multi-component Fluid Mix- 
tures, I: Statistical Order-Disorder Analy- 
sis in Multi-component Fiuid Mixtures, 
by Zevi W. Salsburg, Peter J. Wojtowicz, 
and John G. Kirkwood, J. Chem. Phys., 
vol. 26, June 1957, pp. 1533-1541. 


Entropy of Chemisorbed Hydrogen, 
by Tetsuo Takaishi and Akio Kobayashi, 
J. Chem. Phys., vol. 26, June 1957, pp. 
1542-1543. 

Infrared Spectra and Force Constants 
of Cyanoacetylene, by George C. Turrell, 
Walter D. Jones and Arthur Maki, J 
Chem. Phys., vol. 26, June 1957, pp. 1544— 
1548. 

Electronic Spectrum and Dissociation 
Energy of Fluorine, by A. L. G. Rees, 
J. Chem. Phys., vol. 26, June 1957, pp. 
1567-1571. 

Mechanism of the Gas Phase, Thermal 
Decomposition of Ozone, by Sidney W. 
Benson and Arthur E. Axworthy, Jr., 
J. Chem. Phys., vol. 26, June 1957, pp. 
1718-1719. 

Reaction of Carbon Monoxide and 
Ozone, by P. Harteck and S. Dondes, 
J. Chem. Phys., vol. 26, June 1957, pp. 
1734-1735. 

Preliminary Investigation of Propane 
Combustion in a 3-inch-diameter Duct at 
Inlet-air Temperatures of 1400° to 1600°F, 


“MONOBALL” 


Self-Aligning Bearings 


CHARACTERISTI¢ 


ANALYSIS 


" 3 Bronze Race and ew operating under normal loads 
with minimum friction requirements. 


Chrome Steel Ball 


RECOMMENDED USE 


Stainless Steel For types operating under high temper- 
Ball and Race ee (800-1200 degrees F.). . 
Chrome Alloy {Fee types operating under high radial 
Steel Ball and Race ultimate loads (3000-893,000 Ibs.). 


by Erwin A. Lezberg, NACA TN 4028, 
July 1957, 19 pp. 

Reaction of Fluorine with Carbon as a 
Means of Fluorine Disposal, by Harold 
W. Schmidt, NACA RM E5702, July 
1957, 17 pp. 

Propellant Vaporization as a Criterion 
for Rocket Engine Design; Calculations 
of Chamber Length to Vaporize a Single 
n-heptane drop, by Richard J. Priem, 
NACA TN 3985, July 1957, 41 pp. 

Combustion Instability in Liquid Rocket 
Propellant Motors, 20th Quarterly Prog- 
ress Report, Feb. 1-April 30, 1957, by 
Jerry Grey, Princeton Univ., Dept. Aeron. 
Engng., Rep. 216-t, June 1, 1957, 25 pp. 

Investigation of Jet Flame-holders, II, 
by A. W. Noon, Wright Air Dev. Center, 
TN 56-316, Part II, April 1957, 30 pp. 

Drop Burning Rates of Hydrocarbon 
and Nonhydrocarbon Fuels, by Arthur L. 
Smith and Charles C. Graves, NACA 
RM E57F11, Aug. 1957, 22 pp. 

Streamlined Method for Determining 
the Available Diphenylamine Content of 
Propellants, by W. Dewey Tullis and Mil- 
ton Roth, Picatinny Arsenal Rep. 2432, 
July 1957, 9 pp. 

Effect of Fuel Volatility on Altitude 
Starting Limits of a Turbo-jet Engine, 
by H. D. Wilstead and J. C. Armstrong, 
NACA RM E50G10, Sept. 11, 1950, 
27 pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 118, 
p. 7, 8/16/57.) 

Energy Limitations of Liquid Rocket 
Propellants, by J. F. Tormey, Aircraft 
Engng., vol. 29, Aug. 1957, pp. 248-251. 

Hydrogen Cyanide-Fluorine-Oxygen 
Flame, by Charles S. Stokes and A. V. 
Crosse, Jr., Ind. Engng. Chem., vol. 49, 


Aug. 1957, pp. 1311-1314. es 

Solid Fuel Industry Roundup, Missiles 
and Rockets, vol. 2, Aug. 1957, pp. 67-73. 

Swing to Solid Propellants, by John 
Wilson, Missiles and Rockets, vol. 2, Aug. 
1957, pp. 92- 94. 

Machining for Solid Propellant Rockets, 
by Joachim H. Kauffmann, Missiles and 
Rockets, vol. 2, Aug. 1957, pp. 105-111. 

Processing for Solids at Thiokol, 
Missiles and Rockets, vol. 2, Aug. 1957, 
pp. 112-113. 

The Hunt is on for New Sources of Pro- 
pulsion Energy, by Hugh Harvey, SAK 
J., vol. 65, Aug. 1957, pp. 17-20. 

An Unclassified Literature Survey on 
Combustion Pressure Oscillations 
Liquid Propellant Rocket Motors, by J. 
R. Osborn and M. J. Zucrow, Purdie 
Univ., Rocket Lat., Tech. Mem. 57-\, 
June 1957, 134 pp., 64 ref. (Contra:t 
N7-onr 39418). 

Experimental and Calculated Histories 
of Vaporizing Fuel Drops, by R. J. Priei, 
G. L. Borman, M. M. El Wakil, O. 
Uyehara and P. S. Myers, NACA 7.\ 
3988, Aug. 1957, 66 pp. 

Stability Limits and Burning Velocities 
for Some Laminar and Turbulent Pro- 
= and Hydrogen Flames at Reduced 

ressure, by Burton Fine, NACA 7. 
4031, Aug. 1957, 49 pp. 

Simulated High Altitude Combustion 
Research, \//7', Fuels Res. Lab., Quarter! 
Progress Rep. 6 (Contract AF33(6163037 
Sept. 1-Nov. 31, 1956, Jan. 11, 1957, 1) 
pp. (ASITA AD-122,931), (DSR-7361 

Investigation of Ionization in Flames, 
I. A Method for the Extraction of the 
Ions from a Flame, by A. Van Tiggelen, 
L. Vanreusal, and J. Deckers. II. Meas- 


GRIEVE-HENDRY 


Standard line of 1000 F. 
RECIRCULATING 


For Solution 
Heat Treating 
and Similar 

Processes 


With potentiom- 
eter control will 
meet MIL Speci- 
fications. 


Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


Full range of sizes 
for gas or 
electric heat. 


OTHER 
STANDARD MODELS 


s1105° 


and up 


Write for literature 


GRIEVE-HENDRY CO., INC. 


1410 W. Carroll Avenue, Chicago 7, Illinois 
Export Dept., 10406 S. Western Ave., Chicago 43, Illinois 
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three to make 
‘ready for outer space. 


AMMONIUM 
PERCHLORATE 


in th ec h emi st ry and AMMONIUM PERCHLORATE— American Potash & Chemical 


Chemical Engineering Achievement Award for pioneering work 


an d SO | j d p ro pel | ant in the atomic age metals, Trona continues in the forefront of basic 


ee Ox dants , compounds for high energy and missile applications. 


ms 8 Se 3000 West Sixth Street, Los Angeles 54, California 
, 99 Park Avenue, New York 16, New York 


: vat Offices: LOS ANGELES, NEW YORK, CHICAGO, SAN FRANCISCO, 
me PORTLAND (Ore.), ATLANTA, COLUMBUS (0.), SHREVEPORT 


* 
TRADENAME AND TRADEMARK OF AP&CC 


Wherever men make ready for outer space... from barren arctic IGY 
laboratories to tropical ballistic missile blast-offs...Trona chemicals are 
— ‘ putting the punch into high energy fuels and other propulsion systems. 
TRONA Is bas IC As the only basic producer of these three important raw materials used in 
high energy fuels and as solid propellant oxidizers—BORON, LITHIUM 


of hi g h ene rgy fuels Corporation has a vital stake in the space age. Recipient of the 14th ee 


research and technical development of these vital chemicals and their 


Boot American Potash & Chemical Corporation 


Plants: LOS ANGELES AND TRONA, CALIFORNIA; HENDERSON, NEVADA; SAN ANTONIO, TEXAS 
(American Lithium Chemicals, Inc. & San Antonio Chemicals, Inc.); WEST HANOVER, MASSACHUSETTS 


Recipient of 14th Annual Award for Chemical Engineering Achievement in BORON and LITHIUM 


Other Trona chemicals for high energy propulsion: POTASSIUM PERCHLORATE ® LITHIUM NITRATE ® LITHIUM PERCHLORATE ® LITHIUM METAL 
RUBIDIUM & CESIUM ® TRIMETHYL BORATE ® TRIMETHOXYBOROXINE © DECABORANE ¢ ELEMENTAL BORON ® BORON TRICHLORIDE ® BORON TRIBROMIDE 
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At the level of investigation in progress at 
the Flight Propulsion Laboratory, so many 
technical disciplines merge and overlap, that 
engineering management here can and will 


_ RE-DESIGN JOB SPECIFICATIONS 


TO UTILIZE YOUR TECHNICAL 


BACKGROUND TO THE FULL 
...where your experience is substantial. 


Opportunities now exist for highly quali- 
fied people to work with the Laboratory on 
advanced concepts and engineering princi- 
ples applicable to jet and rocket propulsion 
systems that are planned to operate up to 
Mach 5 and beyond. FPL is interested in you 
if you have had 8 to 5 years or more in the 
following areas: 


_ ADVANCED APPLIED MECHANICS « 
COMBUSTION e MECHANICAL DESIGN 
& ANALYSIS NUMERICAL ANALYSIS 


INTERESTED IN a os LURGICAL THEORY & PROCESSING 
e AERO-THERMO ANALYSIS 


_ (Also a number of positions for engineers and 
_ scientists with 1 to 3 years’ experience.) 


Write Mr. Mark Peters, 7. 
Building 100 


Department 34-MP 
Flight Propulsion Lab. 


PROPUL 
LABORATORY. 


GENERAL @@ ELECTRIC 


CINCINNATI 15, OHIO 


“Thrust and Progress’’"— 

showing creative de- 

velopment et FPL—will 
sent you On request. 


Female Coupling 
& Dust Cap 


Male Coupling 
& Dust Plug 


ho}man LIQUID OXYGEN 


QUICK COUPLINGS ey 
POSITIVE SEAL—NO GASKETS 


Hofman Quick Couplings are designed to utilize mechanical CAM ACTION giving 
metal to metal medium pressure contact in place of gaskets for a POSITIVE seal 
while requiring only % turn to open or close. Fabricated from high quality 
bronze, these units will give long, trouble-free life. Hofman Quick Couplings 
are available in 1”, 142”, 2”, 242”, 3” and 4” sizes. Armored transfer hoses 
with couplings installed can also be furnished. 


Send for our new 16 page catalog LOW TEMPERATURE APPARATUS — 


Laboratories, Inc., 
Dept. F, 5 Evans Terminal, Hillside, N. J. 


BLAIR-MARTIN CO., INC. 1010 Fair Oaks, So. Pasadena, California (Calif., Ariz., N. M.) 


WESTERN 
THE DARLING CO., P. O. Box 277, Wheatridge, Colorado (Colorado, Wyoming, Utah, Nevada) 
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urement of Ionization in Flames by a 


Conductivity Method, by J. Poncelet, 
Lowain Univ., Nov. 1956, 25 pp. (Con- 
tract AF 61(514)-922). 

Spark Ignition of Flowing Gases, hy 
Clyde C. Swett, Jr. NACA Rep. 1287, 
1956, 18 pp. 

Method for Combustion of Metals, hy 
A. V. Grosse, U. S. Patent. No. 2,764,109, 
Sept. 25, 1956, 9 pp. 

Starting Characteristics and Combus- 
tion Performance of Magnesium Slurry 
in 6/5-inch-diameter Ram-jet Engine 
Mounted in Connected-pipe Facility, 
by James B. Gibbs, NACA RM E53K05, 
Jan. 1954, 25 pp. (Declassified from Co: - 
fidential by authority of NACA Res. Als. 
116, p. 16, 7/5/57.) 

Performance of Slurries of 50 per cen 
Boron in JP-4 Fuel in 5-inch Ram-j 
Burner, by Thaine W. Reynolds and Dona 
P. Haas, NACA RM E54D07, June 
1954, 31 pp. (Declassified from Con- 
fidential by authority of NACA Rk: 
Abs. 116, p. 17, 7/5/57.) 

Combustion of Elemental Boron, Qua»- 
terly Summary Report for March-M: 
1957, Experiment Inc., TM-952, Jw 
1957, 11 pp. 


Blow-out Velocities of Various Petroleum 
Slurry and Hydride Fuels in a 1 7/8 Inch 
Diameter Combustor, by Preston N. Cook, 
Jr., Albert M. Lord, and Samuel Kay: 
NACA RM ¥E54A28, April 1954, 20 p; 
(Declassified from Confidential by autho: 
ity of NACA Res. Abs. 116, p. 17, 7/5/57 

Stabilization of 50-percent Magnesium- 
JP-4 Slurries with Some Aluminum Soaps 
of C, Acids, by Robert M. Caves, NACA 
RM E54C10, May 1954, 48 pp. (Declas- 
sified from Confidential by authority oi 
NACA Res. Abs. 116, p. 17, 7/5/57.) 

Preparation and Properties of Concen- 
trated Boron-Hydrocarbon Slurry Fuels, 
by Irving A. Goodman and Virginia O. 
Fenn, NACA RM E54F 18a, Aug. 1954, 
35 pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 116, 
p. 18, 7/5/57.) 

Effect of Surface-Active Additives on 
Physical Behavior of 50-Percent Slurries 
of 1.5-Micron Magnesium in n-Decane, 
by Murray L. Pinns, NACA RM E54K22a, 
Feb. 1955, 54 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs 
116, p. 18, 7/5/57.) 

Blow-Out Velocities of Several Slurry 

and Liquid Fuels in a 1 7/8 Inch Diameter 
Combustor, by James F. Morris, Robert 
M. Caves, and Albert M. Lord, NACA 
RM E54L27a, Feb. 1955, 12 pp. (De- 
classified from Confidential by authorits 
of NACA Res. Abs. 116, p. 18, 7/5/57.) 
“ Performance of Pure Fuels in Single 
J33 Combustors. II. Hydrocarbon and 
Nonhydrocarbon Fuels, by Arthur L. 
Smith and Jerrold D. Wear, NACA RM 
E55B02, 63 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
116, p. 19, 7/5/57.) 

NACA Research on Slurry Fuels 
Through 1954, by Walter T. Olson and 
Roland Breitweiser, NACA RM E55B14, 
April 1955, 30 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abs. 116, p. 19, 7/5/57.) 

Preparation and Handling of Magne- 
sium-Hydrocarbon Slurries for Jet-Engine 
Applications, by Henry C. Barnett, A. M. 
Lord, and P. H. Wise, NACA RM E55D01, 
June 1955, 60 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
116, p. 19, 7/5/57.) 

Preparation of 50 Percent Boron- 
Hydrocarbon Slurries Using Combinations 
of Glycerol Sorbitan Laurate With Var- 
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Seventy-live Years... 


Ten MILES an hour was “speed” in 1883 when produced by the most advanced techniques of the 
Wyman-Gordon started to make forgings for the high- day. 

wheel bicycle. Through the 75 intervening years Wyman-Gordon is proud of its achievements in _ 
forgings have made important contributions to the these fields and, as the largest producer of automotive _ 
phenomenal advances in propulsion. Progress from and aircraft forgings, is prepared to accept the chal- 
the first “horseless carriages”. ..from the early “‘fly- lenge of the future. Today, as for 75 years, there 
ing machines”. ..to the supersonic speeds of today is no substitute for Wyman-Gordon quality and 
...would not have been possible without forgings experience. 


WYMAN-GoRDON COMPANY 


Established 1883 


i FORGINGS OF ALUMINUM @ MAGNESIUM @ STEEL @© TITANIUM 
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ARNOUX 
LOX TEMPERATURE PROBES 


Semiconductor Series 7X XXL provide output in 
volts* for a temperature span of only 20 F in the 
liquid oxygen (LOX) temperature range...also 
can be used with liquid nitrogen, liquid argon, 
and almost any liquified gas from —240 F to 
—320 F range. 

Applications include temperature-stratification 
measurements (temperature gradients through a 
LOX tank), sensing liquid levels in tanks, and 
pipeline measurements during LOX-fueling 
operations. 


*From 0 to 5 vde, used with Arnoux TME System a 
or similar equipment. 


Nominal resistance: 12,000 ohms at —-77 F 

Calibration accuracy: 0.1 F 
Repeatability and hysteresis: within 

calibration accurac 

Calibration points: —240 F and —320 F; 1 

increments available between —275 © and 


ARNOUX foremost in TEMPERATURE MEASUREMENT EQUIPMENT... 
TRANSDUCERS, ACCESSORIES AND CUSTOM ENGINEERED SYSTEMS 


PORATION 


11924 WEST WASHINGTON BLV . AN 


YEARS OF EXPERIENCE 
in the manufacture of 

| SMALL MACHINE PARTS 

amt, OF HIGHEST PRECISION 


The equipment you design and produce demands the finest quality in each 


of its components. For those machine parts requiring fine finishes, difficult 
machining operations, and closest tolerances specify LaVezzi and be certain 
of parts that meet approval. Your quotation requests will receive prompt 


: and careful consideration. Illustrated brochure will be mailed upon request. 


WEST LAKE ST., CHICAGO, I 


ious Thickeners, by Irving A. Goodman 
and Virginia O. Fenn, NACA RM E55E17, 
July 1955, 30 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
116, p. 19, 7/5/57.) 

Effect of Temperature on Viscosity of 
Slurries of Boron and Magnesium in 
JP-5 Fuel, by Harold F. Hipsher, NAC A 
RM E55E19, July 1955, 20 pp. (Deelss- 
sified from Confidential by authority of 
NACA Res. Abs. 116, p. 19, 7/5/57.) 

Some Dielectric Properties of Suspen- 
sions of Boron Powders in Mineral Oil, 
by Aubrey P. Altshuller, NACA R 
E55F02, Aug. 1955, 21 pp. (Declassificd 
from Confidential by authority of NAC 1 
Res. Abs. 116, p. 19, 7/5/57.) 


Materials of Construction 


The Physical Properties of Several M2- 
terials for Use in Photo-Thermelastic 
Investigations, by Herbert Trampos:h 
and George Gerard, N.Y. Univ., Coi!. 
Engng., Res. Div., Tech. Rep. SM 57-5 
(AFOSR-TN-57-282), May 1957, 27 pp. 

Internal Friction and Shear Modulus of 
Thorium at High Temperatures, by C. |. 
Dixon and H. Hori, Atomic Energy, 
Comm., NAA-SR-1846, Apr. 1957, 12 pp 


The Anode Oxidation of Some Dilute 
Binary Zirconium Alloys, by George |: 
Adams Jr., Curtis E. Borchers and Pierre 
Van Rysselberghe, Atomic Energy Comm., 
AECU-3388, Jan. 1957, 26 pp. 

Properties of Titanium Compounds and 
Related Substances, by Frederick D. 
Rossini, Phyllis A. Cowie, Frank ©. 
Ellison, Clarence C. Browne and William 
C. Arsen, Office Naval Res., Rep. ACR-17, 
Oct. 1956, 448 pp. 

Creep Behavior of Circular Plates, by Lb. 
Venkatraman and P. G. Hodge Jr., 
Brooklyn Polytech. Inst., Dept. Aeron. 
Engng. Appl. Mech., Rep. 369 (AFOSR- 
TN-57-135; ASTIA AD 120491), Apr. 
1957, 24 pp. 

Calculated and Measured Stresses in 
Simple Panels Subject to Intense Random 
Acoustic Loading Including the Near Noise 
Field of a Turbojet Engine, by Leslie W. 
Lassiter and Robert W. Hess, NACA TN 
4076, Sept. 1957, 33 pp. 

Practical Solution of Plastic Deforma- 
tion Problems in Elastic-Plastic Range, b\ 
A. Mendelson and S. 8. Manson, NACA 
TN 4088, Sept. 1957, 52 pp. 

Special Welding Techniques; Final 
Summary Report for 1954, Glenn L 
Martin Co., Atomic Energy Comm., AE- 
CU-3174, Jan. 1955, 41 pp. 

Corrosion Properties of Various Ma- 
terials in High Temperature Waters, b\ 
C. J. Lancaster and W. L. Williams, Naval 
Engng. Expt. Sta., Rep. 4A(23)966870 
(Office Tech. Services, PB 111963), Jan. 
1953, 13 pp. 

Report on Performance of Materials 
Tested in Water at High Temperature, 
by C. J. Lancaster, Naval Engng. Expt 
Sta., Rep. 4A(16)966870 (Office Tech. 
Services, PB 111962), March 1952, 6 pp. 

Anodized Aluminum Wire Takes 1000 
Deg F, Aviation Age, vol. 28, Sept. 1957, 
pp. 64-67. 

Re-usable Hermetic Seal Takes High 
Temperatures, Aviation Age, vol. 28, 
Sept. 1957, pp. 88-91. 

‘Exotic’? Metals Look Good for To- 
morrow’s High Temperature Structures, 
by Irwin Stambler, Aviation Age, vol. 28, 
Sept. 1957, pp. 42-47. 

‘Aluminum Alloys, by R. L. Horst, 
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Indust. Engng. Chem., vol. 49, Sept. 1957, 
Pt. 2, pp. 1578-1583. 

Carbon and Graphite, by W. M. Gay- 
lord, Indust. Engng. Chem., vol. 49, Sept. 
1957, Pt. 2, pp. 1584-1587. 

Chat by John W. Lennon, /ndust. 

ng. Chem., vol. 49, Sept. 1957, Pt. 2, 
1592. 

Less Common Metals, by E. M. Sher- 
wood, Indust. Engng. Chem., vol. 49, Sept. 
1957, Pt. 2, pp. 1612-1617. 

Nickel Including High-Nickel Alloys, 
by R. M. Fuller, Indust. Engng. Chem., 
vol. 49, Sept. 1957, Pt. 2, pp. 1618-1628. 

Protective Coatings, by Francis Sco- 
lield, Indust. Engng. Chem., vol. 49, Sept. 
957, Pt. 2, pp. 1639-1642. 

Stainless Steels Including Other Fer- 
rous Alloys, by Walter A. Luce, Indust. 
Engng. Chem., vol. 49, Sept. 1957, Pt. 2 
pp. 1643-1652. 

Titanium, by H. B. Bomberger, /ndust. 
Engng. Chem., vol. 49, Sept. 1957, Pt. 2, 
pp. 1658-1662. 

Non-metallics for Missiles, by John H. 
Lux and Robert L. Noland, Missiles and 
Rockets, vol. 2, Sept. 1957, pp. 137-138, 
(41-142, 144. 

Cooperative Investigation of Relation- 
ship between Static and Fatigue Proper- 
ties of Wrought N-155 Alloy at Elevated 
Temperatures, NACA Subcommittee on 
Power Plant Materials, NACA _ Rep. 
1288, 1956, 35 pp. (supersedes NACA T'N 
3216). 

Generalized Master Curves for Creep 
and Rupture, by George J. Heimer] and 
Arthur J. McEvily Jr., NACA TN 4112, 
Oct. 1957, 31 pp. 

High-Temperature Oxidation and Ig- 
nition of Metals, by Paul R. Hill, David 
Adamson, Douglas H. Foland and Wal- 
ter E. Bressette, NACA RM 1L55L23b, 
March 1956, 12 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstracts 120, p. 8, 10/3/57.) 

Investigation of a Cermet Gas-turbine- 
blade Material of Titanium Carbide In- 
filtrated with Hastalloy C, by Charles A. 
Hoffman, NACA RM E55H12, Jan. 1956, 
23pp. (Declassified from Confidential by 
authority of NACA Res. Abstracts 120, p. 
8, 10/3/57.) 

Investigation of Some Mechanical Prop- 
erties of Thermenol Compressor Blades, 
by Donald F. Johnson, NACA TN 4097, 
Oct. 1957, 14 pp. 

Experimental Investigation of Cermet 
Turbine Blades in an Axial-Flow Turbojet 
Engine, by William C. Morgan and 
George C. Deutsch, NACA TN 4030, Oct. 
1957, 20 pp. 

Theoretical Study of Decomposable 

Varnishes for Combustion Chamber Wall 
Protection, by H. Gelly, Fusées et Re- 
cherche Aeron., vol. 2, no. 13, July 1957, 
pp. 219-228 (in French). 
Thermal Behavior of Titanium and Its 
Alloys, by R. Syre, Fusées et Recherche 
Aeron., vol. 2, no. 3, July 1957, pp. 229- 
235 (in French). 

Uses of Glass Cloth Reinforced Plastics 
in Aviation, Missiles and Rockets, by J. 
Duflos, Fusées et Recherche Aeron., vol. 2, 
no. 3, July 1957, pp. 237-240 (in French). 


Instrumentation, Data 
Recording, Telemetering 


Time Constants and Frequency Re- 
sponse of Coated Hot Wires as Turbu- 
lence Sensing Elements, by Avis Borden, 
David Taylor Model Basin, Rep. 952, June 
1957, 48 pp. 
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ROCKETDYNE F2 


quiz for 


Rocket Propulsion Engineers 


(Time Limit: 10 years) 


1. What is the maximum theoretical specific impulse obtainable 
for a rocket motor with zero combustion volume? 

2. What is the basic mechanism involved in jet separation in 
nozzles, and how can the separation point be predicted as a func- 
tion of nozzle divergence angle? 

3. What is the effect on turbine blade design of a working gas 
which is still reacting and changing its composition? 

4. What are the dynamic interrelations among the various com- 
ponents of a rocket engine, and what control systems will best 
solve the rocket engine control problems? 

5. What is the flame holding mechanism in a rocket motor and 
how can fundamental understanding of such a mechanism be 
used to increase reaction ratio? 

6. What types of fluid transport systems offer simplicity, effi- 
ciency, and light weight (other than the conventional 
turbopump)? 

7. What are some workable methods of translating chemical 
energy into thrust — other than by the standard combustion at 
high pressure and then conversion to velocity? 

If you would like to work in these areas and have an outstand- 
ing technical background, an inquiring mind, and the creative 
ability to visualize new solutions, possess a M.S. or Ph.D. degree 
in Physics, Chemistry, M.E., A.E., Ch.E., E.E., or a B.S. degree 
in the above with a thorough background in applied research, we 
have a place for you at Rocketdyne. 

Please write to Mr. A. D. Jamieson, Engineering Personnel, 
6633 Canoga Avenue, Canoga Park, California. 


A DIVISION OF NORTH AMERICAN AVIATION, INC 


BUILDERS OF POWER FOR OUTER SPACE 
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Pressure Vessel 


Technology 


and now Titant 


Spheres of titanium alloy represent 
the optimum configuration and 
material of maximum strength/weight 
for the storage of gases and liquids 
under extremes of pressure and 
temperature. 

Rheem Aircraft, presently supplying 


ai 

7 pressure units to current ICBM pro- 
af grams, is directing its years of experi- 
a ence to the development and pro- 
duction of all types of vessels and 
propellant chambers. 

RHEEM MANUFACTURING CO. 
5 AIRCRAFT DIVISION 

7 11711 woodruff avenue, downey, california 


Engineers: Join Rheem in challenging technical work. 


Project Vanguard Report 18. Minitrack. 
Report No. 1-Phase Measurement, by C. 
A. Schroeder, C. H. Looney Jr. and H. E. 
Carpenter Jr., Naval Res. Lab. Rep. 
4995, July 1957, 28pp. 

Programming and Coding for ORDVAC, 
by Tadeusz Leser and Michael Romanelli, 
Aberdeen Proving Ground, Ballistic Res. 
Labs., Rep. 997, Oct. 1956, 278 pp. 

What You Should Know for Reliable 
Missile Circuitry, by D. D. Zimmerman, 
Aviation Age, vol. 28, Sept. 1957, pp. 50-55. 

A Digital System for Collecting and 
Processing Flight Test Data, by H. W. 
Royce, SAE Preprint 192, Oct. 1957, 3 pp. 

How Badly Are High Temperature Elec- 
tronic Systems Really Needed? by Rob- 
ert R. Janssen, SAE Preprint 230, Oct. 
1957, 6 pp. 

North Atlantic Treaty Organization, 
Advisory Group for Aeronautical Research 
and Development, Flight Test Manual, 
vol. 4: Instrumentation Systems, 1957. 

Hot Wire Apparatus for Measurement 
of Transpiration Effects on Heat Transfer 
Coefficients, by John M. Lenoir, Calif. 
Inst. Tech., Jet Propulsion Lab., Progress 
Rep. 20-300, Nov. 1956, 30 pp. 


Guidance Systems and 
Components 


Systems Analysis Approach to Firebee 
Flight Control Systems Development, by 
C.D. Corden and H. J. Hansen Jr., SAE 
Prepr. 230, Oct. 1957, 8 pp. 

Aids to Inertial Navigation, by Fred- 
erick Stevens and Frank W. Lynch, Aeron. 
Engng. Rev., vol. 16, no. 11, Nov. 1957, pp. 
43-47. 

_ Missile-Tracking Camera, by George L. 
Christian, Aviation Week, vol. 67, no. 17, 
Oct. 28, 1957, pp. 110-111, 114. 

Radius of Action of a Guided Missile; 
Its Dependence Upon Certain Arbitrary 
Parameters, by H. R. Voellmy, Fusees et 


Recherche Aeron., vol. 2, no. 3, July 1957, 


pp. 213-218 (in French). 


Servomechanisms 
and Controls 


Requirements, Parameters and Design 
Considerations for Pneumatic Inlet Con- 
trol Systems, by R.W. Matzdorff and C.F. 
ee SAE Preprint. 236, Oct. 1957, 

lpp. 

Packaged Hydraulic Systems for Mis- 
siles, by A. L. Stone and J. W. Woodward, 
SAE Prepr. 214, Oct. 1957, 6 pp. 

Why Hate High Pressure Pneumatics? 
by H. E. Wright, SAE Prepr. 237, Oct. 
1957, 7 pp. 

Hydraulic Pumps and Motors for High 
Velocities, by C. R. Himmler, Fusées et 
Recherche Aeron., vol. 2, no. 3, July 1957, 
pp. 248-250 (in French). 


Flight Vehicle Design 
and Testing 


Proceedings, Symposium on Environ- 
mental Testing of Aircraft and Missiles, 
SAE, Special Publication SP141, Dec. 5-7, 
1955. 

Guided Missile Sensation at Farnbor- 
ough, by Miles, [nteravia, vol. 12, no. 10, 
Oct. 1957, p. 1004. 

Helium Comes into Its Own in Missile 


| Technology, by Kurt R. Stehling, Avia- 
| tion Age, vol. 28, Oct. 1957, pp. 86-92. 
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Jupiter-Thupiter-Thorpiter-Thor? Mis- 


siles and Rockets, vol. 11, no. 10, Oct. 1957, 
p. 45. 

Missiles and Rockets Around the Globe, 
Missiles and Rockets, vol. 11, no. 10, Oci. 
1957, pp. 71-106. 

Missiles for NATO, by Normal |., 
Baker, Missiles and Rockets, vol. 11, no. 
10, Oct. 1957, pp. 107. 

French Missile Production, by Jean- 
Marie Riche, Missiles and Rockets, voi. 
11, no. 10, Oct. 1957, p. 112. 

Air-to-Air Rockets, by F. I. Ordway 
III and R. C. Wakeford, Ordnance, vel 
38, no. 224, Sept.-Oct. 1957, p. 239. 

Guided Missiles of the U. S. Armed 
Forces, Ordnance, vol. 38, no. 224, Sept.- 
Oct. 1957, p. 264. 

British Guided Missiles, by James Hay 
Stevens, Aeronautics, vol. 37, Sept. 1957, 
pp. : 
Maximum Ranges of Intercontinentz! 
Missiles, by Derek F. Lawden, Aerov. 
Quarterly, vol. 8, Aug. 1957, pp. 269-27% 

Guided Weapon Engineering, by Brig::- 
dier J. Clemow, Aircraft Engng., vol. 2°, 
Sept. 1957, pp. 258-267. 

French Test Ramjet-Powered Missiles, 
by Robert Roux, Aviation Age, vol. 2%, 
Sept. 1957, pp. 62-65. 

> 
Space Flight 

P 

The Simulation of Gravity, by Derek F. 
Lawden, J. Brit. Interplanet. Soc., vol. 16, 
no. 3, July-Sept. 1957, pp. 134-140. 

On the Equation of Rocket Motion, b» 
G. Leitmann, J. Brit. Interplanet, Soc 
vol. 16, no. 3, July-Sept. 1957, pp. 141-147 

The Use of Probe Rockets, by C. A. 
Cross, J. Brit. Interplanet. Soc., vol. 16, no 
3, July-Sept. 1957, pp. 148-162. 

Moon Cities Planned by Red Civil De- 
fense Chief, Missiles and Rockets, vol. 11, 
no. 10, Oct. 1957, p. 45. 

Project Far Side, Missiles and Rockets, 
vol. 11, no. 10, Oct. 1957, pp. 120-128. 

Sectional Satellites, by George W. 
Hoover, Missiles and Rockets, vol. 11, 
no. 10, Oct. 1957, pp. 135-137. 

Summing Up the Satellite, The Aero- 
plane, vol. 93, no. 2407, Oct. 18, 1957, pp. 
581-582. 

Bridgehead in Space, by O. C. Winzen, 
Interavia, vol. 12, no. 10, Oct. 1957, pp. 
1040-1041. 

Orbital Behavior of Earth Satellites, by 
Robert E. Roberson, J. Franklin Inst., 
vol. 264, no. 3, Sept. 1957, pp. 181-202. 

Orbital Behavior of Earth Satellites, II, 
by Robert E. Roberson, J. Franklin Inst., 
vol. 264, Oct. 1957, pp. 269-285. 

Spiral and Elliptical Orbits, by E. F. 
Relf, J. Royal Aeron. Soc., vol. 61, Sept. 
1957, pp. 635-641. 

Air Force Research and Development in 
Space Technology, by Thomas A. Power, 
Aeron. Engng. Rev., vol. 16, Sept. 1957, pp. 
36-39. 

Evidence for Vegetation 
on Mars, by William M. Sinton, Astrophys. 
J., vol. 126, Sept. 1957, pp. 231-239. 


Astrophysics, Aerophysics 


Gas-Dust Interaction in the Accretion of 
Interstellar Dust by the Sun, by Charles 
B. Stephenson, Astrophys. J., vol. 126, 
July 1957, pp. 195-201. 

Cosmic Ray Nuclear Interactions in Ni- 
trogen Gas, by Saburo Miyake, Kensaku 
Hinotani, Itsuo Katsumata and Tatsuno- 
suke Kaneko, J. Physical Soc. Japan, vol. 


12, Aug. 1957, pp. 845-854. 
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Nuclear Propulsion 


An Experiment on the Structure of the 
Electron, I, by G. K. O'Neill and J. A. 
Ball, Atomic Energy Comm., NYO 8015, 
July 1957, 15pp. 

A Series of Thermal, Epithermal and 
Fast Neutron Measurements in the MTR, 
by J. B. Trice, Atomic Energy Comm., CF- 
55-10-140, Oct. 1955, 88pp. 


A Review of the Reaction Kinetics of 
Deuterium and Tritium Compounds, VII: 
Oxidation Reduction Reactions, by Law- 
rence M. Brown, Atomic Energy Comm., 
N BS-4877, Nov. 1956, 25pp. 

Unclassified U. S. Atomic Energy Com- 
mission Reports on Controlled Thermo- 
nuclear Processes; A Literature Search, 
by Sidney F. Lanier, Atomic Energy Comm., 
TID 3506, Jan. 1957, 7pp. 

Preliminary Results on the Kinetic Be- 
havior of Water Boiler Reactors, by D. L. 
Hetrick, J. W. Flora, E. L. Gardner, B. R. 
Moskowitz, D. R. Muller, M. E. Rem- 
ley, and R. K. Stitt, Atomic Energy Comm., 
NAA-SR-1896, Apr. 1957, 20pp. 

A General Method for Comparing Ther- 
mal Performance of Fuel Element Geom- 
etries and Coolants for Non-boiling 
Reactors, by G. A. Freernd and A. L. 
London, Atomic Energy Comm., ANL- 
5589, Jan. 1957, 30 pp. 

Unclassified U. S. Atomic Energy Com- 
mission Reports on Reactor Safety, A Lit- 
erature Search, by Hugh E. Voress, 
Atomic Energy Comm., TID-3503, Nov. 
1956, 9 pp. 

Symposium on Reactor Instrumentation, 
Part I, Alomic Energy Comm., WAPD 20, 
June 1950, 46pp. 

Coulomb Effects in Inner Bremsstrah- 
len, by R. R. Lewis, Jr., and G. W. Ford, 
Physical Rev., vol. 107, Aug. 1, 1957, pp. 
756-765. 

Nuclear-powered Seaplane, by A. D. 
Struble, Jr., SAE J., vol. 65, Sept. 1957, 
pp. 64-66. 

Thermal Conductivity of Insulating 
Materials for Use in Nuclear Reactors, by 
W. R. Morgan and W. G. Baxter, ASME 
Paper 57-F-9, Sept. 1957, 5 pp. 

Direct Interaction Theory of Inelastic 
Scattering, III: Numerical Calculations, 
by Carl A. Levinson and Manoj K. Baner- 
jee, Atomic Energy Comm., NYO-8012, 
Sept. 1957, 22 pp. 

Solid State Physics; Semiannual Prog- 
ress Rep. July-Dec. 1956, North Ameri- 
can Aviation, Inc., Atomic International, 
Atomic Energy Comm., NAA-SR-1885, 
June 1957, 40 pp. 


New Concepts and Materials Needed 
for Nuclear Accessory Systems, by John 
P. Wittry, Aviation Age, vol. 28, Oct. 
1957, pp. 32-35. 

Chemical Effects of Hot Atoms, by 
Louis Monchick, Koichi Funabashi and 
John L. Magee, J. Chem. Phys., vol. 27, 
Sept. 1957, pp. 734-739. 

The Application of Nuclear Power to 
Logistic Aircraft Systems, by Robert W. 
Middlewood and Robert B. Ormsby, Jr., 
SAE Preprint 206, Oct. 1957, 11 pp. 


Economics of the Reactor Fuel Cycle, by 
Harlan W. Nelson, Battelle Tech. Rev., 
vol. 6, Oct. 1957, pp. 9-16. 


Nuclear Planes Break Range Barrier, 
by R. W. Middlewood and R. B. Ormsby 
Jr., SAE J., vol. 65, no. 11, Oct. 1957, pp. 
32-34. 


Power from Fusion Reactions, by Clay- 
ton R. Lewis, SAE J., vol. 65, no. 11, 
Oct. 1957, pp. 48-49. 
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Here’s how Hiller solved 
it with Du Pont Blind 
Expansion Rivets 


The Army’s H-23D Helicopter was de- 
signed to boast a more powerful (250 
h.p.) engine and have longer range 
than its predecessor, the H-23C. 

Consequently, a larger fuel cell had 
to fit in the space under the engine 
deck. No problem there. But then the 
deck had to be secured in place—and 
there was no space left underneath for 
bucking the rivets. 

DuPont Blind Expansion Rivets 
came to the rescue, by permitting Hiller 
Helicopters to do this assembly with- 
out bucking. Everything’s done topside 
and in the open. This convenience, plus 
speed of installation at minimum cost, 
makes these rivets the ideal solution to 
a wide variety of fastening problems. 


Fast Setting—One man, working from 
the head of the rivet, can easily set 
20-25 rivets per minute. 


No Refinishing—There’s no need for 
cutting, trimming, or buffing after the 
rivets are set. 


Liner of fuel cell prevented access for bucking 
rivets under engine deck of H-23D. DuPont 
Blind Expansion Rivets, set by applying heat to 
head of each rivet, came to the rescue. 


Positive Holding Power—The rivet shank | 
expands to fill the hole completely, in- __ 
suring tight joints. : 
Find out how Du Pont Aircraft Blind — 
Expansion Rivets can solve your fas-— 
tening problems. Call your DuPont | 
representative or write: E. I. duPont — 
deNemours & Co. (Inc.), Explosives — 
Department, Wilmington 98, Delaware. 


GU POND 


REG. Pat. OFF 


A Product of Du Pont Research 


DU PONT AIRCRAFT BLIND EXPANSION RIVETS P 


BETTER THINGS FOR BETTER LIVING... 


THROUGH CHEMISTRY 


SEE THE “DU PONT SHOW OF THE MO 


NTH” ON CBS 
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PROPULSION SYSTEM 


PROJECT ENGINEERS 


Advancements in ballistic mis- 
sile and space vehicle design 
poseinteresting problemsinthe 
integration of propulsion sub- 
systems to obtain satisfactory 
performance and reliability. 
Technical coordination with 
rocket engine manufacturers in 
the development of liquid fuel 
propulsion systems presents 
significant opportunities for 
the propulsion engineer. Space 
Technology Laboratories has 
positions in this field for en- 
gineers with five to ten years 
of experience in liquid fuel 
propulsion, preferably rocket 
propulsion. 


Inquiries regarding these 
positions are invited. 
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